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APPLIED  TECHNOLOGY  LABORATORY  POSITION  STATEMENT 

This  report  presents  a method  for  determining  aerodynamic  character- 
istics of  helicopter  shapes  under  the  influence  of  the  main  rotor 
wake.  The  method  is  condidered  to  be  worthy  of  publication  for 
dissemination  of  information  and  the  stimulation  of  further  related 
research.  The  reader  is  cautioned  that  this  method  does  not  predict 
flow  separation  as  the  title  would  imply,  nor  does  the  rotor  wake 
fully  impinge  upon  the  body  of  the  helicopter.  The  method  is  useful,  , 

however,  as  a design  tool  in  determining  rotor-fuselage  aerodynamic 
interference.  * 

< 
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engineer  for  this  effort. 
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An  investigation  has  been  conducted  to  demonstrate  the  use  of  the  Green's  function 
method  to  study  rotor  wake  effects  on  helicopter  hub/pylon  flow.  This  report  consists  of 
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puter program  SHAPES  (Subsonic  Helicopter  Aerodynamics  Program  with  Effects  of  ^ 
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Volume  I presents  the  theoretical  formulation  and  corresponding  numerical  procedure 
for  the  study  of  incompressible  potential  aerodynamics  with  separated  flow.  While  the 
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formulation  is  valid  for  fully  unsteady  aerodynamics,  this  report  is  mainly  concerned  with 
rotor  aerodynamics.  A potential  flow  aerodynamic  program,  SHAPES,  with  suitable  rotor 
wake  representation  was  developed  to  predict  the  separation  characteristics  of  arbitrary 
three-dimensional  helicopter  configurations.  In  particular,  the  effect  of  the  rotor  blade 
wake,  blade  shank  wake,  and  hub  wake  in  the  separation  of  the  flow  over  a lifting  heli- 
copter in  forward  flight  is  analyzed. 

The  present  method  has  potential  application  in  the  design  of  helicopters  because  it 
provides  an  analytical  capability  which  can  be  used  to  develop  low -drag  profile  as  well 
as  to  explore  problem  areas.  Extensive  numerical  results  obtained  from  Program  SHAPES 
demonstrated  the  flexibility  and  accuracy  of  the  method.  These  results  are  in  excellent 
„ agreement  with  existing  data. 

"^Volume  II,  The  User's  Manual,  describes  the  structure  and  use  of  Program  SHAPES. 
SHAPES  is  written  in  FORTRAN  IV  for  operation  on  the  CDC  6600  digital  computer  system 
at  the  Applied  Technology  Laboratory,  U.S.  Army  Research  and  Technology  Laboratories 
(AVRADCOM),  Fort  Eustis,  VA.  The  User's  Manual  (Volume  II)  contains  detailed  informa- 
tion for  preparing  an  input  data  deck  and  interpreting  the  computed  data;  a discussion  of 
various  subroutines,  flow  charts,  common  storage  and  definition  of  FORTRAN  variables; 
sample  cases  to  illustrate  the  program  output;  and  a FORTRAN  listing  of  the  program. 
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This  report,  prepared  by  Aerospace  Systems,  Inc.  (ASI),  Burlington, 
Massachusetts,  for  the  U.S.  Army  under  Contract  DAAJ02-75-C-0041 , documents 
the  results  of  research  performed  during  the  period  June  1975  to  September  1977. 

The  study  was  sponsored  by  the  Eustis  Directorate,  U.S.  Army  Air  Mobility  Research 
and  Development  Laboratory,*  Fort  Eustis,  Virginia,  with  Mr.  Frederick  A.  Raitch 
serving  as  Technical  Monitor. 

The  effort  was  directed  by  Mr.  John  Zvara,  President  and  Technical 
Director  of  ASI.  Mr.  Paul  Soohoo  served  as  Principal  Technical  Staff  Member  under 
the  supervision  of  Mr.  Richard  B.  Noil,  Vice  President  of  ASI . Dr.  Luigi  Morino, 
Director  of  Computational  Continuum-Mechanics  Program  at  Boston  University,  and 
Dr.  Norman  D.  Ham,  Director  of  the  V/STOL  Technology  Laboratory  at  MIT,  con- 
tributed to  the  study  as  Principal  Consultants. 


*Redesignated  Applied  Technology  Laboratory,  U.S.  Army  Research  and 
Technology  Laboratories  (AVRADCOM),  effective  September  1,  1977. 
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SECTION  1 
INTRODUCTION 

This  volume  describes  the  structure  and  use  of  a digital  computer  program, 
SHAPES,  Subsonic  Helicopter  Aerodynamic  Program  with  Effects  of  Separation,  which 
was  developed  by  Aerospace  Systems,  Inc.  (ASI),  to  analyze  the  incompressible 
potential  aerodynamics  of  helicopter  configuration  with  rotor  wakes  and  separation 
effects.  Volume  I of  this  report  contains  a complete  description  of  the  theoretical 
formulation  and  corresponding  numerical  procedure  for  the  aerodynamic  method  for 
use  in  the  computer  program.  It  also  includes  extensive  numerical  results  showing  the 
flexibility  and  accuracy  of  the  method  as  well  as  comparison  with  several  existing 
results.  Proper  use  of  SHAPES  is  predicated  on  the  assumption  that  the  user  is  familiar 
with  the  techniques  and  limitations  set  forth  in  Volume  I. 

SHAPES  is  written  almost  entirely  in  FORTRAN  IV  for  operation  on  the 
CDC  6600  digital  computer  system  at  the  Applied  Technology  Laboratory,  U.S.  Army 
Research  and  Technology  Laboratories  (AVRADCOM),  Fort  Eustis,  Virginia.  The  pro- 
gram was  developed  with  a highly  modular  structure  for  ease  of  program  checkout,  to 
simplify  the  user's  understanding  of  the  program,  and  to  facilitate  any  modifications 
that  might  be  required  for  future  applications. 

Sections  2 and  3 contain  programming  details  of  SHAPES:  functions  of 
the  various  routines,  flow  charts,  common  storage,  and  definition  of  FORTRAN 
variables.  The  use  of  the  program  is  presented  in  Section  4.  Sample  cases  are  pre- 
sented in  Section  5 to  illustrate  the  output  of  Program  SHAPES.  Section  6 contains 
the  FORTRAN  listing  of  the  program. 
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SECTION  2 

COMPUTER  PROGRAM  SHAPES 


Computer  program  SOUSSA  (Steady  Oscillatory  Unsteady  Subsonic  and 
Supersonic  Aerodynamics)  discussed  in  Subsection  2.2  of  Volume  I has  been  modified 
to  include  rotor  dynamics  and  separation  from  hub/pylon  components.  The  revised  pro- 
gram  is  called  SHAPES,  for  Subsonic  Helicopter  Aerodynamic  Program  with  Effects  of 
Separation.  Geometry  preprocessor  for  single-rotor  helicopter  configurations  is 
included  in  Program  SHAPES. 

In  Program  SHAPES  the  user  need  not  be  familiar  with  the  aerodynamic  por- 
tion of  the  program,  unlike  other  sophisticated  aerodynamic  programs  in  which  the 
choice  for  the  combination  of  various  aerodynamic  functions  (sources,  doublets,  vortices) 


is  an  art  which  requires  considerable  understanding  of  the  method.  Another  advantage 
of  SHAPES  is  that  the  paneling  used  for  the  aerodynamics  is  completely  arbitrary  and, 
therefore,  may  coincide  with  the  one  used  for  structural  analysis.  A typical  finite 
element  representation  of  a helicopter  fuselage  configuration  is  shown  in  Figure  1 . 

The  aerodynamic  paneling  can  be  accomplished  manually  by  inputting  the  coordinates 
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Figure  1 . Finite  Element  Representation  of  Arbitrary  Helicopter  Fuselage. 
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of  the  corner  points  of  all  selected  aerodynamic  panels  or  by  using  a geometry  pre- 
processor. A geometry  preprocessor  for  conventional  helicopter  configurations  (see 
Figure  2)  is  incorporated  in  SHAPES.  The  fuselage  is  assumed  to  have  an  elliptical 
center  section  with  the  nose  and  tail  approximated  by  elliptical  paraboloids.  An 
ellipsoidal  element  is  used  to  approximate  the  rotor  hub  and  pylon  sections.  In  addi- 
tion, the  rotor  shaft  and  blade  shank  are  represented  by  cylindrical  sections.  Given 
the  helicopter  configuration  geometry,  the  preprocessor  computes  the  corner  point 
locations  of  required  aerodynamic  panels  for  the  complete  helicopter  configuration. 
Required  inputs  include: 

• Shape  and  dimensions  of  fuselage,  pylon,  shaft,  and  hub 

• Number  of  blades 

• Rotor  diameter 

• Extent  of  root  cutout  or  shank 

• Cross-section  of  shank 

• Rotor  blade  airfoil  designation 

• Rotor  blade  chord  distribution 

• Rotor  blade  thickness  distribution 

• Rotor  blade  twist  distribution 

• Rotor  tip  speed 

• Freestream  velocity 

• Rotor  wake  geometry 

This  capability  eliminates  the  requirement  for  the  selection  of  the  aerodynamic  panel- 
ing by  the  user. 
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SECTION  3 

PROGRAM  SHAPES  DESCRIPTION 

This  section  contains  a brief  description  of  the  organization  of  Program 
SHAPES  and  its  subroutines,  a flow  diagram,  a table  of  common  block  storage,  and  a 
list  of  FORTRAN  variables. 

3.1  PROGRAM  SHAPES  ORGANIZATION 

The  present  computer  program  for  incompressible  potential  aerodynamics  with 
separated  flow  for  a complete  helicopter,  i.e.,  fuselage,  pylon,  and  rotor,  contains 
three  main  parts:  (1)  geometric  preprocessor,  (2)  coefficient  matrix,  and  (3)  pressure 
distribution  and  force.  Based  on  data  provided,  the  geometric  preprocessor  automatically 
generates  the  aerodynamic  panels  for  the  helicopter  as  well  as  the  numbering  and  loca- 
tion of  the  panel  corner  points  in  the  reference  Cartesian  (global)  coordinate  system. 

The  matrix  coefficients  are  determined  by  evaluating  the  doublet  and  source  integrals 
over  the  surface  of  the  aerodynamic  panels.  This  yields  N simultaneous  linear  equations 
with  N unknown  velocity  potentials  at  the  centroids  of  N elements.  A standard  IBM 
subroutine  GELG  (Gaussian  elimination  method)  is  used  to  solve  for  the  unknown 
potentials.  To  obtain  a continuous  distribution  for  the  velocity  potential,  an  averaging 
scheme  is  introduced.  Hence,  the  perturbation  velocities  and  the  pressure  at  the  cen- 
troid of  each  element  can  be  evaluated.  Finally,  the  aerodynamic  coefficients,  i.e., 
lift  and  induced  drag,  are  computed. 

3.1.1  GEOMETRIC  PREPROCESSOR 

Presently,  five  main  subroutines  form  the  geometric  preprocessor: 

• DATA 

• PREPRO 
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• COODPT 

• GEOMET 

• VEC123 

Basic  geometric  inputs  of  the  problem  are  defined  by  the  user  in  DATA.  This  subroutine 
automatically  generates  the  global  subsurface  numbering  of  the  complete  helicopter 
geometry;  i.e.,  fuselage,  pylon,  shaft,  hub,  shank,  blade,  and  vertical  tail  sections. 
Figure  3 shows  the  subsurface  global  numbering  for  the  helicopter  fuselage  and  rotor 
configurations.  One  important  input  for  DATA  is  the  number  of  the  elements  in  the  x 
and  y directions  for  each  subsurface,  IS.  This  input  controls  the  total  number  of  aero- 
dynamic panels  used  to  model  the  geometry  of  the  problem. 

The  variation  of  the  element  sizes  (uniform  or  nonuniform  finite  element 
breakup)  is  another  basic  input  parameter  to  the  program.  In  the  present  investigation, 
particular  attention  must  be  paid  to  the  flow  field  in  the  vicinity  of  the  rotor  shaft, 
pylon,  hub,  and  upper  fuselage  elements.  Hence,  small  elements  are  prescribed  for 
the  hub/pylon  region,  whereas  larger  elements  can  be  used  on  the  lower  fuselage 
section.  Moreover,  the  local  curvilinear  coordinates  5 and  T1  are  automatically  defined 
so  that  the  normal  to  each  surface  panel  is  always  directed  outward. 

Finally,  DATA  automatically  defines  the  parameter  KWAKES,  which  has 
the  value  1 if  the  subsurface,  IS,  generates  a wake  and  a value  of  0 otherwise.  For 
example,  KWAKES  = 1 on  the  pylon.  A vortex-layer  wake  is  assumed  to  emanate  from 
the  separation  point  and  this  is  represented  as  a doublet  layer. 

Subroutine  PREPRO  generates  a nodal  function  NOFCT  (IXY,  IS),  which 
relates  the  local  nodal  numbering,  IXY,  with  the  required  global  nodal  numbering. 
Next,  the  Cartesian  coordinates  of  these  nodes  is  obtained  by  COODPT.  Basic  inputs 
required  for  COODPT  are  the  shape  and  overall  dimensions  of  the  fuselage,  pylon,  and 
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hub  sections,  blade-twist  distribution,  blade  thickness,  chord  distribution,  number  of 
rotor  blades,  rotor  diameter,  and  extent  of  root  cutout  (shank).  Subroutine  GEOMET 
automatically  generates  a parameter  NODE  as  a function  of  the  corner  point  ICORNR 
for  each  aerodynamic  panel  of  the  helicopter.  Finally,  VEC123  completes  the  geometry 
of  the  aerodynamic  panel  (hyperboloidal  element)  that  contains  the  four  corner  points 
of  the  element.  This  results  in  continuity  of  the  complete  helicopter  geometry.  Note 
that  the  equation  that  represents  a hyperboloidal  element  is  given  by 

p = Po  + P]?  + P211  + P3?T1  (1) 

where 


Ppp'  ppmr  pmp'  and  pmm  denote  Cartesian  location  of  the  four  corner  points  and 
Pq  9've$  the  centroid  of  the  element.  This  subroutine  also  performs  the  Prandtl-Glauert 
transformation  and  the  rotation  of  the  axis  due  to  angle  of  attack  and  angle  of 
sideslip.  PRINTA  writes  the  specifications  of  the  problem,  node  coordinates,  the 
centroid  of  each  element,  nodal  numbering  for  each  surface,  and  the  nodal  numbering 
of  the  corner  points  of  the  element. 

Subroutine  CHECK  verifies  if  the  maximum  permissible  number  of  elements 
along  the  x and  y directions  on  each  subsurface  and  the  total  number  of  elements  are 
exceeded.  Several  compatibility  conditions  exist  between  the  data.  If  an  inconpatible 
relationship  is  present,  an  error  code  will  be  printed  and  the  execution  of  the  program 
terminated. 
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3.1.2  COEFFICIENT  MATRIX 


Subroutine  COEFF  forms  the  matrix  coefficients  by  evaluating  the  source 
and  doublet  integrals  over  the  surface  of  the  aerodynamic  panels.  The  effect  of  the 
rotor  wake  and  the  presence  of  separation  are  automatically  included.  COEFF  yields 
a system  of  N linear  equations  with  N unknown  velocity  potentials  at  the  centroid 
of  N elements.  SOLUTN  calls  GELG  to  solve  N simultaneous  linear  equations.  To 
obtain  the  velocity  potential  at  the  node  of  each  element,  an  averaging  scheme  is 
employed  by  AVERAG.  Using  the  value  of  the  velocity  potential  at  the  nodes,  PHI 
provides  a continuous  distribution  of  the  velocity  potential.  Note  that  the  equation 
that  represents  this  continuous  distribution  is  given  by 

cp  = 9q  + CP2 T)  + CP3 § T1  (3) 

where 


co  , cp  , cp  , and  cp  denote  the  velocity  potential  at  the  nodal  points  and  cpn 
Ypp  Tpm'  Tmp'  Tmm  r u 

gives  the  value  of  cp  at  the  centroid  of  the  element. 

3.1.3  PRESSURE  DISTRIBUTION  AND  FORCE 

VELXYZ  evaluates  the  perturbation  velocity  (Cartesian  coordinate  system) 
at  the  centroid  of  each  aerodynamic  panel.  Subroutine  CPUNR  uses  Bernoulli's  equa- 
tion to  compute  the  pressure  distribution  at  the  centroid.  PRINTB  writes  the  distribution 
of  sources,  doublets,  velocity  potential , velocity,  and  pressure  distribution.  Finally, 


FORCE  evaluates  the  total  lift  coefficient  and  induced  drag  on  the  helicopter 
configuration. 
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A brief  description  of  the  main  program  and  subroutines  is  given  in  Table  1 . 

3.2  FLOW  CHART 

In  order  to  provide  a general  understanding  of  the  overall  program,  the  flow 
chart  for  SHAPES  is  shown  in  Figure  4. 

3.3  COMMON  BLOCK  STORAGE 

Most  FORTRAN -related  variables  in  Program  SHAPES  used  by  more  than  one 
subroutine  are  organized  into  a number  of  common  blocks  as  shown  in  Table  2. 

3.4  FORTRAN  VARIABLES 

Table  3 presents  definitions  of  all  principal  FORTRAN  variables  used  in 
Program  SHAPES.  Where  appropriate,  mathematical  definitions  are  also  indicated. 

The  units  of  each  variable  are  those  used  internally  by  SHAPES,  and  occasionally  differ 
from  the  input  units. 
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Table  1 . Description  of  Main  Program  and  Subroutines. 


Preprocessor 

MAIN 

Controls  the  logical  flow  of  information  supplied  by  the 
subroutines. 

DATA 

Reads  input  data. 

PREPRO 

Defines  nodal  numbering  for  the  helicopter.  (The  surface  of  the 
helicopter  configuration  is  divided  into  34  subsurfaces  with  a 
maximum  of  four  rotor  blades  permitted.) 

COODPT 

Defines  and/or  reads  in  the  Cartesian  coordinates  of  the  nodes 
on  the  surface  of  the  helicopter  configuration. 

GEOMET 

Defines  the  four  nodal  corners  for  each  aerodynamic  panel. 

VEC123 

Defines  the  equation  of  each  hyperboloidal  surface  (i.e., 
aerodynamic  panel). 

PRINTA 

Writes  specification  of  the  problem,  nodal  numbering  for  the 
helicopter  configuration,  nodal  coordinates,  and  the  centroid 
of  the  elements. 

CHECK 

Verifies  if  the  maximum  defined  in  the  main  program  is  exceeded. 

DEBUG 

Writes  error  message. 

Coefficient  Matrix 

COEFF 

Forms  the  coefficient  matrix. 

SOLUTN 

Calls  GELG  to  obtain  the  perturbation  aerodynamic  potential. 

GELG 

Solves  system  of  general  equations  by  Gauss  elimination. 

AVERAG 

Obtains  potential  distribution  at  the  nodes  by  using  an  averaging 
scheme. 

PHI 

Defines  a continuous  distribution  of  velocity  potential. 

Pressure  Distribution  and  Force 

VELXYZ 

Evaluates  the  velocity  (perturbation)  distribution  at  the  centroid 
of  each  aerodynamic  panel. 
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CPLINR 

Evaluates  the  pressure  distribution  at  the  centroid  of  each 
aerodynamic  panel  (Bernoulli's  equation). 

FORCE 

Evaluates  the  lift  and  induced  drag  coefficients  on  the  heli- 
copter configuration. 

PRINTB 

Writes  the  distribution  of  source,  doublet,  wake,  velocity 
potential,  perturbation  velocities,  pressure,  lift,  and 
induced  drag . 

Miscellaneous 

ASINH(X) 

Evaluates  sinh  ^(X)  . 

ATANP(X) 

Evaluates  principal  part  of  tan  ^ (X)  . 

LOG(X) 

Evaluates  ln(X)  . 

Reads  in  Specification  of  the  Helicopter 
(DATA) 


Defines  Aerodynamic  Panel  Nodal  Numbering  for 
Conventional  Single  Rotor  Helicopter 
in  Local  System 
(PREPRO) 


Detines  Nodal  Point  Locations  in  Cartesian 
(Global)  Coordinates 
(COODPT) 


Defines  Corner  Location  of  Each 
Element  in  Global  System 
(GEOMET) 


Evaluates  pg,  P2,  and  £3* 
Performs  Prandtl-Glauert  Transformation; 
Incorporates  Rotation  of  Axis  due  to 
Angle  of  Attack  and  Sideslip 
(VEC123) 


Evaluates  Source  and  Doublet  Integrals 
and  Forms  Coefficient  Matrix 
(COEFF) 


Includes  Rotor  Wake  Doublet 
Integrals  in  Coefficient  Matrix 
(COEFF) 


Figure  4.  Flow  Chart  for  SHAPES  Program. 
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Figure  4.  Flow  Chart  for  SHAPES  Program  (Concluded). 


Table  2.  Program  SHAPES  Common  Block  Organization. 
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CPLINR 


Table  3.  FORTRAN  Variables. 


Variable 

Name 

Maximum 

Dimension 

- 

Definition 

AA 

(NESQ) 

A coefficient  matrix  (stored  columnwise). 

AVG 

(NNODE) 

The  number  of  elements  surrounding  the  node  INODE. 

BC 

(NELEM) 

Downwash  on  the  centroid  of  lELEMth  element. 

CP 

(NELEM) 

Pressure  coefficient  at  the  centroid  of  the  lELEMth 
element. 

HCSI 

01) 

Defines  the  normalized  nodal  line  locations  along  the 
x direction  for  surface  IS  • 

HETA 

(ID 

Defines  the  normalized  nodal  line  locations  along  the 
y direction  for  surface  IS  . 

ISFACE 

(NS) 

Defines  the  global  numbering  of  the  surface  IS. 

KNORML 

(NS) 

Indicates  the  direction  of  the  unit  normal  on  surface  IS. 

= 1 outward 
= -1  inward 

KROTOR 

(NELEM) 

= 1 if  the  element  is  on  the  rotor 
= 0 otherwise 

KROTORS 

(NS) 

= 1 if  the  subsurface  is  part  of  the  rotor  assembly 
= 0 otherwise 

KWAKE 

(NELEM) 

= 1 if  the  element  is  in  contact  with  a wake,  i.e., 
rotor  blade,  hub,  pylon,  etc. 

= 0 otherwise 

KWAKES 

(NS) 

= 1 if  the  surface  IS  generates  a wake 
= 0 if  otherwise 

NELEM 

Total  number  of  elements. 

NNODE 

Total  number  of  nodes  on  the  surfaces  considered. 

NODE 

(4,  NELEM) 

Defines  the  nodal  numbering  at  the  four  corners  of 
the  element  IELEM. 

NOFCT 

(NXYMP,34) 

Defines  the  nodal  numbering  of  the  nodes  on  the  surface 
IS. 

NS 

Total  number  of  surfaces  considered  . 
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Table  3.  FORTRAN  Variables  (Continued). 


Variable 

Name 

Maximum 

Dimension 

Definition 

NT 

(34) 

Subtotal  of  the  number  of  elements  up  to,  but  not 
including,  the  surface  IS. 

NTMAX 

Maximum  number  of  elements  permissible  on  the 
helicopter. 

NX 

(34) 

Number  of  elements  along  the  x direction  of  the 
surface  IS. 

NXMAX 

Maximum  number  of  elements  permissible  along  the 
x direction  on  any  subsurface  . 

NXY 

(34) 

Total  number  of  elements  on  the  surface  IS 
[ NXY(IS)  = NX(IS)  * NY(IS)] 

NXYMP 

Defines  the  maximum  number  of  elements  permitted  on 
a surface 

= NXMAX*NYMAX 

NY 

(34) 

Number  of  elements  along  the  y-direction  of  the 
surface  IS. 

NYMAX 

Maximum  number  of  elements  pe  missible  along  the 
y direction  on  any  subsurface. 

PC 

(3,NELEM) 

x,y,z  coordinates  of  centroid  of  elements  in  the 
reference  Cartesian  coordinate  system. 

PI 

(3,NELEM) 

A vector  along  the  ? direction. 

P2 

(3,  NELEM) 

A vector  along  the  T)  direction. 

P3 

(3,NELEM) 

A vector  normal  to  the  element  IELEM  . 

PHIC 

(NELEM) 

Velocity  potential  at  the  center  of  the  element. 

PHI  1 

(NELEM) 

Interpolation  coefficient  - the  local  variation  of 
velocity  potential  along  the  § direction. 

PHI2 

(NELEM) 

Interpolation  coefficient  - the  local  variation  of 
velocity  potential  along  the  T)  direction. 
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Table  3.  FORTRAN  Variables  (Concluded). 


Variable 

Name 

Maximum 

Dimension 

Definition 

PH  13 

(NELEM) 

Interpolation  coefficient  - the  local  variation  of 
downwash . 

SOR 

(NNODE) 

Velocity  potential  at  the  nodes. 

SOURCE 

(NELEM) 

Source  distributions  at  the  centroids  of  elements. 

After  calling  subroutine  GELG,  this  variable  stores 
the  velocity  potential  distributions  at  the  centroids 
of  the  elements. 

VELX 

(NELEM) 

x component  of  the  velocity  at  the  centroids  of  the 
elements. 

VELY 

(NELEM) 

y component  of  the  velocity  at  the  centroids  of  the 
elements. 

VELZ 

(NELEM) 

z component  of  the  velocity  at  the  centroids  of  the 
elements. 

XK 

(3,  NNODE) 

Coordinates  of  the  nodes  in  the  reference  Cartesian 
coordinate  system. 

25 


1 


SECTION  4 

PROGRAM  SHAPES  USAGE 

This  section  contains  a description  of  the  Program  SHAPES  data  input  in- 
cluding program  options  and  a description  of  error  codes. 

4.1  INPUT  DESCRIPTION 

All  SHAPES  data  input  is  accomplished  via  punched  cards.  The  input 
variables,  required  formats,  and  options  are  presented  in  Table  4.  Some  inputs  may 
require  more  than  one  card.  The  number  of  input  data  cards  required  varies  from  case 
to  case  depending  on  the  type  of  problem  considered.  Several  compatibility  condi- 
tions exist  among  some  of  the  data.  If  an  incompatible  relationship  is  encountered, 
an  error  code  will  be  printed  and  the  execution  terminated.  A typical  data  setup  for 
running  multiple  cases  is  shown  in  Figure  5. 


1 

l 
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Data  1: 

Data  2, 

Data  5: 

Data  6: 


Data  7: 


Table  4.  Input  Description. 


FORMAT  (15)  (Main  Program) 

NCASE  Number  of  cases  to  be  run 

3,  and  4:  FORMAT  (20A4)  (Subroutine  DATA) 

GG(20)  Name/  location/  date,  descrip- 

tion of  job,  and  remarks 


FORMAT  (2A4,  1015)  (Subroutine  DATA) 

GG(2)  Remarks 

KREAD  Control  code  for  the  definition 

of  nodal  coordinates 
= 0 the  coordinates  will  be 
defined  automatically 
= 1 the  coordinates  will  be 
read  in 

FORMAT  (2A4, 1015)  (Subroutine  DATA) 

GG(2)  Remarks 

NS  Number  of  surfaces  considered 


KSYMMY  Symmetry  condition  between  left- 

and  right-hand  side  of  helicopter 
= -1  if  geometry  is  antisymmetric 
= 0 if  geometry  has  no  symmetry 

= 1 if  geometry  is  symmetric 

KSYMMZ  Symmetry  condition  between  upper 

and  lower  surfaces  of  helicopter 
= -1  if  geometry  is  antisymmetric 
= 0 if  geometry  has  no  symmetry 

= 1 if  geometry  is  symmetric 

FORMAT  (2A4/10I5)  (Subroutine  DATA) 

GG(2)  Remarks 


NPYLON 


Pylon  is  considered  for  analysis 
= 1 yes 
= 0 no 


NBODY1  Fuselage  nose  section  is  considered  for  analysis 

= 1 yes 
= 0 no 
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Table  4.  Input  Description  (Continued). 


NBODY2 

NBODY3 

NVTAIL 

NSHAFT 

NHUB 

NSHANK 

NBLADE 


Fuselage  midsection  is  considered 
for  analysis 
= 1 yes 
= 0 no 

Fuselage  tail  section  is  considered 
for  analysis 
= 1 yes 
= 0 no 

Vertical  tail  is  considered  for  analysis 
= 1 yes 
= 0 no 

Rotor  shaft  is  considered  for  analysis 
= 1 yes 
= 0 no 

Rotor  hub  is  considered  for  analysis 
= 1 yes 
= 0 no 

Blade  shank  is  considered  for  analysis 
= 1 yes 
= 0 no 

Rotor  blade  is  considered  for  analysis 
= 1 yes 
= 0 no 


Data  8:  FORMAT  (2A4, 1015) 

GG(2)  Remarks 


(Subroutine  DATA) 


Data  9: 


KPYL1  Number  of  elements  in  the  x direction 

of  the  mid-fuselage  before  the  pylon 

KPYL2  Number  of  elements  in  the  x direction 

of  the  mid-fuselage  after  the  pylon 


FORMAT  (2A4,7F8.3) 
GG(2) 

UMACH 

OMEGA 

AREA 


Remarks 

Freestream  Mach  number 

Rotor  rotational  speed  (qom) 

Reference  area  for  evaluating 
aerodynamic  coefficient 


(Subroutine  DATA) 
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Table  4.  Input  Description  (Continued). 

Data  10:  FORMAT  (2A4,7F8. 3) 

(Subroutine  DATA) 

GG(2) 

Remarks 

ALFA 

Angle  of  attack 
(degrees) 

BETA 

Angle  of  sideslip 
(degrees) 

If  NPYLON  - 1,  use  Data  11;  otherwise,  skip. 

Data  11:  FORMAT  (2A4,7F8.3) 

(Subroutine  DATA) 

GG(2) 

Remarks 

XPYCTR  ’ 

YPYCTR 

Coordinates  of 
j pylon's  center 

ZPYCTR 

RXPYL  ] 

( 

RYPYL 

Radii  of 
pylon  section 

RZPYL  ] 

If  NVTAIL  = 1,  use 

Data  12;  otherwise,  skip. 

Data  12:  FORMAT  (2A4,7F8.3) 

(Subroutine  DATA) 

GG(2) 

Remarks 

VSPAN 

Span  length  of  vertical  tail 

XLEZV 

x coordinate  of  vertical  tail  leading  edge 

XTEZV 

x coordinate  of  vertical  tail  trailing  edge 

TANLEV 

Tangent  of  leading  edge  sweep 

angle 

TANTEV 

Tangent  of  trailing  edge  sweep  angle 

TAUV 

Thickness  ratio  of  vertical  tail 

ZPV 

Height  of  root  chord  of  vertical  tail  with 
respect  to  global  Cartesian  coordinate 

system 
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Table  4.  Input  Description  (Continued). 
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See  Figure  6 


See  Figure  6 


If  (NBODY1  + NBODY2  + NBODY3)  f 0,  use  DATA  13-16; 
Data  13:  FORMAT  (2A4,  7F8.3) 

GG(2)  Remarks 

XNOSE 
XBDJ 
XBD2 
XTAIL 

Data  14:  FORMAT  (2A4,7F8.3) 

GG(2)  Remarks 

YNOSE 
YBD1 
YBD2 
YTAIL 

Data  15:  FORMAT  (2A4,7F8.3) 

GG(2)  Remarks 

ZNOSE 
ZBD1 
ZBD2 
ZTAIL 

Data  16:  FORMAT  (2A4,7F8.3) 

GG(2)  Remarks 

RYBD1 
RZBD1 
RYBD2 

RZBD2  ' See  Figure  6 


See  Figure  6 


Radii  of 
fuselage  section 


otherwise,  skip. 
(Subroutine  DATA) 


(Subroutine  DATA) 


(Subroutine  DATA) 


(Subroutine  DATA) 
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Table  4.  Input-  Description  (Continued). 


Coordinates  of  hub  center 


If  NSHAFT  = 1 , use  Data  17;  otherwise,  skip. 

Data  17:  FORMAT  (2A4,7F8.3) 

GG(2)  Remarks 

RSHAFT  Radius  of  rotor  shaft 

L SHAFT  Length  of  rotor  shaft 

If  NHUB  = 1,  use  Data  18;  otherwise,  skip. 

Data  18:  FORMAT  (2A4,7F8.3) 

GG(2)  Remarks 

XHUBCR 
YHUBCR 
ZHUBCR 
RXHUB 
RYHUB 
RZHUB 

If  NSHANK  =1,  use  Data  19;  otherwise,  skip. 

Data  19:  FORMAT  (2A4,7F8.3) 

GG(2)  Remarks 

RSHANK  Radius  of  blade  shank 

LSHANK  Length  of  blade  shank 

If  (NSHAFT  + NHUB  + NSHANK  +NBLADE)  / 0 , use  Data  20 
Data  20:  FORMAT  (2A4,7F8.3) 

GG(2)  Remarks 

RROTOR  Blade  tip  radius 

YCUT  Blade  root  cutout 

XBLE 


Radii  of  hub  section 


XBTE 


x coordinate  of  root  chord 
leading  edge 

x coordinate  of  root  chord 
trailing  edge 


(Subroutine  DATA) 


(Subroutine  DATA) 


(Subroutine  DATA) 


■23;  otherwise,  skip. 
(Subroutine  DATA) 


Reference  blade  at 
azimuth  angle  of  90° 
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Table  4,  Input  Description  (Continued). 


Data  21 : 


Data  22: 


Data  23: 


Data  24: 


TANBLE 

Tangent  of  leading  edge  sweep  angle 

TANBTE 

Tangent  of  trailing  edge  sweep  angle 

TAUBL 

Thickness  ratio  of  blade  airfoil 

FORMAT  (2A4,7F8.3) 

(Subroutine 

GG(2) 

Remarks 

THET75 

Blade  collective  pitch  at  three-quarter 
radius 

THET1C 

Blade  lateral  cyclic  pitch 

THETIS 

Blade  longitudinal  cyclic  pitch 

CONING 

Blade  coning  angle 

AZIMUTH 

Blade  azimuth  angle  (degrees) 

FORMAT  (2A4,7F8.3) 

(Subroutine 

GG(2) 

Remarks 

RPITCH 

Collective  pitch  angle  at  blade 
root  (degrees) 

TWIST 

Blade  twist  (degrees) 

FORMAT  (2A4,10|5) 

(Subroutine 

GG(  2) 

Remarks 

K BLADE 

Number  of  blades 

1 TWIST 

Defines  blade  twist  distribution 

FORMAT  (2A4,10|5) 

(Subroutine 

GG(2) 

Remarks 

NWAKPY 

Flow  separation  considered 
for  analysis 
= 1 yes 
= 0 no 
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Table  4.  Input  Description  (Continued). 


NWAKHB 

Rotor  hub  wake  is  generated  for  analysis 
= 1 yes 
= 0 no 

NWAKSK 

Blade  shank  wake  is  generated  for  analysis 
= 1 yes 
= 0 no 

NWAKBL 

Rotor  blade  wake  is  generated  tor  analysis 
= 1 yes 
= 0 no 

If  (NWAKHB  + NWAKSK  + NWAKBL)  /0,  use  Data  25;  otherwise,  skip. 

Data  25:  FORMAT  (2A4,15,2F8. 3)  (Subroutine  DATA) 

GG(2)  Remarks 

N SPIRAL  Number  of  elements  along  one  rotor 

wake  spiral 

SPIRAL  Number  of  rotor  wake  spirals 

UWAKE  Induced  rotor  velocity  = V 1/2  Cy  DR 

If  NWAKPY  = 1,  use  Data  26-27;  otherwise,  skip. 

Data  26:  FORMAT  (2A4, 1015)  (Subroutine  DATA) 

GG(2)  Remarks 

NSTAG  Vortex  layer  from  separation  line 

is  considered  for  analysis 
= 1 yes 
= 0 no 

NVORT  Isolated  vortex  (branch  wake)  is 

considered  for  analysis 
= 1 yes 
= 0 no 

Data  27:  FORMAT  (2A4,7F8.3)  (Subroutine  DATA) 

GG(2)  Remarks 

CSTAG  Intensity  of  the  vortex  layer 

from  separation  line 

CVORT  Intensity  of  the  isolated  vortex 

(branch  wake) 
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Table  4.  Input  Description  (Continued). 


Data  28:  FORMAT  (2A4, 1015)  (Subroutine  DATA) 

GG(2)  Remarks 

KPRINT(I)  Output  control  code 

= 1 yes 

= 0 no 

I = 1 specification  of  the  problem 
I = 2 nodal  numbering  of  surfaces  and  corner 
nodal  numbering  of  elements 
I = 3 Cartesian  coordinates  of  the  nodes 

I = 4 Cartesian  coordinates  of  the  centroids 

of  the  element 

I = 5 coefficient  matrix  AA 

I = 6 source  integrals 

I = 7 velocity  potential  distributions 

I = 8 perturbation  velocity  distributions 

I = 9 pressure  coefficients 

I = 10  lift  and  induced  drag  coefficients 

If  NBODY1  = 1,  use  Data  29;  otherwise,  skip. 


Data  29:  FORMAT  (2A4, 1015) 


(Subroutine  DATA) 


GG(2)  Remarks 

NX(KS)  Number  of  elements  in  x direction 

NY(KS)  Number  of  elements  in  y direction 

KNSELE  Define  the  variation  of  element  size 

= 0 define  by  input 
= 1 uniform  along  x direction 
= 2 quadratic  along  x direction 

KNSSHP  Defines  nose  shape 

= 1;  r = R 
= 2;  r = R (SI1''3 
= 3;r  = R [5l1/4 

KNSTYP  Defines  nose  type 

= 1 circular  cross  section 
= 2 elliptical  cross  section 
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Table  4.  Input  Description  (Continued). 


If  NBODY2  = 1,  use  Data  30;  otherwise,  skip. 


Data  30:  FORMAT  (2A4,10|5) 


(Subroutine  DATA) 


GG(2)  Remarks 

NX(KS)  Number  of  elements  in  x direction 

NY(KS)  Number  of  elements  in  y direction 

KBDELE  Defines  the  variation  of  element  size 

= 0 define  by  input 
= 1 uniform  along  x direction 
= 2 quadratic  along  x direction 


KBDSHP  Defines  the  fuselage  shape 

= 1 cylindrical 

KBDTYP  Defines  the  fuselage  type 

= 1 circular  cross  section 
= 2 elliptical  cross  section 

If  NBODY3  = 1,  use  Data  31;  otherwise,  skip. 

Data  31:  FORMAT  (2A4,10I5)  (Subroutine  DATA) 

GG(2)  Remarks 

NX(KS)  Number  of  elements  in  x direction 


NY(KS)  Number  of  elements  in  y direction 

KTNELE  Defines  the  variation  of  element  size 

= 0 define  by  input 
= 1 uniform  along  x direction 
= 2 quadratic  along  x direction 

KTNSHP  Defines  aft-body  shape 

= 1;  r = R/?"" 

= 2;r  = R [?]1/3 
= 3;  r = R [ ?]1/4 

KTNTYP  Defines  aft-body  type 

= 1 circular  cross  section 
= 2 elliptical  cross  section 
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Table  4.  Input  Description  (Continued) 


If  NPYLON  = 1,  use  Data  32;  otherwise,  skip. 
Data  32:  FORMAT  (2A4, 1015) 


(Subroutine  DATA) 


GG(2) 

NX(KS) 

NY(KS) 

KPYELE 


KPYSHP 

KPYTYP 


Remarks 

Number  of  elements  in  radial 
direction 

Number  of  elements  in  circumferential 
direction 

Defines  the  variation  of  element  size 
= 0 define  by  input 
= 1 uniform  element  distribution  along 
radial  and  circumferential  direction 
= 2 quadratic  along  radial  direction 

with  uniform  element  distribution  in 
circumferential  direction 

Defines  pylon  shape 
= 1 elliptical 

Defines  pylon  type 
= 1 elliptical 


If  NVTAIL  = 1 , use  Data  33;  otherwise,  skip. 
Data  33:  FORMAT  (2A4, 1015) 


(Subroutine  DATA) 


GG(2) 

NX(KS) 


Remarks 

Number  of  elements  in  x direction 


NY(KS)  Number  of  elements  in  y direction 

KVTELE  Defines  the  variation  of  element  size 

= 0 define  by  input 
= 1 uniform  element  distribution  along 
axial  and  spanwise  direction 
= 2 nonuniform  (quadratic)  element 

distribution  along  axial  and  span- 
wise  direction 

= 3 quadratic  along  spanwise  direction 
with  uniform  element  distribution 
in  axial  direction 

= 4 quadratic  along  axial  direction 
with  uniform  element  distribution 
in  spanwise  direction 
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KVTSHP 


Table  4.  Input  Description  (Continued). 


Defines  vertical  tail  shape 
— 1 circular  biconvex  airfoil 


2 define  by  following  equation 


- 3 define  by  following  equation 


Table  4.  Input  Description  (Continued). 


* 


KVTSHP  (Cont.)  where 

c = chord  length 

XLE  = x components  of  sectional 
leading  edge 

XyE  = x components  of  sectional 
trailing  edge 

y = spanwise  location  of  the  section 
S = half  span 

r = t (iyr\(2c) 

max  \4  / ° 

C = root  chord 
o 

KVTTYP  Define  vertical  tail  type 

= 1 

If  NSHAFT  = 1,  use  Data  34;  otherwise,  skip. 

Data  34:  FORMAT  (2A4, 1015)  (Subroutine  DATA) 

Remarks 

Number  of  elements  in  shaft 
direction 

Number  of  elements  in  circum- 
ferential direction 

Defines  the  variation  of  element  size 
= 0 define  by  input 
= 1 uniform  element  distribution 

Defines  rotor  shaft  shape 
= 1 cylindrical 

Defines  rotor  shaft  type 
= 1 circular  cross  section 


GG(2) 

NX(KS) 

NY(KS) 

KSHELE 

KSHSHP 

KSHTYP 
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Table  4.  Input  Description  (Continued) 


If  NHUB  = 1,  use  Data  35;  otherwise,  skip. 


Data  35:  FORMAT  (2A4, 1 0|5) 
GG(2) 

NX(KS) 

NY(KS) 


KHBELE 


KHBSHP 


KHBTYP 


(Subroutine  DATA) 


Remarks 


Number  of  elements  in  x direction 

Number  of  elements  in  circumferential 
direction 

Defines  the  variation  of  element  size 
= 0 define  by  input 
= 1 uniform  element  distribution 

Defines  rotor  hub  shape 
= 1 elliptical 

Defines  rotor  hub  type 
= 1 elliptical 


If  NSHANK  = 1,  use  Data  36;  otherwise,  skip. 


Data  36:  FORMAT  (2A4, 1015) 
GG(2) 

NX(KS) 

NY(KS) 


(Subroutine  DATA) 


Remarks 


KSKELE 


KSKSHP 


KSKTYP 


Number  of  elements  in  shank  direction 

Number  of  elements  in  circumferential 
direction 

Defines  the  variation  of  element  size 
= 0 define  by  input 
= 1 uniform  element  distribution 

Defines  blade  shank  shape 
= 1 cylindrical 

Defines  blade  shank  type 
= 1 circular  cross  section 


If  NBLADE  = 1,  use  Data  37;  otherwise,  skip. 


Data  37:  FORMAT  (2A4, 1 015) 
GG(2) 

NX(KS) 


(Subroutine  DATA) 


Remarks 


Number  of  elements  in  x direction 


r — ■ — i 


Table  4,  Input  Description  (Continued). 


NY(KS) 

Number  of  elements  in  y direction 

KBLELE 

Defines  the  variation  of  element  size 
= 0 define  by  input 
= 1 uniform  element  distribution 

- 2 nonuniform  (quadratic)  element 

distribution 

KBLSHP 

Defines  rotor  blade  shape 
= 1 circular  biconvex  airfoil 

- 2 define  by  following  equation 


- 3 define  by  following  equation 
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Table  4.  Input  Description  (Continued). 


KBLSHP  (Cont.)  where 

c = blade  chord 

><LE  = x components  of  blade 
leading  edge 

Xj£  = x components  of  blade 
trailing  edge 

r = radial  location 

R = rotor  radius 

; - w (3-  A)  ocj 

C = blade  root  chord 
o 

KBLTYP  Defines  rotor  blade  type 

If  NPYLON  = 1,  use  Data  38-39;  otherwise,  skip. 

Data  38:  FORMAT  (10F8. 3)  (Subroutine  COODPT) 

HCSI(IX)  Normalized  x direction  nodal  line  coordinates 

[ Note  that  HCSI(IX),  IX  = 1,  NXP  is  implied 
with  NXP  denoting  the  number  of  nodal  lines 
along  the  x direction] 

Data  39:  FORMAT  (10F8.3)  (Subroutine  COODPT) 

HETA(IY)  Normalized  y direction  nodal  line  coordinates 

I Note  that  HETA(IY),  = 1 , NYP  is  implied 
with  NYP  denoting  the  number  of  nodal  lines 
along  the  y direction] 

If  NVTAIL  = 1,  use  Data  40-41;  otherwise,  skip. 

Data  40:  FORMAT  (10F8. 3)  (Subroutine  COODPT) 

HCSI(IX)  Normalized  x direction  nodal  line  coordinates 

[ Note  that  HCSI(IX),  IX  = 1,  NXP  is  implied 
with  NXP  denoting  the  number  of  nodal  lines 
along  the  x direction] 
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Table  4,  Input  Description  (Continued). 


Data  41:  FORMAT  (10F8. 3) 
HETA(IY) 


(Subroutine  COODPT) 

Normalized  y direction  nodal  line  coordinates 
| Note  that  HETA(IY),  IY  = 1,  NYP  is  implied 
with  NYP  denoting  the  number  of  nodal  lines 
along  the  y direction] 


If  NSHAFT  = 1,  use  Data  42-43;  otherwise,  skip. 


Data  42:  FORMAT  (1GF8. 3) 
HCSI(IX) 


Data  43:  FORMAT  (10F8.3) 
HETA(IY) 


(Subroutine  COODPT) 


Normalized  x direction  nodal  line  coordinates 
[ Note  that  HCSI(IX),  IX  = 1 , NXP  is  implied 
with  NXP  denoting  the  number  of  nodal  lines 
along  the  x direction) 

(Subroutine  COODPT) 

Normalized  y direction  nodal  line  coordinates 
( Note  that  HETA(IY),  IY  = 1,  NYP  is  implied 
with  NYP  denoting  the  number  of  nodal  lines 
along  the  y direction] 


If  NHUB  = I,  use  Data  44-45;  otherwise,  skip. 


Data  44:  FORMAT  (10F8. 3) 
HCSI(IX) 


Data  45:  FORMAT  (10F8. 3) 
HETA(IY) 


(Subroutine  COODPT) 


Normalized  x direction  nodal  line  coordinates 
( Note  that  HCSI(IX),  IX  = 1,  NXP  is  implied 
with  NXP  denoting  the  number  of  nodal  lines 
along  the  x direction] 

(Subroutine  COODPT) 

Normalized  y direction  nodal  line  coordinates 
[ Note  that  HETA(IY),  IY  = 1,  NYP  is  implied 
with  NYP  denoting  the  number  of  nodal  lines 
along  the  y direction] 


If  NSHANK  = I,  use  Data  46-47;  otherwise,  skip. 
Data  46:  FORMAT  (10F8.3) 


(Subroutine  COODPT) 


HCSI(IX) 


Normalized  x direction  nodal  line  coordinates 
[ Note  that  HCSI(IX),  IX  = 1,  NXP  is  implied 
with  NXP  denoting  the  number  of  nodal  lines 
along  the  x direction] 
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Table  4.  Input  Description  (Continued). 


Data  47:  FORMAT  (10F8. 3) 
HETA(IY) 


(Subroutine  COODPT) 

Normalized  / direction  nodal  line  coordinates 
[ Note  that  HETA(IY),  IY  = 1,  NYP  is  irr^lied 
with  NYP  denoting  the  number  of  nodal  lines 
along  the  y direction] 


If  NBLADE  = 1,  use  Data  48-49;  otherwise,  skip. 
Data  48:  FORMAT  (10F8. 3) 

HCSI(IX) 


(Subroutine  COODPT) 


Data  49:  FORMAT  (10F8. 3) 
HETA(IY) 


Normalized  x direction  nodal  line  coordinates 
[ Note  that  HCSI(IX),  IX  = 1,  NXP  is  implied 
with  NXP  denoting  the  number  of  nodal  lines 
along  the  x direction] 

(Subroutine  COODPT) 

Normalized  y direction  nodal  line  coordinates 
[ Note  that  HETA(IY),  IY  = 1,  NYP  is  implied 
with  NYP  denoting  the  number  of  nodal  lines 
along  the  y direction] 


If  KREAD  = 1,  use  Data  50;  otherwise,  skip. 
Data  50:  FORMAT  (3E1 2.4) 

XK(K,  INODE) 


(Subroutine  COODPT) 


The  nodal  coordinates  in  global  Cartesian  system 
( Note  that  (XK(K,INODE),  K = 1,3),  INODE  = 1, 
NNODE)  is  implied] 

If  NWAKPY  = 1 and  NSTAG  = 1,  use  Data  51-54;  otherwise,  skip. 

(Subroutine  COODPT) 


The  separation  wake  nodal  coordinates 
in  global  Cartesian  system  [ Note  that 
K=1 ,3  is  implied] 


Data  51: 

FORMAT  (3E12.4) 

YPP(K) 

Data  52: 

FORMAT  (3E12.4) 

YPM(K) 

Data  53: 

FORMAT  (3E12.4) 

YMP(K) 

Data  54; 

FORMAT  (3E12.4) 

YMM(K) 
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Table  4.  Input1  Description  (Concluded). 

If  NWAKPY  = 1 and  NVORT  = 1,  use  Data  55-58;  otherwise,  skip. 

(Subroutine  COEFF) 

The  isolated  vortex  branch  wake  nodal 
coordinates  in  global  Cartesian  system 
[ Note  that  K=l,3  is  implied] 


YMM(K) 


Data  55: 

FORMAT  (3E12.4) 

YPP(K) 

Data  56: 

FORMAT  (3E12.4) 

YPM(K) 

Data  57: 

FORMAT  (3E12.4) 

YMP(K) 

Data  58: 

FORMAT  (3E12.4) 
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4.2  ERROR  CODES 

To  insure  proper  data  are  read  in  and  to  minimize  unnecessary  computer 
time,  SHAPES  contains  several  debugging  statements  at  specified  checkpoints.  If  an 
inconsistency  among  input  data  or  violation  of  the  input  specification  is  detected,  an 
error  code  will  be  printed  and  the  computer  run  terminated.  Table  5 summarizes 
these  error  codes. 

Table  5.  Error  Codes  Summary. 

Code  Number  Description 

100  Mach  number  is  greater  than  1 (Present  SOUSSA 

is  a subsonic  program) 

(DATA) 

200  NS  / (NPYLON  + NVTAIL)*MULTY 

+ (NBODY1  + NBODY2  + 

(N  BOD  Y3)*MULT*MULTY 
+ (NSHAFT*MULTY 
+ (NHUB*MULT*MULTY* 

+ K BLADE  *(N  SHANK  + NBLADE*MULT) 

(DATA) 

300  Number  of  elements  in  the  x direction  of  a sub- 

surface exceeds  the  maximum  permissible  value.  This 
is  limited  by  the  storage  capacity  of  the  computer. 

The  user  has  the  option  of  changing  the  value  of 
NXMAX  and  NYMAX  defined  in  the  main  program. 
(DATA) 

400  Same  as  code  number  300  except  for  the  y direction. 

(DATA) 

500  Total  number  of  elements  required  for  analysis  ex- 

ceeds the  maximum  limit  specified  in  the  main 
program . 

(DATA) 

IER  IBM  subroutine  GELG  provides  error  code 

= 0 no  error 

= 1 no  result  is  obtained  because  the  number  of 
equations  is  less  than  1 or  the  pivoting  element 
at  any  elimination  step  equals  0 
= N warning  is  indicated  because  of  a possible  loss 
of  significant  figures  at  elimination  step  N+l 
where  the  pivoting  element  is  less  than  equal 
to  the  specified  tolerance  times  the  magnitude 
of  the  greatest  element  of  matrix  A 
(DATA) 
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SECTION  5 


PROGRAM  OUTPUT 

In  Program  SHAPES  several  output  options  are  available  and  these  are  briefly 
summarized  here.  The  amount  of  output  required  to  be  printed  out  is  controlled  by 
Data  28,  described  in  Subsection  4.1 . Among  the  information  available  is  the  specifi- 
cation of  the  problem  defined  by  the  user  and  the  basic  geometric  inputs  to  the  problem, 
i.e.,  overall  dimensions  and  shapes  used  to  model  a helicopter  configuration.  In 
addition,  the  nodal  numbering  of  the  aerodynamic  breakup  as  well  as  the  Cartesian 
coordinate  location  of  the  nodal  points  and  the  centroid  of  each  aerodynamic  panel 
can  be  requested.  Furthermore,  the  coefficient  matrix  and  source  distribution  can  be 
output.  Also,  the  perturbation  potential,  velocity,  and  the  pressure  coefficient  are 
available  as  outputs.  Finally,  the  lift  and  induced  drag  can  be  requested. 
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SECTION  6 


PROGRAM  LISTING 


SAMPLE  CASE  1: 
SINGLE -BLADE  ROTOR 


Data  Input 


9 

PWMI 

mm  rnim  nm  ipmiiiw 

SUKLi  JLSUBt  J&JLBL -FLOBLi* 

mss  • 

! t 1 1 

— . _ 

— 

•if 

t • 

• • 

• • 

MIC* 

Hfl 

MM 

. A.M1 

9ft. || 9 

ifiw 

s».m 

1.91* 

••Ml 

rntto 

ts.tlf 

•.III 

MMRC 

• • 

S 1 

MINI 

JWKJ10* 

i i 

s ? 

1 1 

t 1 

• 1 

1 

• • I 


1.1*9  I.Mfe  I.IM 

l*^«  ... 


till 


IU1 


NtiM  • ir 


IS*  I 1 KTI  P N*f 

is*  — t.  tum. — i mi  ? m*i 


?1  isrtce*  t MV*  « KHO*H  • 1 1 «««Tf*9*  | 

a ISMCC*  t «!•  II  JUUN.I  -1  KVAKCS*  I RMTOKi  4 


Specifications  of  the  Problem 
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X,  Y,  and  Z Coordinates  of  the  Nodes 
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X,  Y,  and  Z Coordinates  of  the  Centroid  of  Aerodynamic  Panels 
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Nodal  Numbering  for  Surfaces 
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The  Distribution  of  the  Source 
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Pressure  Distribution 
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SAMPLE  CASE  2: 

51A  HELICOPTER  ROTOR 
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Specifications  of  the  Problem 


POP  FART  24 
NIX=  5 
NY=  f 


POP  PART  25 
NX=  £ 

NY  = 5 


FOR  FART  27 
NX=  5 
NY=  « 


FOR  PART  28 
NX=  ? 

NY=  f 


FOR  P 6 FT  30 

NX  = £ 

N Y=  5 


FOR  PART  31 
NX=  «■ 

NY=  f 


FOR  PART  33 
NX=  5 
NY=  f 


FOR  PART  34 
NX=  c 
NY=  £ 


NELEM=2or 

KSYMt«Y=  0 
KSYMPZ=  0 
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6 8 -7.79120 
•9 — -17.  0J8  SO- 
TO -12.03000 
7i — -17.  osroo 
72  -12.03800 


-.19991 

o.oooor 

-.00300 

-.0r757 

-.0972* 


-15.07  230 

-16.89323 

-16.89320 

-16.89  370  ' 

-16.89320 

-ir;89-S70_ 

-1S. 89323 

-17.50  COO 

-17.50000 

-17-.5000T 

-17.50003 


1.0  6496 

.ooerr- 

.0  93  35 
—.17300 
• 3 67  a 1 

TO  8077— 

1.06996 
. 00300 
.04256 
-.17037 
.36321 


51 — 

-mTtrtnn 

—.33177 

- •- 

92 

-1  7.50  000 

l • 06**6 

- 

-r.  jrtnro 

.ovrsi 

' - 

9* 

-2.33  000 

.16763 

* 

9F 

-7.  33C33  " 

• 3 78  o 2 

* 

56 

-2.33000 

• 67667 

* 

97" 

-7".T9  120 

. 06181 

— 

58 

-7.79120 

• 16763 

- 

<59 

-7.79123 

. 3 78  32 

“ 

100 

-7.79120 

. 676A7 

* 

101 

-ii.i)36d(J 

• 0618 1 

“ 

1 C2 

-12.03880 

.16763 

* 

-12.03880 

-15.07780 

-15.07280 

-15.07  280 
•16.89320 
-16.89320 

-16.89373 

-16.89320 


.6766  7 

-.nnr 

.16763 
. 378  67 
.67667 
.09181 
.16763 


-.19991 

o.  oocor 

-.  00366 
-.01756 
- • 09776 
--.10309'  ~ 
-.19991 
O.  00000 
-.  0079  6 
-. T3190 
-. C7175 


. 01707 
. 04623 


. G9672 
.15215 


-.15215 

-.01707 

-.04623 

- * 15215 
-.01207 
-.04623 


-.09726 
--V10S09- 
-. 19991 
O • 00000 
-.00306 


1(2 

-.38761 

-7.79120 

-.09728 

— trs — 

- * . 68587“ 

—7 . 7912* 

- • 10909 

1 2* 

-1.06446 

-7.79120 

-. 19991 

~ 125  - 

'-.00900' 

-17.  *3389 

0.  00000  - 

126 

-.0-335 

•12.03880 

-.00388 

127 

-■--;i730i 

--12. 03330 

’ *7 01758 

128 

-.387*1 

•12.03860 

-.04726 

— 175 

-.883*7 

— 17.0  3500 

-.10307 

120 

-1.0644* 

-12.03880 

-.19941 

— xsr  - 

--.:oroo 

-15.07737 

0.00000 

132 

-.09  335 

-15.0778* 

-.  00386 

133 

-.17  0*4 

“4  XT’.  *77*9 

139 

-.3*7*1 

-15.077** 

-.04  726 

— iss- 

—  .*33*7 

-17.-977*9 

-.10399  - 

ue 

-1.0644* 

-15.17280 

-.19991 

- 137 

-.005*9- 

-1*. 99379 

TT.  00  090 

1!S 

-.04335 

-16.89320 

-.00388 

IJT 

- • l7 300 

-16 .193  79 

- . ai  25  c 

160 

-.38781 

-16.39320 

-.0*726 
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216 

2 06 

20 T 

200 

“20T 

210 

nr 

212 

217 

216 


17.50000 

750000 

17.50000 

r.-nnr 

2.23000 
-2733 
2. 23000 


7.79120 

7779120- 

7.79120 

tr.rstoo-- 

12.03680 


.06(23 
70963 
.15215 
.Ol  207  - 
.06(23 
09632 
.1521! 
.01207— 
.06(23 
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161 

—"=7  63587“ 

-12755321-“ 

—=770759 

142 

-1.06448 

-14.89320 

-.199=1 

163 

-.03  0 03 

-17.50000 

“07 0 0 0 0 C 

166 

-.06253 

-17.50000 

-.00796 

— r=3~ 

-717032— 

-=T77Tr3Tt 

-.03190 

16b 

-.38121 

-17.50000 

-. 07179 

167- 

6.68  127  " 

-177  50  005 

“=  .72562“" 

148 

-1.06448 

-17.50000 

-.19961 

169 

-.06181 

—=2757000 

“--701257 

150 

-.16763 

-2.33000 

-.C6623 

ITT- 

-.37  862 

*Z. 33000 

' =.09632 

152 

-.67667 

-2.33000 

-.15215 

1T3 

—=706181 

=7.791 28 

- • 0IZ07 

t!6 

-.16  763 

-7.79120 

-.06(23 

159 

— *378 6? 

•7.79170 

- • 09632 

156 

-.67667 

-7.79120 

-.15215 

— rrr“ 

— -.mm 

-12.038811 

0120  7 

1S8 

-.16  763 

-12.03880 

-.0=623 

159 

-.37(52 

-1Z.338BTJ 

-riJ95TZ 

160 

-.67667 

-12.03380 

-.15215 

1(1 

-7017181  - 

-1V757230" 

=751207 

162 

-.18763 

-15.07260 

-.06823 

— rn— 

-.37  5b2 

-15.J72S0 

-.09632 

166 

-.67667 

-15.07280 

-.15215 

— irr 

-706181  “ 

--15. 89527“ 

—=.51207  — 

1(6 

-.16763 

-16.89320 

-.06(23 

167— 

-.37862 

“=16759325 

09632 

168 

-.67667 

*15.89320 

-.15215 

1(9“ 

2.37000 

- • 00  0 0 0 3 . 00  0 0 0 

170 

2.  33  000 

-.0=335 

-.00388 

171 

— 2.33000 

" -.17300' 

-.01756 

172 

2.33000 

-.38231 

-.06726 

173 

2.33000 

- -.66567 

“*. 10309 

176 

2.33000 

-1.06448 

-.19941 

175 

7.79120 

-.00000 

0.00000 

176 

7.79120 

-.0=335 

-. 00388 

177 

7779125 

— *.17305— 

“=717758 

178 

7.79120 

-.36761 

-.  06  726 

179 

7779 171T 

.689TT 

—-.70309 

180 

7.79120 

•1.066=6 

-.19961 

— m— 

— 127  028  90— 

* . 3 TO  OX 

07  COOOC 

1(2 

12.03(30 

-.04335 

-.00388 

— tor 

-T2.0S860" 

.17700 

— - . 01758— 

1(4 

12.03880 

-.38781 

-.06726 

105" 

-12.03880“ 

*7  68587 

- . TO 309" * 

1(6 

12.03(80 

-1.06668 

-.19961 

107— 

to.orzao— 

— ■ .tmtns — 

— 5.000  00 

188 

15.07280 

-.04335 

-. 00388 

1(9- 

"19.07280- 

■— • trjoo  ~ 

. 0175  e — - 

190 

15.07280 

-.36/81 

-.16726 

191— 

“15.  07280 . 08587 =.10309 

192 

15.07280 

-1.064498 

-.  19941 

193— 

18. 89320 

«.3D0  08 

— 5.00000 

196 

16.89320 

-.0=335 

-. 00388 

zis  i T.trstso-- 

2»  12.03*60 

- tit — T*-.-~rrrz — 

218  15.07780 

219“  15  • 37280 

220  15.07200 

221  16.89720  ‘ 

222  16.8*329 


■*.1786? 

-.676*7 

--rrtr»m — 

-.16753 

-.37862“ 

-.67667 

-.06181 

-.16783 


-. 09632 
-.15215 
--~ci2or- 

-.06623 
t.  09632 
-.15215 
-.01207 
-.  06623 
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X,  Y,  and  Z Coordinates  of  the  Centroid  of  Aerodynamic  Panels 


1 .02167 

~2 71*817“ 

3 • 280**  1 

% ~ .536e* 

5 .d7«?!7 


5.06060 

"T.TJliOBO' 

5.06060 

*.060611 

5.06060 


-.0019*. 

- • mm 

-.0  32  <.2 

-.07rir 

-.15125 


“6 .‘  02157 T.  9 1700 - . 00190“ 

7 .10617  9.91500  -.010  73 

* .28147 — 9-.9150: . 03292 

9 .53664  9.91500  -.  075  17 

10  .67517“  9.91500 -.  151 25 

11  . 02167  13.5556  0 -.  00194 

12  .10617  13.55560  -.01073 

13  .26041  13.55580  -.03242 

14  .53684 13755560 -.07517 

15  .67517  13.55580  -.15125 


.10817  15.98300 

-.261151 15.96300 

.53684  15.98300 


.02148  17.19660 
.10731 — lT.r96S0 
.27859  17.19660 


.67403  17.19660  -.15736 


.10472  5.06060  -.02915 

• 52765  5.06060  -.12423 

.8705  8 5.05080 *.17376 

.02091  9.91505  -.0  C6C4 

.10472 9791500—  -‘.0  2915 

.27313  9.91500  -.07127 

.52755 9.01530 -.124  23 

.87058  9.91500  -.17578 


.104  7 2 1 3.5  5 58  0 -.0  29  15 

.27313  13.35580  •*-.07127 

• 52705  13.55560  -.12423 

.8715  8 13.33532  - .17578 

.02091  15.98300  -.  00604 

.27313  15.98300  -.07127 

. 52765 IT.  96330 *712473 

.87158  15.98330  -.17578 

. 021  nr ‘it.  r96so — -.  ootoi — 

.10559  17.19660  -.02455 


49  .52994  17.19660 

53 ; 87T73  17719650" 

51  -5.06060  .02167 

32 — -r.  Ob  060 710017— 

S3  -5.06060  .28041 

‘53 — *57  0616  3 .53683 

55  -5.06060  .87517 

56  “ *9. 91500 702167 

57  -9.91500  .10817 

56— -r.  91730 .28  041“ 

59  -9.91500  .53684 


-.11197 

-716963 

-.00194 
•701073' — 
-.03242 

■ . U7517 

-.15125 
--.0OT9r — 
-.01073 
*703237  “ 
-.07517 


61  -13.55580 

62  -13.  55570 

63  -13.5556 C 
63  -137  55560“ 
65  -13.55580 
63  - 15. 96300 
67  -15.98300 
Sr  - IT.  98300— 
69  -15.90300 


.02167  -.00194 

-710617 -7  01073“  ‘ 

.28631  -.03232 


.10817 

.28031- 

.53664 


-.1512? 
- . 00193 
-.01073 
-.33237 
-.07517 
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021*6 


IIOZ96 


71  -17.19660 

.021*8 

-.0  0296 

1 1 *1  / • l'J&GU 

.rU  731 

-.0153* 

73  -17.19660 

.27859 

-.0*213 

^ "l^o  l9bb  D 
75  -17. 19660 

.37403 

-.15738 

76  “-53  0606  0 - 

.02091“ 

• .HOE  04" 

77  -5.06060 

.10*72 

-.  0 291* 

70  -5.06065 

.2^313' 

-.07127 

79  -5.  0606  0 

.52765 

-.124*23 

80  -5.06060 

• 37  Q*8 

-.175  72 

dl  -9.915CC 

.02091 

-.CC6C4 

82  -9.9150  0 

310*72 

~ - * 0 29  15 

03  -9.91500 

.27313 

-.07127 

34  -T.  915C0 — 

752765 

-712*23“ 

85  -9.91500 

.87058 

-.17578 

86  -13.  5556  0 - 

~ .02091 

- . o 06  tm 

07  -13.  55580 

.10*72 

- . 0 29  15 

55  -133  55580 

-“.27313 

-•=.07127 

09  -13.555c  0 

.52765 

-.124  23 

90“~ - 13.55580 — 

•67C5S 

— *.173  78 

91  -15.98300 

.32091 

- . 0 06  04 

92-35.98300- 

710*72 

—*.02913 

93  -15.98300 

.27313 

-.07127 

9*-- 15.-90300  - 

i 5 2765 

- * . 12*23 

95  -15.96300 

• 3 7 05d 

-.17578 

• ix  j 

* • 0 03  01 

97  -17.19660 

.10559 

-.02*55 

90--I79 1966  0- 

727499 

—=.06158- 

99  -17.19660 

.5299* 

-.11197 

tOO"  -32.19660 

•67 173 

- .1696* 

101  -.02167 

-5.06060 

-.0  0194 

10£  -.105x7 

“,r9»  OS  089 

- • 0 73 

103  -.280*1 

-5.06060 

-.032*2 

tor — -.5368T- 

-9.0A090  - 

.079 17 

105  -.87517 

-5.06060 

- . 1 51  25 

106  -.32167 

-9.91500 

— -.0  0194 

107  -.  10817 

-9.91500 

-.01073 

- itra — <■.  2 os  * i — ~=939roo: «irw 

109  -.5366*  -9.91500  -.0T317 

tin  -.ar«ir  -».9i500 — -.istw 

111  -.02162  - 13.59560  -.0019* 

112  -'.  I0817“-rT.555SO -.  OIOTT 

113  -.263*1  -13.59903  -.032*2 

"TH -.53589-  *15753380 *T07317- 

115  -.87512  -13.55500  -.15125 

115  ».  02167-13.90300 -.00195 

117  -.10817  -15.90300  -.010  73 

119  -•  5366*  -15.90300  -Io75  17 

“175 -.87517  -15.98350 *0  191 29 

121  -.  021<-0  -17.19660  -.0  0296 

123 -.27859_-  17.19660  -.0*213 

125  -*07*0  3 - 17*19660  -.15738 

“125 =302591 -5.-0605  0 -'.006  5* 

127  -.10*  72  - 5.0  6 0 6 0 -.0  29  15 

128  -3  2731 3 — =5.06  060 = . 07177- 

129  -.52765  -5.06060  -.  12*23 

13  0 • • 67055  "TolB  OSD  • • l / 5 


131 

-.02391 

-9.9133C 

-.0C6C* 

172— 

1 -.10*72 

-9.-91300 

-302913 

133 



-.27313 

— • »7TK« 

-9.91500 
— sro  ai  m 

-.07127 

135 

-.07058 

•’lo’llOOO 

-9.91500 

-.17371 

136  - 

—■-'.02591 

-13.53310 

— -“.oorrs 

137 

-.10*72 

•13.55360 

-.0  29  15 

- ~ . 7i  n 3 

'- 13.53310 

*.'07127 

139 

-.32765 

•13.53360 

-.12*23 

ISO- 

- .87038  *13.33380 

==-.17378 

191 

-.02091 

-13.96300 

-.006  0* 

1*2— 

■ -.10*72 

~ 13.9T300- 

-YT2315 

1*3 

-.27313 

•13.96360 

-.0712r 
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-.52765  -15*98300 
-.87058  -15.9o30: 

-.  m-ro  - ir-.  i96SO~ 

-.10559  -17.19660 
- ,r79°S  - ir.  19668- ~ 
-.52999  -17.19660 
. 8TI73  -17.1966-0 — 
5.06360  -.02167 


5.  06060“ 
5 . O0O6  0 
5.  06363“ 
5.  06060 
“9.915-00— 
9.91300 
9.  91500 
9.91S00 


-.10817 

-.25091 

-.37386 

-.3751/ 

—202757- 

-.10817 
-.2309  7“ 
-.53688 


“9 . 91500 -167517“ 

13.  55580  -.02167 


-.12923 

-.17578 

-.00557 

-.02955 

“-.06156 

-.11197 

“-1T89K9 

-. 00199 

-.01073 

-.03292 

-.07517- 
- . 151 2* 

~-7  00159 

-.010  73 

-“.0  3297 

- . 075  17 

-.  15125 

-.0  0199 


163  13. 56569 

“166  -Ij;565e0“ 

165  13.55580 

166  - 752  98300 — 

167  15.98300 


-.2*061 

9.53686 

-.87517 

-.02167 

-.10817 


-.132  62 
-9.07517 
-.15125 
— . 00156 
-.01073 


169  15.98300 

170“  “15.95300 

171  17.19660 

172  17. 19660 

173  17.19660 


-.53689 

-.87517“ 

-.02198 

-.11737“ 

-.27859 


-.  075  17 
=115125 
-.00296 
=-.01536 
-.06213 


175  17.19660 
176- ' 5.06060- 
177  5.06060 

176  5.06065“ 
179  5.06060 


- . 157  38 
• . 006  00 
- . 0 29  15 
1 07127 
-.12623 


9.91500 


9.91500 
156  9.9 1500 
185  9.91500 


187  13.  55580  -.10972  -.  02915 

188  TT.  55560 -127313 =■ . 07127 

139  13.55580  -.52765  -.12923 

190  13.5555  0 -.37055 ~ -1 17576 

191  15.  98300  -.02091  -.  006  09 

192  15795TOO =110172 . O 29  15 

193  15.98300  -.27313  -.07127 

199  15.98300  -152765 -.12923 

195  15.98300  -.37058  -.17578 

196  17.  19660  -.02110  - . 0 05  Cl 

197  17.  19660  -.10559  -.  02955 

198  1711566  0 ==".27195 =116116 

159  17.1946C  -.52996  -.11197 

200  17.19660  -.87173  --.16SC9 
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Nodal  Numbering  for  Surfaces 


FCR 

SURFACE 

24 

1 

7 

13 

19 

25 

31 

2 

8 

14 

20 

25 

32 

3 

9 

15 

21 

27 

33 

% 

10 

15 

22 

28 

3% 

5 

11 

17 

23 

29 

35 

s rz r» 25 — — nr if 


* X 

T 

TT 

— rr 

“29“ 

TT 

37 

41 

45 

49 

53 

32 

35 

4? 

06 

50' 

50" 

3 3 

39 

43 

47 

SI 

55 

34 

45 

-4 

4 8 

52 

56 

35 

6 

12 

18 

24 

30 

36 

FCR 

SURFACE 

27 

57 

63 

69 

75 

81 

87 

j 8 

64 

TO" "" 

73 

62 

59 

6? 

71 

77 

33 

89 

6 8 

66 

TT 

78 

84 

90 

61 

67 

73 

79 

85 

91 

FOR  SURFACE 28 


57— 

-b3~ 

65  TT 

TT 

9T 

93 

97 

101 

105 

109 

38 

94 

98 

107 

106 

117 

85 

95 

99 

103 

107 

111 

90 

' 96 

Tiro — 

tuo- 

ro8 

112 

91 

62 

68 

74 

80 

86 

92 

FCR 

SURFACE 

30 

113 

119 

125 

131 

137 

143 

rro — 

120 

176" 

132 

138 

115 

121 

127 

133 

139 

145 

116 

122 

128 

130 

146 

146 

117 

123 

129 

135 

141 

147 

118 

17% 

I JO Too 

t%2“ 

1%6 

— 113 

119 

- 1 29 

— 131 

TS7 

T43 

149 

15  3 

157 

l&L 

165 

144 

150 

154 

— 158” 

— 162 — 

166  — 

145 — 

151 

155 

159 

163 

167 

146 

IT? 

156 

i60 

15% — 

1 66 

1%7 

118 

124 

130 

136 

1*2 

148 

FOR 

SURFACE 

33 

169 

175 

181 

187 

193 

199 

170 

If  6 

l 87' 

168 

19% 

20U 

171 

177 

1 83 

189 

195 

201 

172  173 

T84‘ 

- 190 

196 

27  r 

173 

179 

185 

191 

197 

203 

174 

188 

- 186 

192 

195  ” 

204 

FCITTURFACE 

“ 34 

169 

175 

l 81 

187 

193 

199 

205 

209 

213 

217 

221 

203 

206 

210 

214 

213 

2 22~ 

201 

207 

211 

215 

219 

223 

202 

20  e 

212 

215 

223 

224 

2"3 

174 

loO 

185 

192 

196 

204 

64 


'Mi* 


105 

12* 

118 

123 

117 

— ior- 

125 — ms — 

175“ 

IT9 

107 

127 

121 

12$ 

120 

108~ 

128“ 

172 

127 

" 1 21 

109 

129 

123 

128 

122 

rxo 

ITU 

"124 

~T29 

“"123 

in 

132 

12$ 

131 

125 

nr 

13T  “ 

127 

132 

113 

134 

128 

133 

127 

114 

135  — 

" 1 25 

134“ 

178 

115 

136 

130 

135 

129 

1X6 

' 138 

"172 

137 

131 

117 

139 

133 

138 

132 

— ur 

tm — m — 

X39- 

X3T 

119 

141 

135 

140 

134 

I2ff 

~mz 

1 36 

141 T35 

121 

144 

138 

143 

137 

122 

145 

139 

144 

136 

123 

146 

140 

145 

139 

12*  - 

t*7“ 

-m — 

“146 — r*n 

12$ 

146 

142 

14  7 

1*1 

126 

149 

193 

~ X13 

1T9 

127 

ISO 

154 

149 

153 

129 

152 

156 

151  155 

TIT 

TX8“ 

-12*— 

152 T56 

131 

153 

157 

119  125 

132 

— T68 

“1*3  r 57 

133 

155 

159 

15*  158 

1*1  161  IS*  131  ur 


The  distribution  of  the  Source 


FCR  SUBSURFACE 

24 

-.67898E*C1 

-.129657*02 

-.176017*02 

-.202377*02 

-.11J51E*02 

-.6<*9ft4E»01 

-•12450E ♦ 02 

-•lb917E*02 

- • 19  516E  *02 

• • 1 09b*E  ♦ 0 Z 

-.367C5£>01 

" #898267*01 

-.995921*01 

”-.10  0797  * 02 

-.620317*01 

• 2 0676E+01 

• 4 0 146E  ♦ 0 1 

.546357*01 

• 63 976E * 01 

.336527*01 

- .761457*01- 

-~axm>E*V2 

.201377*02  " 

- .232617*  C2 

.127507*02 

FOR  SUBSURFACE  25 


-.70504E+01  — 

-.  135217  *02 

- .193577*02 

- -.211717*02 

- .110647*02 

-.75670E*01 

-.147847*02 

-.202207*0  2 

- .23463E*02 

-.131027*02 

- *5^701 t ♦Ol 

- . 113927*  02 

-;1B2C3E»TJT  ’ 

- • 18  9BBE>02 

• • IT-  B64E+D  2 

-.8F8O5E^00 

- • 20  23CE ♦ 0 1 

-.29820E*01 

- . 37  337E*  01 

-.22b39E*0l 

•603%3E401 

•115S6S»02 

1156  637*02  - 

.179  337  * 02 

.986577*01 

FOR  SUBSURFACE 

27 

-.67b967*01 

-.129857*02 

-•11351E»02 

-•l7RQiE*02 

• *20  237E*02 

• •6*'  984F  ♦ 01 

-.12450E^02 

-.169177*02 

- . 19  51  €£♦  0 2 

-.10964E4-02 

-.36705E+01 

- .69626t^  0 i 

-.965927*01 

-.10 0797*02 

-.620317*01 

• 2067bE  *01 

• 40 146E ♦O  1 

• 54635E^  01 

• €3  976E*01 

• 336  527*01 

.771457*01 

•I 4 5D  BE  ♦ IT? 

.201377*02-"" 

T23261E»02 

.127507*02 

FOR  SUBSURFACE  28 

-*.7  85847*91 -.135217*82 


-.7Sb70E*01  -.14784E*02  -.202207*02  -.23  463E*02 


* .597017*01 -Ttl7927»02 


-.058957*00  -.20  23  0E*C 1 -.290207*01 

•603437*01  .115067*92 .156637*07 1179337*02 


-.118447*02 

-.13  1 827  * 02 
— - .106647*92 

- . 226397*  0 i 
*986577* 01 


* 
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F.CR  SUBSURFACE  30 

-.67#90E*01  -.1Z965E*0Z  -.17601E»0Z  -.Z0Z37E*0Z  -.11351E»0Z 

-•6».<»«-E  *01  -.1Z»50E*0Z  -.16917E*0Z  -.195l6E*0Z  -.1C964E»0Z 

-.36T05E*0l  - • 6932EE  + 0 1 — -79459ZE*0I  -.10879E»0Z  -.6Z031E*01 

«ZB676E»gl «40146£+01  .5463&E+01 .E3976E«01  . 336S2E*H 

— .7M4YEH1I .148060*77 .Z0137E*OZ .7J741E*02 .'1Z750E*0Z 

FCR  SUBSURFACE  31 

~-;705BVE*01 175Z1F*B? .I83B7EMTZ -.Z11TEE*07 -mB4I.E*tlZ 

-.75670E*01  -t»W»8*K  -.Z0ZZ0£*0Z  -.ZJ463E*07  -.1310ZE*OZ 

-.59701E*01  -.1I79ZE*0Z -.16Z03E*1IZ -.18990E+0Z  -.10664E»BZ 

- . 85895c  ♦ 00  -.Z0Z36E*01  -.Z98Z0E*01  -.37337E*01  -,ZZ639E*01 

.SC343F+01  .115S6E*T>Z .15863E*0Z  .17933E*0Z  .986E7E*01 

FOR  SUBSURFACE  33 


6769EE>01 

-.lZ95Sr»07 

-’•17601E+02 — 

~-.ZBZ37E*BZ 

-.11351E+0Z 

6i9S<»Em 

- .1Z456E*07 

-.16917E*0Z 

1951EE*07 

-.10964E»0Z 

367C$E»tl 

-.698Z6E*01 

-.9499ZE»01 

-.10B79E*0Z 

- *670  31 E*  0 1 

?Cb7ftF4?l 

• »*0 148E  ♦ 01 

• 94»fi35£4  01 

• S3  97br»01 

• 3365 ZE ♦ 0 1 

76145E*C1 

•143C5E+0Z 

* ?T)  137E+  0? 

• 73Z41E*  OZ 

• 17750E*  0 Z 

FOR  SUBSURFACE  34 

-«.Twrr*rm  — *vi35zir»rz — *.i8387F*c7 *rznrsE*77  — -.no  me»oz 

-.78670F*01  -«14784E»02  -.ZOZZOE»OZ  -.Z3463E*0Z  -.1310ZEOZ 

* .S970IT*OI «.1I79ZE»0Z ■*  ,l«7irSt»r7 » 7H9S0E»0Z * .1 C689E*TZ 

-.8S898E*0#  -.Z0ZS8E»01  -,7S8Z0E*01  -.37337E*C1  -.ZZ639E»01 

-;M54SFMrt  .ll?8et*07 .17933F*C7 .»8»9rE»Pl 
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The  Distribution  of  the  Velocity  Potential 


FOR  SUB SURF  ACE 

24 

•mo2r*o3 

.?0640r»t)3 

.2tTJ2E*03 

.100b0t*03 

•l5bt8E*03 

•178B6E*C3 

.21  339F^ C 3 

•2l796E*C3 

« 1*  379E*  C3 

• 167t6E*0  3 

.15350E»03 

.E?3E?E*0J 

• ?3153E*0 3 

.20  9E>1E  ♦ 03 

• 17963E*  03 

.lft<W<*E*03 

.233StE*C3 

.2b566E«C3 

• 22  568E  ♦ C 3 

. If 938E.03 

.16643E*0J  — 

.23073E *03 

• 183b6E*0  3 

«24566f*03 

' • 2Z257E*C3 

FOR  SUBSURFACE 

25 

.U0HE*0J 

•175B3E*03 

Vr7fl*0E*03  ~ 

— 2136192* 03 

■ 112 12E  * 0 3 

.1E183E*0S 

.lfefclBE*03 

.l*58bE*03 

« 11  8?bE  * 03 

• 99b  COE* 0 2 

• 10  69bE  *03 

• 8275b E* 0 2 

• l % W *0J 

• 15  72Jt ♦ U3 

• 145B3E* 03 

•lb380£*03 

.15562E*03 

.lbb0b£*03 

• 1057lE*03 

• 7 8885  f *02 

• 1 *•  1 35E  *03  -- 

.197982*03 

Jl*b51E*03 

• 10  988E  * 03 

.8208bE*02 

F GR  SUBSURFACE 

27 

.17602E*03 

.?0  6b0t*03 

.35  065E ♦OS 

• 194MC*  0 3 

. 17866E ♦ 03 

.213J9E*03 

.21790E*03 

.193792*03 

• 167lbf  * 0 3 



V2?322E>03 

V73l53£»0T  — 

“ V20961E»03 

- n7963E*03 

■189?9E*0J 

.23  3S6E*  03 

.2b568E*®3 

,2256tt»03 

* 1893bE*0  3 

• zaor  Jt ♦ o s 

•353%*C*TJ  3 

FOR  SUBSURFACE 

28 

.l«081f*03  — 

• 17B83T*t3 

— n7«wt»fl — 

• 138T9E*0S 

— ni212f*03 

.151B3E *03 

.16619E*03 

.159062*01 

.UIT««IS 

■ 95bOOE*B2 

nrortm — 

.imimr 

1 r0696E*03  * 

• B?  7882*07  “ 

ib380E*03 


15S6  2E*03 


imueoi 


1I571E*03 


Minnii 


FOR  SUBSURFACE  30 


.1760?F*03 

.zr6<.r£»c3 

.zurieE*o3 

180  686*03 

“ .156606.03 

.17886E«03 

.213396*03 

.217986*03 

.193796*03 

. 167 146*  0 3 

•IB JBBF^UJ 

.189296*03 

.233566*03 

.21.566  6*03 

.225666*03 

.189366*03 

.1S346E*03 


.112126*03 

.951.  006*  0 Z 

•87754£*02 

.788856*02 

.820  646*02 


• 151.68  6*0  3 

.167146*03 

.18380£*03  .273226*03 .-231636*1)3 rtttltE»t3  ~ .179636*03 

•169296*03  .2  3 3566*03  .245666*03  . 22568E*03  .189346*03 


-*306436*83  .Z30rsf*03 .743*66*  03 ,ZZZi«»M .183446*0  3 


FOR  SUBSURF ACE 

34* 

»17W3E^HJ  •lrv'l 

IOC*U3  o 13D. 

•151B3F  *03  .165196*03  .155866*03  . 11874E*03  .954016*02 

.827546*9  2 

.r*««iE*oz 

.141356  *03 .157566*93  rH.»51E»T>3 .10  9886  * 03  -6870646*0? 


.145576*03 .137786*03  HWIMI .10694E*03  — 

• 143"Gf*C3  .155626*  C 3 .144046*03  .105716*03 


• 1 B 6*»3E  ♦ 03 

.230736*03 

.745666*03 

.772576*03 

FOR  SUBSURFACE 

31 

.160*16*03 

.178836*03 

.170986*03 

.136196*03 

« 151B3E ♦ 03 

.166196 *03 

.155866*03 

. 116746*03 

.145576*03 

~ .157286*03 

.145636*03 

.10  6946*03 

.163606*03 

.15562E  * 03 

.144046*03 

. 10571E*  03 

•l*13Sf *03 

.157566*03 

.344516*03 

- — .109066*03 

FOR  SUBSURFACE 

33 

.176026*03 

.?3^%0E403 

.707376*03 

.180686*03 

•17B86F401 

.21 3396*03 

.717986*03 

. 193796*03 
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Pressure  Distribution 


FCR  SU4SUSF8C6  24 


-.10193E»01 

-.185405*01 

- • 290556*0 1 

-.36  553E*C1 

- «16660E*0 1 

Hinnm 

-.301306*011 

-.5*5576*00 

- .CIT06  0EMI0 

-.116336*00 

••215502*00 

334752*00 

- • 38  244E+  0 0 

- . 21574£*0  0 

■.wnw-n- 

• *1 0783C»OO 

-•“.156016  *00 

* E* 01 

• 27  048  E ♦ 00 

.<*5  0V5E»O0 

•666232*00 

• 79597E*00 

•4l766E*00 

— re*  sv8su**>ce  23 


.123322*01 

.207396*01 

.306706*01 

.360606*01 

.214466*01 

— •7%W3e*0«- — 

— .371406  * 00 

—.012216*00- 

;H«636  * 00 

.359  776  * 00 

.771556-01 

.1107*6*00 

.144196*00 

.178626*00 

•42484E-01 

-.293016-0? 

-.113236-01 

- • 10274E- 0 t 

• 9? 

-.255242*00 

-.450146*00 

-.665766*00 

- .824486*  00 

- • 4 1975  E*  0 0 

for  subsurface 

27  ' * 

- .22063E^QJ 

-.660356*00 

-.126276*01 

- .179036*01 

-.107216*01 

-.242196-01 

•.9*38SE-tTX 

•-.22  3056*00  - 

* .330  376  *00 

-.392456  * 00 

-.146656-01 

-.60-666-01 

-.153116*00 

-.22  7446  * 00 

-.42499EMJ0 

.654492 -Or  

-•.194566-02 — — 

--.51913C-01 

-.992336-01 

-.392496*00 

•71613E-01 

.179536*00 

.261456*00 

.329196*00 

-.220286-01 

PUR 

28 

.26  237E*00 

.755156*00 

.132746*01 

.176906*01 

.107346*01 

.58637T-01 

--.1TS1 26*00 

1100706*00 

• 2%  2 56E  ♦ 0 C 

.297586*00 

.206456-01 

.404756-01 

.376366-01 

.317676-01 

.217746*00 

.614686-02 

-.363556-02 

- . 267816-0 1 - 

-.60  2126-01 

.196156*70 

-.507832-01 

-.  156626  *00 

-.319316*00 

-.44820E*00 

-.149596*00 
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3< 


•562032*00 

.512962*00 

•36899E»00 

• 16555E*  0 0 

• 16  3 22E-01 

•03755S*01 

..2i5,i-«i 

. 289J0r»0I 

— .995532-8 J 

•66417E-01 

•S9976E-01 

• 6239  0E-01 

. 16  581E-01 

-.316052-02 

— .tbomt-oi 

— .rrrB9r-tn — 

;TS9T3E-ai — 

.32  6 79E-02 

-.1O3S6E-01 

-.168662*00 

-. 12  36  OE*  0 0 

-.85126E-01 

-.389272-01 

-.149402-01 

31 

-.b7996E*00 

-.S7177E*00 

-.39208E*00 

-.17259E*00 

-.107412-01 

*«13386t*0(r  — 

*.ir27*E*ra — 

-.28  222E-01 

— -.262232-0  2 

-.62657E-01 

-.30662E-01 

-.18796E-01 

-.895822-02 

•20696E-02 

*.tOT?JT-at 

•.mjjt'tr 

• .67929C-03 — 

-*-.tlt78E-0? 

.614492-02 

.1%1?1E>00 

• 1269  2c  * 00 

.0600SE-01 

. 365852-0 1 

.653532-02 

33 

-.23&nJE*90 

-.S83SBE*00 

-.12»39E*01 

169952*  01 

• *79798  2*0  0 

-.43758E-0r~ — 

•nn«e«ti — 

~*.TB396E*8tr' 

— *.I633B£*B0 

.106002*00 

-.37067E-91 

-.9J03SE-01 

• • 1 39  26E  ♦ 0 0 

-.  13860E*00 

.206072*00 

~ *.I7294F*0I 

-.»I69Tr-0I — 

•*.626872*83- 

350662-82 

- • 315792*0  0 

.50636E-01 

.16732E*00 

.29965E*00 

• 931  BSE*  00 

.625332*00 

FOR  SU8SUPE8CE 

34 

•27039E*00 

•76693E*00 

.13879E*0l 

. 16b64E*01 

.105252*01 

.5C350F-01" 

.13Z53T4ff0 

•267292*00 

.345032*00 

.59768E-01 

•13853E-01 

• 39S66E-91 

•878CCE-11 

.13789E*JC 

-. 133192*70 

.66585E«ar - 

•30078E-B2 

• 2T352F-81  ~ 

• 477  B7E-01 

* *200 28E*  0 0 

-.632U1E-01 

-.l3356E»0C 

-.28P65E*C0 

- .36l7CE*0B 

-.263622*10 
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SAMPLE  CASE  3: 


BO- 105  FUSELAGE  (NO  PYLON,  NO  ROTOR) 
WITH  EFFECTS  OF  SEPARATION 


Data  Input 


3 

PROGRAM  SHAPES 

HUi3  PYLON  FLOW  S£PAkATIDM 


fuselage 

- WITH  3S 

FAR4TIGN  e 

FFECTS 

KREAD 

0 

* 

12  a 

0 

* * 

0 1 

1 1 

0 0 

0 0 

L 1 

UAHOH 

. 2 04 

0.000 

- 0 . u 0 0 

ALFA 

0.003 

ij  . C 0 0 

N300Y 

a.  o jo 

12.  0 0 j 

36. 000 

41.00  0 

N300Y 

0.000 

J.  000 

„ . r U 

■ ■ . CO  c 

N800Y 

0.00] 

0.  000 

0.  000 

0.  00  0 

N30DY 

7.  25  J 

o.  25  C 

7 . 2C  0 

6. 25  3 

NWAKE 

1 0 

0 0 

NPYLON 

1 1 

CSTAG 

-2.. 

1 • w J 

KPkINT 

1 1 

1 1 

0 1 

1 1 

N300Y1 

4 -» 

2 1 

1 

N80U Y2 

4 4 

1 1 

1 

N80UY3 

•4  4 

2 1 

1 
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Specifications  of  the  Problem 


t 


..  FOR  PAP.T  . i_.  „ 

MX*  4 
•iY*.  L 

FOR  FART 
NX*  4 

2 

NY*  u 

...  fCA^fARI 

, 

NX*  A 

NY*  4 

. . .. 

FOR  PART 
NX*  A 

NY*  « 

5 

FCr.  PA^T 

NX*  4 

NY*  4 

c 

FCR  FART 

NX*  4 

NY*  « 

7 

FOR  FAF7 
NX*  4 

NY*  4 

R 

FC»  PAr*T 
NX*  4 

NY*  4 

10 

FOR  CAPT_ 
NX*  4 

NY*  .4 

11_ 

FOR  PART 
. NX*  4 

13 

NY*  4 

FOR.  PART 
NX*  4 

_Nl*  4.  ... 

14 

FOR  PA  FT 
NX*  4 

NY*  4 

15 

NEL£-*1}2 

<synmy*  o 

<$YM*Z*  P 
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X,  Y,  and  Z Coordinates  of  the  Nodes 


fiOCx  * 

3* 

» 

2 

1 

c . ■:  3 0 : 3 

3.  j 0 J 3 0 _ 

0 . C 1 a a a „ _ 

2 

• 7 e 0 3 3 

-. 0C0-G 

2.06250 

3 

3*00000 

* • 0 0 j . c 

-.  12500 

4 

f .79003 

-.ooojc 

1.  Ia2**0 

5 

12.00000 

. 30  3 Q& 

6.2:000 

& 

.75  000 

. a9JM 

v.canKo 

r 

3.0000) 

1.  J6723 

3.61100 

9 

c.75  003 

’.0305. 

5.21630 

9 

12. 00  DO 3 

.....  i. 

2..  6 ? ?.91 

ic 

*f9?o: 

l.iilM 

1.45*41 

u 

3.00 300 

2 .56326 



12 

o.  75  0CJ 

3 • 3 V * o 9 

-.37522 

jl3_ 

__12.»_00.P  09 

5 . 126£2_ 

5. 83363 

14 

.79000 

1.37453 

.76923 

15 

3.00000 

3.3*936 

It  £735  7 

16 

6.75000 

5.02J59 

2.367f5 

17 

12.00003 

a. 39313 

3.15  714 

16 

.75000 

l .81250 

0. 00000 

19__ 

.3, 00  0 03. 

3.325JQ- 

a .JULIAS 

<J 

6 . 7 5 m Z 

a.42750 

0 C 2 

.1  _ 

12,00000 

r.is ooo. 

i..oo.oAa_  _ 

J2 

16.50030 

-.00030 

6.25000 

13 

2? . oonoa 

- . 03000 

5. 25000 

ih 

31.53  000 

-.36308 

3.25Q0Q 

....  . <5  _ 

_ j <■ . 

...  . ^ .... 

. 

26 

1 c • '»  ) C 0 3 

i.rrthi 

7.1 2201 

27 

>c .0000) 

2.72.-!- 

7 . » 2201 

2c 

31 . 33  LQ0 

2.27445 

/.I  2201 

29 

J6,0uC9) 

2.77--S 

. .'.12201 

30 

1 fc  . 3 0 0 ) J 

..  . 12  i . 02 

5.  cllt.3 

31 

21.000:) 

5. 126?  2. 

5t_C.lTc3 

32 

31.50:03 

0.12a:2 

5.833*3 

33 

Sr . 0000) 

a.l 20. 2 

_ .3.  c 3363 

3s 

18. 5009) 

a . a 9 3 1 3 

3. 15 '1- 

3J 

2J.900Q) 

5. j’atS 

3,1:  715  . 

36 

31.5000) 

a • 39a 13 

3.  16  71. 

. . .37 

36.00000. 

a . a Ol  1 JL 

J..JCJJL- 

3a 

lu .50030 

7 .25000 

0. GO  000 

39 

2 r . 0 : " ' 

2.2 

« J. 

40 

31 • SO  C 9 3 

2 .25000 

3. 0Q0G0 

<*1 

J6.3GCQ) 

.'.2-000 

- a.aQocj 

42 

39.31293 

-.00300 

6. 1875C 

s3 

. * 0 . 25  0 0 a 

_ -.acooi. 

>OZ5 

*4 

4 0.31 2 5 J 

-.00300 

7. 047*) 

-5 

>*1.00  300 

0 . 0 L 0 0 G 

_ a.Q99QC 

~6 

39.  31  253 

2.08334 

c.7i*5o 

47 

4G.7JC9) 

1.3 1023 

i^ilXQa  — 

46 

4c . U75-. 

. 193  ,1 

1 , c ; 55 

. 49.. 

39. 31253 

3 • a 4 4 : -i. 

it.  *T"*99 

50 

st.  25C0J 

2 . Of  126 

7.  91*- 67 

51 

40.51290 

1.2:163 

_ U4524J— 

52 

39.3125) 

a .02353 

7. 39743 

S3 

41 .2i 093 

3.3433S 

54 

« J . 3 1 25  ) 

1 • 6 2*  a 3 

• 76978 

- 55 

J9*.  31.25 1. 

n^nonon 

56 

4b. 25  000 



3. 62a  00 

0.00000 

57 

4C .3122- 

. 1.J12J- 

^.—12  0 •*  ' 

56 

• 7 5 Q 9 0 

-.00000 

-2.0525  0 ~ 

19.— 

3.3Q20 3 

-.10  0,11. 

-..125110 

(0 

6.75033 

-.00000 

-fi.  13  7*0 

tl 

„lZ.Q0fliU 

«.aco  oil. 

-a. 66. i a: 

12 

.73  003 

.33331 

-1.90550 

*3 

3. Q3C00 

1.36723. 

-3*iiiAa. ._ 

14 

c • * 3 0C3 

2.0808- 

-5.71650 

f5 

12.00003 

2.774-: 

-7.fa2  201 

04 

./;  :: 

1.2*1  .3 

-1.4.0  361 

17 

3.00300 

. 2.ia32i 

- ,r2Jali02_  _ 

(0 

- .75  00  9 

3. a-. >9 

--. 3’*  22 
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THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 


i'i 
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1 5 

-.701672*00 

-.M-25£»  30 

-.7T709E*00 

- . 8J *2J£  * 00 

-.2?233-»::  _ 

- .2551/£^GD 

-.695186-01 

-.72807E-01 

-.763966-01 

- . 75  37  9£-  0 1 

^z.XlMSi QiX 

_-juLL.901i.-i.l-- 

-U.:i27---i_L... ... 

.i^93  713£r(U— 
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AEROSPACE  SYSTEMS  INC  BURLINGTON  MASS  F/G  20/4 
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ASI-TR-76-38-V0L-2 


USARTL-TR-7B-1B 
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SAMPLE  CASE  4: 


MODEL  1 HELICOPTER  CONFIGURATION 
WITH  ROTOR  WAKE  AND  SEPARATION  EFFECTS 


Data  Input 


l 

o°0»RA,'l. SHAPES  .... 


HUB  •’YL'IN  FLTH  SEPARATION 


FUSELAGi 

- SINGLE 

ROTOR  (TWO 

9LA0ES) 

KR£\0 

♦ 

0 

16  0 

0 

** 

**» 

(1  1 
0 ..  0 

1 1 

0 0 

0 0 

1 

IJHACR 
ALFA  . 

.067 

0.000 

3600.000 

0.000 

-0.000 

NRO  IT 
N900Y 

0.000 

0.000 

9.000 
0.000  . 

18.000 
0.J10  0. 

36.000 

0.000 

N90TY 
_NS09T.  .. 

0.000 

.1.450 

0.000 

..  3.250. 

0.000 

1.450 

0.000 

2. 900. 

ROTOR 
_LH  El  7 5 . . 

20  000 
0.000 

S.  040 
0.000 

11.550 

0 .0  00 

15.050 
0.000  . 

0.000  0.000  .12) 
90.000 

RPITCH 

K9HQE 

0.000 
? 1 

0.000 

NHA<E 
_SPI  RAJ.  ... 

1 0 

0 1 
00.Q_.450.0  0 0. 

NPYLON 
CSIAG  ._ 

1 1 
. -60.000 

. . 1.000 

KPRINT 

NBOJY 

1 1 

3_  3 

1 1 

2 1 _ 

0 1 
1 

1 1 

1 1 

NBOOY 

.NBOOY 

3 3 

3 . 3 

1 1 

2 . . 1 

1 

1 

NPOTnR 

3 3 

2 1 

1 
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Specifications  of  the  Problem 


por  pari.  1 

NX  = 3 
NY=  3 


PC”  PARI  ? 
NX=  3 
NY=  3 


FOX  PARI  3 
NX  = •* 

NY=  3 


^0°  PA»I  5 
NX=  3 
NY=  7 


FOR  FAPI  6 
NX=  ’ 

NY=  3 


FO°  P A P I 7 

NX  = S . _ 

NY  = 3 


FOR  P AP  I 9 

NX=  3 
NY=  3 


POR  PARI  IT 
NX=  3 
NY  ■=  7 


FOR  PARI  11 
NX  = 7 
NY=  3 


FOR  PAPI  13 

tlX=..3 ... 

NY=  3 


FOR  PARI  14 
NX=  3 
NY=  3 


FCR  FART  IP 
NX=  3 
NY=  3 


FOR  PART  24 
NX  = 3 
NY=  3 


FOP  PAPl  25 
NX=  3 
NY=  7 


FOR  PARI  27 
NX  = 3 
NY=  3 


Fop  PA° I ?P 
NX=  3 
NY=  3 


NELr1=l4^ 

KFYIMyr  0 
KSY  H 7=  o 
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X,  Y,  and  Z Coordinates  of  the  Nodes 


NOOE- 

X 

Y 

Z 

i 

3.00000 

0.0003  0 

0,00000 

2 

1.10000 

-.00000 

1.00333 

3 

■^TTTTB 

-.00009 

4 

3.00000 

.00000 

3.25000 

« » 

KWlFTH 

■\TfTTI 

6 

6.30000 

.60339 

1.87639 

r 

9.00000 

■nmn 

6 

l.aooio 

.61058 

.56167 

9 

4.00009 

■imm 

10 

9.00000 

1.25576 

1.62500 

11 

^bWTTTT^H 

■nm 

-.30000 

12 

6.00000 

.96667 

-. 00000 

13 

3.00000 

■nuTTl 

s 1 1 ] *.55 

14 

12.00000 

-.00000 

3.13333 

15 

■InilriB 

-.00000 

3.01667 

16 

19.00000 

.00000 

2.90000 

. 1* 

u*nn  1 1 ■ 

^■YTXTT 

IS 

15.00000 

.72500 

2.61251 

19 

RTRTim 

Mill  1 1 ■ 

■Aims 

2* 

12.09000 

1.25576 

1.56667 

. 21 

■in  I II T ■ 

■WrTTB 

22 

18.00000 

1.25576 

1.65008 

23 

mwAiim 

-.00000 

26 

15.00000 

1.65000 

-.00000 

26 

■irhmbi 

1.65000 

-.00080 

26 

29.00000 

-.  00000 

1.93333 

17 

tv  ^ m 

-.00900 

?0 

36.00000 

0.00000 

0.00000 

29 

28.00000 

Ml  III  ■ 

30 

36.00000 

.26167 

.03716 

11 

29.00000 

WTtTTa 

32 

36.00000 

.61650 

.60333 

31 

28.00000 

•96667 

-.00000 

34 

9 6.00000 

.68333 

-.80000 

35 

1.90000 

-.  00000 

- 1.  06333 

36 

4.00000 

-.00000 

-2.16667 

37 

hr  - - m 

■II 1 1 1 1 ■ 

-3.25000 

38 

1.00000 

.26167 

-.93019 

39 

6.00000 

.68333 

-1.8/639 

40 

9.30009 

.72  509 

-2.01650 

41 

■TWTTTt™ 

^1  ' 1 1 1 ■ 

-.56167 

62 

6.30000 

.93716 

-1.90333 

63 

3.00000 

■ IfUU* 

-1.62500 

44 

12.00000 

00008 

-3.13333 

65 

15.00000 

-.00000 

-3.01667 

46 

10.00000 

.00000 

-2.90000 

47 

■In  II  IT  ■ 

Ml  HUM 

48 

15.00000 

.72508 

•2.61251 
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1 


98 

36.00000 

-.41050 

-.40833 

li.S5  0flfl 

j m 7300  ft 

6.67000 

100 

11*33009 

3.75000 

6.70662 

101 

13^10336 

3.73000 

6.75662 

10  2 

19*15000 

3.75000 

6.67000 

1Q3 

11.33000 

12.77778 

6.67000 

10  V 

It* 93009 

12.77778 

6.70462 

iQ3 

13.1  ft  336 

12.77770 

6.73642 

136 

15.05000 

12.77778 

6.67000 

iot 

11.33000 

10. 196k  % 

6.6700ft 

103 

11.93009 

16.19%%% 

6.70662 

ioq 

1 3.1  0336 

10. 19%%% 

6.75647 

110 

15.05010 

10.19666 

6.67000 

ti  i 

11.33000 

70. 00000 

6.67000 

112 

11.43009 

20.00000 

6.67000 

in 

13.10336 

70.  00000 

6.67000 

11 V 

15.05000 

20.00000 

6.67000 

m 

11.13003 

3. 75000 

6.63538 

116 

13.10556 

3.75000 

6.58358 

If  7 

1 t .330  0 Q 

17.22778 

6.63538 

11  0 

1 3.10556 

12.77778 

8.58358 

11-9 

11.13003 

18.  19%%% 

6.635  38 

120 

13.10556 

10.19%%% 

6.58358 

I2t 

,n  .33HOQ 

-3.75008 

6.67800 

12? 

-11.93009 

-3.75000 

6.70662 

i?f 

• 1 3.1  A 336 

-3.75600 

6.75667 

12% 

-15.05006 

-J.  25800 

6.67010 

t Z3 

-11.55000 

-12.22728 

6.67000 

126 

•11.61809 

-12.27278 

6.70462 

i>7 

-11. f 0556 

-1 7.7777  ft 

6.75642 

12  8 

-15.05000 

-12.77778 

6.67000 

igq 

-11.33000 

-18.19666 

6.67000 

130 

-tl. 13009 

-18.19644 

6.70462 

131 

-13.10556 

-18.19646 

6.73642 

13? 

-15.05000 

-18.19644 

6.67000 

13  3 

-11.35000 

-28.  ooaaa 

6.67000 

13% 

-11.93009 

•28.08800 

6.67000 

133 

-13.10556 

-20.00000 

6.67000 

136 

-15.05000 

-20.0000  0 

6.67000 

13  7 

- 11.33009 

-3.75000 

6.63538 

1*8 

•13.10556 

-3.75000 

6.50350 

133 

-11.93689 

-12.77770 

6.63530 

l%0 

-11. 10556 

-12.77778 

6.50350 

1%1 

-11.93089 

•10.19446 

6.63530 

1%? 

•13.10556 

-10.1946  6 

6.50356 
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X,  Y,  and  Z Coordinates  of  the  Centroid  of  Aerodynamic  Panels 


YPC 

z°c 

1 - 

,5005S 

. Q6Q42 

mcrnrm 

2 

2.5)277 

.10125 

1.51615 

1 

ft. *0714 

— Ill  ■ 

% 

.S00S5 

.16506 

. 36997 

5 

■jfTTTTfl 

■vTTTtb 

6 

6.51719 

.42531 

1.86983 

7 

■nTTYl 

9 

2 > 1 277 

.67663 

.60625 

9 

mi  mil 

10 

1 0.51162 

.36250 

2.97707 

11 

1 1.41444 

1 111 

2m A&901 

12 

16.51029 

. 36250 

2.76016 

13 

10.51162 

Ml  IIIIB 

14 

1 3.51696 

.9903  7 

2.10  026 

IS 

rr  *r  •- , v 

.99037 

^ 1 A 

16 

1 0.51162 

1.35207 

.79792 

17 

13.51696 

■WTTTT1 

Ml  'll  ■ 

10 

16.51029 

1.35207 

.73990 

14 

rnrrnr* 

M0  Mill 

20 

31.03630 

.10125 

1.35287 

71 

39. li  072 

■rrYl 

i^LVI  6 ' M 

22 

23.42969 

.82931 

1.65061 

23 

rnrnv* 

.49S16 

.94037 

26 

35.13072 

.16506 

.33012 

24 

73. 12964 

26 

31.43630 

. 67663 

. 36250 

2/ 

S ' ~ 

^mTTTTI 

20 

.5)095 

.06062 

-.50530 

24 

■fTTTTl 

-X.  51615 

30 

6.54719 

.30200 

-2.52691 

11 

■firm 

M 

32 

2.50277 

.69910 

-1.10990 

11 

6.50  719 

Ml  1 1 I 1 ■ 

BBirii 

36 

.50 059 

.22560 

-.13562 

14 

mT7tT1 

via 

36 

6.50719 

1.12739 

-.67708 

17 

1 Q. SI  142 

^■ITTTT® 

-2. 977A7 

30 

13.51696 

.36250 

•2.06901 

14 

Ml  Mill 

-2.  76  016 

60 

10.51162 

.99037 

-2.17995 

41 

13  91696 

.99037 

bhi  rui 

62 

16.»1029 

.99037 

-2.42050 

63 

1 0.51162 

1.35207 

-.79792 

*4 

13.51696 

1.35257 

-.76075 

4§ 

lft.41A>4 

Mtzrm 

46 

23. 42  965 

.30208 

-2.25670 

4 7 

11.144*11 

- t A 1 24 

BmtITI 

46 

35.03072 

.06062 

-.65096 

44 

■OTTXTTH 

■ -1.6806X. 
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■50  310J5I8 


59518 


9903* 


[ 


5 2 

28.8  2565 

1.12739 

-.60617 

S3 

31.03610 

.67663 

— * 36960 

56 

35.03 S72 

•22568 

-.  12063 

ss 

.«hik 

.60  61/1 

56 

2.50  277 

-.18125 

1.51615 

57 

6.  sana 

-.30906 

9 . 69  691 

5* 

.50055 

-.16506 

• 36997 

<0 

2*  SB  777 

-.69516 

1.10990 

60 

6.50  719 

*.82531 

1.86983 

61 

**0066 

-.22. 9. 

.13669 

62 

2.50277 

-.67663 

•60625 

mt 

6*iU719 

-1.19739 

.67788 

66 

10.51162 

-.36250 

2.97707 

M 

U.itMk 

— .3696  6 

2.86981 

66 

16.51825 

-.36258 

2.76016 

67 

10-61 16? 

-.99837 

9.17999 

68 

13  »t  696 

-.99037 

2.10026 

66 

16**1  *75 

-.99837 

9.89068 

70 

10.5116 2 

-1.35257 

.79792 

71 

1 S..lMk 

-1.3S987 

.76  876 

72 

16.51875 

-1.35282 

. 73958 

73 

93*0  9 666 

-.30288 

2.751.78 

76 

31.J3610 

-.18125 

1.15287 

7% 

16*1  3*77 

-.0601.9 

. 66096 

76 

23.02565 

-.82531 

1.65  061 

77 

3i  i iu  in 

69si i 

.99037 

78 

15.13872 

-.16506 

.13012 

79 

93.U9S6S 

-1.19719 

. 60  61  7 

88 

31.036  30 

-.67663 

.36250 

/ii 

IS.  1.11179 

..  99666 

. 19  083  __ 

87 

•5085S 

-.06062 

-.50530 

63 

_2**fl  9 77 

-. 1619S 

-1*51616 

St 

6.5U  719 

-.30208 

-2.52691 

*s 

-*0*66 

-.16606 

-*36997 

86 

2. >0277 

-.69518 

-1.10990 

07 

fc.iil7^ 

-.89631 

-1*86981 

88 

• 50059 

-.22568 

-.13562 

?.>ll  977 

-.676^3 

_ -.68696 

98 

6.59  719 

-1.  12739 

-.  67  708 

•9.07707 

92 

13.51696 

-.36258 

-2.06901 

nt 

-.16968 

-2. 76016 

99 

10.51162 

-.99837 

•2.17995 

-.99037 

.7* 16096 

96 

16.91825 

-.99037 

-2.02058 

qy 

1 «-•»<  (19 

•1.369*7 

•*79709 

96 

IS.  *1999 

-1.35287 

*.75875 

141  2 3.12545 


30206  *2.29473 


tai — 31*134.10 =*16125 -1*39267 


102  39.13872 

— UL3 73.47949, 


-.08042 

.-*07931 


-.49  098 
•1-83  i)8t 


104  31.03430 

163 33.03372- 


-.49918 

=*-16566 


99037 

=*33012 


108  23.02949 

167 31.03436 


-1.12739 

-.67643 


-.60417 
— 16730 


106  39.03872 

109.  11.79744 

110  12.93608 

111  14. 19313 


-.22940 

0.26309 

0.26389 

6*26369 


-.  12  083 
6.68731 
6.73092 
6.71321 


112  11.79744 

113  12.33686 


19.46611 

15.46611 


6.66731 
6.  73052 


114  14.09339 

115  11.7J744 


19.48611 

19. 8972 2 


6.71321 

6.67666 


116  12.93600 

117  14.14339 


19.09722 

19.09722 


6.70  026 
6. 69161 


118  11.79744 

—119 12*43616- 


6.26369 

6.26369 


6.89269 

8*611943 


120  14.19  335 

121  U. 75  746. 

122  12.93600 

-123 14*09316 


0.26309 

19*44611 


19.46611 

-19*61611 


6.  82  679 
6.65269 


6.68  948 
6.62  679 


124  11.79744 

L25.  12.13611 


19.69722 

19.69722 


6.66134 

4-63974 


126  14.1*339 

127  -1  1.78  744 
126. -12. *3608 
129  -14.0  J139 


19.69722 

-1.26319 


-8.26389 

-.1.26319 


6.  64639 

6.68731 


6.73  092 

6.71321 


130  -11.79744 

131  -12.5  3616 


-19.46611 

•15.46611 


6.68731 

6*72052 


132  -14.09439 

L31.  -11.75744- 

114  -12.93606 
139  -14.01119 


•19.40611 

-19.09722 


-19.09722 
-19. 09722 


6.71321 

6..67066 


6.70  026 
6.69161 


136  -11.79744 

137  -12.23410 


136  -14. "9135 

133-711.75744- 


-0.26309 

_7.0l*263A1_ 


•0.26399 

-19.46611 


6.69269 
-6*613.41- 
6.  62679 
6*65269 


140  -12  93608 

141  -14.09  318 


15.46611 

-=15*41611— 


142  -11.73744  -19.09722 


6.  60949 
-6.62679 
6.66134 


144  -14.09  315  -19.097  2 2 6.  64  0 39 
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Nodal  Numbering  for  Surfaces 


102  106  110  116 

FOR  SURFACE  25 

"Ti  Tor  ioT  uT 


126  128  132  136 

'OR  SURFACE  28 
Til  125  129  13? 


The  Distribution  of  the  Source 


Ti,,riE*02  -.rt439E»02  -.7i»2S2S»02 


,liriw*gt  .2KWC»I1  .3S99tE»«l 


-.J9793E*92  -.a9i»29E»0Z  -.890ESE*02 


•.  72382EM2  -.M991E»02  -.U299E»02 


929»2E*82  -.*kl(3E»lt  -.MM«H2 


-• 18/8SC»«2  -,19I92E»82  -.18739E»02 
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The  Distribution  of  the  Velocity  Potential 


EOF  SJ8SURF»C£_ 

1 

-.4609  92*0  2 

-.450376*02 

-.456426*02 

- . 4 3515E*  02 

43032E*  02 

- 3»2oae*02 

393146*  02 

-.399312*02 

FtlR  iihsuofoee 2. 


-.2075  16*02 

-,23836E*02 

-. 25503E*  02 

- qfc&4?g»Mi 

--112126*02 

-■  160517*  02 

.594796*01 

•50202E**1 

.196632*11 



- 110161*03 

-.10755E*83 

01939E+  02 

-sai/2F*ai sUialhll 


-.  207 12E*0  3 

-.264946*03 

-.23479E*03 

FS9  >U*l5U?rACP 



-.464492*32 

-.461332*02 

-.457316*02 

-.441076*02 

-.40J63E*02  -.401516*02  -.40099E*02 


»1ii  ^UBSInriCE — 6. 


-•  241286*0  2 -.255256*02  -.250936*02 


aJMMUtlU  --201156*02 


-.179736*0  2 

-.14600E*02 

-.105026*02 

Fn»  SJ2SURF1C6 

7 _ 

-.460172*02 

-.919466*02 

-.60  3326*  92 

-.144906*01 

. 174396*03 


197236*03 


!I7ME»I3 


- ia«Z»'*OZ  •.166trE»*Z  -.!Ji.ZtE*OZ 


. 4MJ6»»0Z  -.SJ5fJE»0Z  '.6?SarE*0Z 


1910HE+93  -.3  66Q«E*-03  31439E»03 


SJ9SIHFECE  1 


%8b0l£-9t  -.68291E-01  -.10»»VE*«0 


.»MtE»ll  .MNIIXI 


.ujme«ei  .riiwMi  .3«ntE*oa 
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FOR  SUBSURFACE 

» 

. 17S?6E»Q  1 

■lgg71E*01 

• 4533BE+00 

-.SB901E-01 

-.gSa>SE-01 

-.037S9F.nl 

cOR  SUBSURFACE 

10 

-.  I6»a3£*ua 

-.Hg73F*«H 

-.  14UJ0EM1O 

13943E»B0 

■az»UE*iio >577hse*qq 


i 


FOR  SUBSUBFACE 

11 

■ 91«96£*H0 

j7S143E*  OIL 

. rzrrzfio 

• 0617OE*  00 

15B14C*B1 

.01Z47E-01 

FOR  subsubfacf 

13 

amaEtno »i2334E4bi 


’•  4411  0E-01 

-.5E030E-01 

- 741413-01 

-■  7Sg31F-01 

'OR  SUBSURFACE 

14 

-.11H/9F4A0 

•.unartOD 

ajioee-oi 


AuaaaEt.Bi. 


-.55051E-61 


--T7ilE:-m 


TU?uiE*ao 


- llft4«E*#0 


1 **77931  ♦ 00 


.S664»E*00 


■ lSOSOF-HO 


■.U6MtE»01 


-• 5E503E-O1 


.r-Zjt/OlErQl 


-.unitm 


77U7E-B1 


FO*  SU99IRFACE  U 


>161«tE»00  ,Z07i4E«-0Q  .29630E»n0 


. 3 7bZ9E*0  0 .33SISE»03  .3J673E*08 


-.377ZHE+03  -.837?9E»D3  -.87S#4i*03 


I 


SECTION  6 

PROGRAM  LISTING 


PROGRAM  MXN(  INPUT. OUTRJT  , TAPEMINPUT,  TAPEMOUTPUT) 


COMMON/ ZZZ3/(PTL0N,N80Wl'»  NBOOTZ.NBOOV  3,  NPTAIL.NSH4FT, MMUB.NSHANX, 

a u mm  Afiff 


COMHON/ZZZA/UHACH, OMEGA,  ALFA, ABCTA,  AREA 


COMMON/  ZZZ6/XN03E,  ESDI  ,(  3 Of,  XT  AIL 


[*T.  .(•rT#*>4/rrTiV' jn  iimli'i 


COMMON/  ZZZ8/ZNOSE  , ZB  Ot ,?  9 0*  , ZTAIL 


VbbbATtA  l T|TP 


CONMON/ZZZIO/  RSHAFT,  LSHArT  ,RSHANK,LSMAM( 


CONHOH/ZZZ12/RROTOR,BCKi*0,TAU8L,AirA3 


• T ' 111  TbbbAb.  'A  4 ' * bA 


COHMON/ZZZ1A/K8LAOE, TANL i B > T ANTEB.XBLt ,X8TE,KR0T0RS ( 3A) 


riT'F###lLfJ'LT  T CM  Jlvi1! 


C OHMON/ ZZZ 16/  NN  A HP  T , NM  K i B > NMA  <S<  , N MA  K 9L 


C0HN0N/ZZZ10/HANGPY,MAIG(B,MANGS<,NANG8L 


■H  •*  \T*T  F#  t r J 1^1  I J li  I llil 


C OHMON/ ZZZZB/ PI 


COMMON/ ZZZZ2/KNSELE?  <NS3 (P.KNSTTP 


» b 0 A m Vi  J-Ti  (■ 


C0MH0N/Z7Z2VKTNELE,  KTN3(P,<TNTVP 


■Mil.  I.  ill  ri  CCIJ 


COHN ON /ZZZ26/KV TELE,  KVRIP.RVTTTP 


COMHON/ZZZZS/KHBEIE.XMBXNP.KMSTYP 


Tb  0 Al'.'Ji-n  11^4^  rfl  I 


C OHMON/ ZZZ30/ <8t£L E , K3LS -IP  i K81TVP 


rii:  :nrr#>*j.trj”k-*/T-wrriT  . 


C 0MN0N/ZZZ3Z/  C STAG  , C VORT 


C OHMON /ZZZ3W  ▼ CUT 


COHNON/ZZZ80/XCTR,VCT4,'  JTR,RX»RY*RZ 


bbbl  < Vi’  71 


OIHCNSION  *Cl T(280l  , TEL’  (280)  ,PHII(  250  >,KROTDR(280> 


DIMENSION  PC(I,28«>,P1(3,Z80I,P2(I,  288*,PI(3,288I  ,(<(3,288) 


DIMENSION  PMIC  (288  ) , PM II  ( 280) , PHIZ  ( 280  ) , BC  (280 1 


WYHAX*6 


NVORTaO 
JIWHW 
NMAK9L  *0 


REA0CNREA0,2>  NCASE 
— tHITF  1MMBTTF, 01MCA9F 
i FORMA  T (1019) 


A FORMAT  (2A6,2I9) 

00  999  ICASE«1,NCASE 


WRITE (6,60) 


1 CALL  FRINTA(W£LEN,NXYM=,  XX,NMOOE,PC,NOFCT,NOOE,t) 

C 

IF(XPRIMTCtl.EQ.l) 

(CALI  PRTMTAmFl EH.Mif  »mp- « ».  »i>innrf»r.  rnnrr:  rf  HOUUZl 

C 

CALL  OfCl?T(l«LrM,MOnFt-.CtB1.R>»Rl,XK.IMOOEI 

c 

IttXBRIMTMl.tg.n 

1C  ALL  PAINT*  ^£LEN,NXYHP,  XX,  NNOOE,PC,NOFCT,NOOE,S> 

C 

IF  (KPRIMT(6I  .EQ.l) 

1CALL  PAINT  A fMClf  M»MXYMA» XX, MMOOE,  AC  ,MO«CT»AMOE,6) 

C 

*6fQ«»6LEMllt 

C 

an. 6 mtMi.MCTQ 

AA(NN)«0.I 

6-  CQWImii 

c 

on  a i»i,MELtM 

NNNaI*(I-l>»NELEN 

xnuarrin.A-ii 

AA(NNN)>1.0 
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9 CONTINUE 

CALL  COEFF (NELEH, PC,  NNOOI  ,XK, NONE, K MAKE,  NESQ,AA,  SOURCE  ,8C 


C 

1 kntxac.nl.  a i 

1 CALL  COEFF«NELEN,PC»NN30E,XK,NOOE,  WAKE, NESO.AA, SOURCE, 8C, 


t CALL  COEFF (NEL  EM,PC ,NN)DE, XK,NOOE,  WA  KE,  NESQ,  AA,  SOURCE ,8C, 


C 

IFlMBLAHF.EQ.1  .AM}.  MMA<  81  «NF.  8) 

t CALL  COEFF<NELEH,PC,NN  OE,XK,NQOE,  WAKE, NESO.AA,  SOURCE, 9C, 


c 


■mnaominf  oata  (nf.iemi 


C BATE  CREATED  t MARCH  31,  1979 


COMMON/ ZZZ1/3X(  SO)  ,NT  (»*  , MXY  (3%)  » KSYH  NT , KSVMNZ.NSYMM  Y ,NSYHM Z 

cnHiinM/zzzz/Mm.NsaonT.  xs.nt  (ia>  . i seacc(  ia> .emorhc  i iai,  aha  rest  tai 

C0NH0N/ZZZJ/NPYU>N»NB00rl,NB0OY2,N80QY3,N¥TAIL.NSHAFT,HHUB,NSMANK, 

1 Hflt.  A OE 

COMMON/ ZZZA/UMACM.OHEOA,  ALFA.ABETA 

C/lMMOM/777B/TBTCTa,JaTCIB.IBTCTB,)M[BYL.>TBTl.B7B¥l  - ■■ 

COMMON/  ZZZSA  NOSE, XSOl,/ 302, XTAIL 

COMMON/  ZaZdXMDSL,  tBOl  ,Y  iflE.YTAII 

COMMON/  ZZZB/ZNOSE , Z801 ,7  302 , ZT  AIL 

coMNQ>u7ZZ8/Braot.Bzaoi.a¥BO»,iu»oa 

C0MNON/ZZZia/RSHAFT,LSI«rT,RSHANA.lSHlNK 

COMMON/  777 1 1/  X MUBca,  YMUA  7 B . 7MUBCA,  B XNU  8.  ATNUB  .BZMU3 

C0MM0N/ZZZ12/RR0T0T  i BCN  2D,TAUBL,ALFA3 

COHMOM/.ZZ1.1.V  TMtTTS^IMg  1 C , INI  US^CONI  MO , A ZZMUXM 

C0MM0N/ZZZ1A/R3LA0E.TAN-  • 8 , TAmTEB,  X8LE  ,x  3TE  .KROT3RS  ( SA) 

COMMON/  17316/  XSOAM,  X06ZB  , XTEIX  .TAMOtH,  TAMTSW.tAU,  XMT  AU, 

C0MM0N/7ZZ16/MMAXPT , NMA<(B,NNAKSK«  NMAKBl 

— COMM  CM/ 7 II 17/ Ml  KOBE,  MAIP  48  , MAKOSK,  MAEO  BA 

COMMON/ ZZZ18/ NANOPV  * MAN;  MB  iMANGSK,  M ANf,  BL 

CflMM0M/.7?31.8/lta«IHl(l8>.  MRCAO.MMBITiyXBEAO 

COMMON/ ZZZ20/PI 

COMMON/ 12321/68*1,  K»T1 

C0NM0N/ZZZZ2/XNSELE,  6NS3  3P,6NSTyP 

CORMOM/777  7T/ranEl.C.XBnB4B,irmTTB 

COMMON/  Z7ZZA/  XT  NEI.E.6TN3  HP,  XTNTYP 

rnHMnM/777>M/MBYri  fr  MBV3  4P  , MPT TTD 

OQNMON/Z7Z28/6TTELC,  6*TS3P  ,6* TTY* 

COMMON/  7ZZ  77/  KSMEL6,ESH;4B,XSHTTB 

3 OHHON/ ZZZZ 8/  K HBELE,  6MB3 HP  , XHBTT® 

cnMMOHzz77Ja/r«tnr-JC«>MB.  esettb 

COMMON/ZZZSO/XSLELE.KBL'MP.KBLTyP 

rnHMnH/777Ti/M<Tir.TMtfnpr.u4PtP4i  -Ycrati 

CONMOM/ZZZ 32/  CSTAGvCVORT 

COMMON/ 7 7Z.3A/TC  III 

COMMON /ZZZ33/UMAKE 

irOMMOM/777AB/YrTB,  YCla.;:TB,ax,B«.A7 

COMNON/ZZZBl/NSPANy  XLEZ-  f XTEZN,TAMLEM,TANTEH,  TAOM * ZP 

COMMON/ 77738/  ZTNr’STygBITCM,  TMI3T 

OINENSION  AC30R0428)  ,MN’  » (2)  ,HNPK2) , MNMP(2)  ,MNMK?) 

OIMKNMOM  06(20), 6SFACK7I) 

OINENSION  TCMOROtlO)  ,T  MIS  (10 » . MCMORO  ( 10  > , VAXIS(IB) 

BEAL  I,  INtAI^lMAMM 

C 

MSTXm  A 

00  19  I>1,3 

Bflf.(HBEAn,)8) 00 

19  MRITE4  NMRITEt  20(00 

BEA01NBEA0.il> 66  <1>  . IT.  ( 2>  .XBEAO 

WRITE (NNRITE, 20) 66(1 1,60(2)  .KREAB 

JJ FOAM!  I (71 A A) 

28  FORMAT ( IX, 29AAI 
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USE  KSYHHZ»0 


AEAO C NREAO, IS ) GGC1),GJ(Z>  .NStRSTHMY, KSYNHZ 


REAOC  NREAOt  13  ) GGC1  > >GS(  Z)  tNPYLON,  NBOOT  l,NB00»2,Na0  0Y3,  NYTAIL, 

uru icy  uuiib  uouiui/  uai  a ■> m 


NRITE  CNMRXTE,2S>GGU>,GGCZ>, HP  YL0N,N8J0Y1,NB00VZ,N83  0Y3, NYTAIL, 


REAOC  NREAO, IS  > GGC1I  ,6,  C Z»  »*CPYL1,KPYL  2 


FORNATCZM.IOIB) 


liili.r  lid  ITVI  ldl 


REAOC NREAO, 1? ) GG  ID  , GG  C Z)  .UMACH.0HE3  A,  AREA 


REAOC NREAO, 12)  GG(1»  ,GG(»>  .ALFA.ASETA 


F0RHATC?A*,7FS.3> 


|T*I.  (lid  T 


IFCNPYION.EQ.!) 


IFCNPYL0N.EQ.1) 


IFCNVTAIL.EQ.  1) 


■ I (-(:#  \inl 


Z tZPV 


INRITE  CNNRITE,  22)  GG  Cl)  ,1G  CZ » , YSPAjC,  XL<£  ZV  . XTEZV,  TANCE  Y , TANTEV,  T AUV 


NROTORsNSHAFY  ♦ NHUB+NSHA.K  tNBLAOE 


IFCNTBO.NE.O) 


IFCNTBO.NE.O) 


• J ■ r r*  ( • ^ ■ 


IFCNTBO.NE.O) 


Hijliinii! 


IFCNTBO.NE.O) 


IFCNTBO.NE.O) 


II  I JiiflMTiliFlril 


IFCNTBO.NE.O) 


II  I IB  i-Ml'"  I 111 


IFCNTBIUNE  .0) 


IFCNTBO.NE.O) 


dJI  HrfCIWfi 


I l.'im  I 


IFCNSHAFT.  EQ.  1)  NRITE  CNN)  ITE,22>  CGC1 ),  GGC2)  , RSMAFT  .LSHAFT 


1 READ  CNREAOylZ) GGCl) ,GG* 2) , XHUBCR.YHUBCR,  ZHU8CR, RXHUB ,RVNU8»RZHU8 

rrcuuna  rn  <t 


1MRITE  C NNRITEt  ZZ>  GGCl)  ,6i  C Z)  .XHUBCRt  TNUBCRf  ZHUBCRt  RXHUB, RYMU0,RZHUB 
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READ  ( NREAO,  II ) GG<1>,G;<2>,NNAK»Y,NNAI«N,NMAKSIC,N.I»I<81 


IFUNWAKMB*NNAKSK*NNAICB.)  .HE.O) 


XF<(NNAKHB*NNAKSONNAKB.)  .HE. 01 


IF(NMAKPY.EQ.l)  REAO  (N!i  AO.lStGG  Cl ),  SWi  > .NSTAG,  NVO*T 


mm  ui 


IFiNNAKPY.EQ.l)  REAOtNREAO,12»GG<ll,GG«2»,CGTAG, CROAT 


lk  FOKMATt  2AA.IB.2F0. 3) 


RCAOINRCAO.IO)  GG(1I  tC;  (2)>(KP*MT  Ul  ,01.101 


IF  INPYLON. HE.  0.  OR.  N*TAI..NE.O>NS«NS»MPYCON*N»TArLf*AULTY 


IF  (NROTOR.NE.O)  NS*NS*<5lAOEMNSNANO«lAOf»NULT» 


(S>K$M 


rv<  n 


READ ( NAEA0.lt)  GG(lk,G»  (2)  .NX ««»  , NY  ( KSt  , KNSELC,  KNSSHP,  KNSTYP 


NUff  (MMRITCi  ZttGC<ll«ta<Z»tM(<KSl  ,M»t  »CS>  .KHSELE.KMSSHP.KMSTYP 

mxyiksi»mx««»»mtikxi 

ZSPfM«KU*l 

«MO>m..XKSl»l 

KMAKES(KS)  »• 

amaMmtsjjq 

THKYmn.yc,  nai  rn  «xn 

NSBOOT«NSBOOVM 

.n»xsnMamrntMannT?«Mr,iYi)*mn.T>nmo>i«MMiAn 

NX(JS)*MX(KS) 

MTUSlAMTXXSl  . .. 

NXY<JS>«MXY<KS> 

i space ijb>«9 

KNAKES(JS>*t 

XMOMMLt.JS)«.l 

KM0T0*S(J3)»B 

ua ronTiMin 

IFCKSYHMZ.NE.  SICO  TO  199 

usaamnsaonmi 

KSL*KS*MBOOYl*N8ODY2»NB5DYJ*NPYL0i« 

MX  KSLUMIIESI 

*YtltSL>»H»«KS» 

MIT  IX  SI  >«MXTfXSI 

ISFACE  IKSL>«9 

EMAEf  SUCH  llfl 

KHO*m.«SL)«-t 

xAOJmasixsx.ua 

c 

IF  lirSTBMY.MC.BXCO  TO  XM 

NSBODY»NSBOOY*t 

Jsi..«m-AiMitnnTiAM.‘ianT2»xionT3i»mn.T*MPTinMAMiiTAri 

NX(JSU>NX(KS) 

MT  UIUMIUttl 

NXYX  JSX.»«NXY(KS> 

T SPACE  I.ISlXxl'T 

KMAKES(JSU»B 

KMnam  i-isi  i.  i 

K MOTORS  I JSL)*B 

X39 COMlIMUi 

C 

C mmmmmmmmmmmmmmm  flltLACt  — ...  mmmmmrn  mm  mmmmmmmmmm  mmm — . 

C 

iriMinnTT.Frt. men  to  ?m 

KS>RM 

MSBOa»»MSBOIW*l 

MFAOXNPEAO.IS)  6C  (1 » ,C(  21  ,NX  (KSX , MT  ( KS1  , KBOEIE,  K90SHP  .xeOTYP 

MBITt  ( BMRITC.  ZAXOCIl  X.C.t  2)  .MX XttX  , MT  IBI  .XBnCI.E.XAQSMB.xanTTB 

NXT<KS)«NX<KS>*MY(KS> 

ISFACPUSXx? 

KNORHUKSIM 

XMOTOBSIXSIbB 

KMAKES(KS) «t 

X 

IFCKSTPMT.ME.  t >60  TO  2P 

MS«anT»MMnnr»i 

JS«KSt(NSOOTl*NSOOVZ»NSOOYJ)*mit.T»NPYlON»Nmil 


NXY<JS>«NXY(KS> 


KNAKES(JS)«t 


X ROTORS  <JS)»t 


IFCKSYHHZ.NE.  (ICO  TO  299 


KSL«KS  »NOOOY1  TN800Y2  *M*3  O Y 3 ♦NFYL  09 


NY(KSll»NY<KSI 


NSB00Y«NS  BOOTH 


NX(JSL)>NX<KS) 


NXV(JSL)>NXV(KS) 


KWAKES  ( JSL I *f 


X ROTO  RSI  JSL  )«t 


REAR CNRCAOt  191  GO  Cl)  . Oil  2)  ,NX  (XSI  , MY  t XS>  tXTNECEtXTNSMPf  XTNTYR 


n * n 


1 1 IjIjN  i 


NXY(KS>*NX<KS)*NY<XS> 


KNORNL  IKS)  •! 


KR0T0RS<KS>*9 


IFIKSVHNY.NC.aieO  TO  379 


JS»KSMNtOOYl*NBOOYt*NBDOYJ>»NULT*NPYLON*NVTAIL 


NY<JSI<NY<KS1 


ISFACE  <JS)  <11 


KNORNL  IJS1<-1 


379  CONUNUE 


NS800Y<NSB00V*t 


KSL»KS»NeoBn«NU0Y2»Natmi*«rrt.0N 

UVIVtl  lallVfliTI 


NY(XSl»»NYOCS) 


ISFACEOCSLW 


KN0RNL(KSU«>1 


i iiUiT  hi  i 


nsbooy.nsoooym 


NX(JSU«NX(KS> 


NXV(JSl)«NXYf  KS> 


1 f 1 LH . i 


REAO<NREAO,11)  GC(1),C(Z>  ,NX(XSI , NY  ( XS>  , XPYELE.XFY SMP , XPYTYP 


MXY(XS)»NX(«)»NY(XS» 


KNORNl  CXSI  •! 


X ROTORS (KS)*Q 


IF<XSY»NV.NE.  0)60  TO  499 


JS*XS* (N800YI ♦n0OOYZ*NBjOY3) •NULTTNPVl  0N4NVTAIL 


NY(JS)  >NV(KS> 


ISFACE  (JS)  *12 


XNORNL  C JS>  *-l 


111  III  MW  1 


499  CONTINUE 


N BODY  S*NBOOYl ♦NBOOYZ ♦N950Y J 


VERTICAL  mi 


NsaooY*NSBoor  *i 


xs»nbooys»nult*nrylon*i 


MRITE (NNRITEt ZAI GS 111 1 6i ( 21 1 NX  (XSI  t NY  < X5)  t XVTELE t XVTSHP,XVTTYP 


ISFACE IXS)*» 


>■< 


f 


r 


-C. 


c 

-saa 

c 

x — i 

c 


c 


699 


K ROTORS (KS)*0 


XFtKSYHNY.NE.  8)60  TO  99) 

MSB00T«MS8QOT*1 

JS«KS  ♦ ( NBODYl  *NBOOY2*N8j  0Y3 ) •HULT»N  PYL  ON*N  VTAU. 

Nxt.m«nxnrs> — 

NV(JS)«NV(KS) 

MXYUS>»MXTm»».  

XSFACE(JS)>16 

mumwtfut — — 

KNORNL  (JS»  «♦! 

YB,0 TORS  <■»>»« — 

CONTINUE 

—■■■■■—  ROTOR  SHIFT  — 

rr  itiCMAFT.rn.fi  1 cn  rn.fcaj 

NS80DY*NSB00Y*1 

KSxAitaOORlAMRaOYZAltaaORU  »nm  TYMIII  TYAltRYi  OmWOMU  1 YMIU  XT. 
KS*KS»1 

at.n.ttc.11,,1 1 ecu  > T e~.  / o i -mu  ire  1 ■ MR  I itci  -xcuKi  t ■ rCMCyn.ITU 
MRITE(NMRITE,ES>G6(ll,6iCE),NX(KSI,NY(ICSI  ,KSHCLE,KSNSNP,KSNTYP 

XXTtKS)  » NYtKS)YNYtXS» — — 

ISFACE(KS)  • 17 

KM0RN1<KS)»1 — 

KM«KES(KS>*0 

KROTORS  USlil — — — 

irnrsYMKY.MC.8na  to.  am 

NS800V«NS800T»1 

JSXSaa — 

NX(JSI>NX(KS) 

MY<J8»«MY«RT>  

NXY(JS)*NXY(KS) 

ISfACCUStMlR 

KMAKESIJS) *0 

WMW8U4 — 

KROTORS (JS)*1 


CONTINUE 


ROTOR  NUB — 


IF<NHU8.EQ.R)60  TO  799 

«»— ■ »./»«■  amamutati  mmtri »mi»  tt 

XS«XS*1 

»tMM.»«MII».l 

RE  AO  ( NRERO  .14  > 6611)  ,6611)  ,NXIXS>  tHYIKS)  ,XM8ELE,XN8SMR,XN8TYF 
M.tT«i»a»iw,»»»a:nircni.»iii«ti.ii»fim  .xmRri f.xnr cmb. kmbttR — 

NXY(KS)  »NX  (KS)  *NY  (KS> 

ISXiCE IRS) »)l  — 

XNORNL  (KSI»1 

muncui — 


KROTORS <KS)«1 


tfcksthht.mf.  n ca  rn  ra 
NSBOOY«NSSOOY*1 

js*K%*nuimritiii 

NX  (JS  ) *NX(  KS) 

MYlISlNMYdfSI 

NXY(JS>»NXY(KS> 

iSFicetJsim1 

KNAKESIJS)«1 

KMO«iai.KlNTl 

KROTORS  ( JS)*1 

TT.  rnMTTMIIf _ 

IFCKSVNNZ.NE.  S>60  TO  T9» 

M<«nnT.y«mnrvn 

KSL«XS*1 

HX(KSl)»MXm> 

NV(KSO»NV(KS> 

MXT(KSL>MMKT(XS1 

ISF ACE  (KSL  ) *20 

KMAKESUSI  >»t 

KNORNL(KSU*-l 

KROTORSUSlUt 

rn. 

NSaonY*hSBOOVM 

■Kl  »-IN«MHUN»MIB  Y 

JSL>KSIHMUS*NUIT 

MX  I JU»H«I»I 

NT(JSU*NT(KS> 

MYYt  iSU«HXY(KS> 

XSFACE  (JSL)«22 

UlUtUJSLl  »1 

KNORNL  (JSL)*t 

KR0T08S(.)SL)»1 

T99  CONTINUE 


c ■ ROTOR  9LA0c  SHANK — — 

£ 

IF(NSHANK.EQ.  0>60  TO  699 

luanmtt  .-1/irmtuanium  mi  run  tyumRyi  Oman  XT  ATI  I f M1U.Z.X-* 

INSHAFT*HULTY*NHU9»HULT*UUTY*1 

1‘  » 66  (A  >.6;  ( 2>  ,NX  (XS>  . MY  > KS>  iKSKELE.KTKSNB,  KSKT-YR 

WRITE  (NMRITE,2S>  66(1  >.G(  2)  ,NX  (KS»  , NY  ( KS)  , KSKELE,  KSK  SMP, KSKTYP 


no  All  K-ltKBlAOE 

MSHPOYMUSRonxMi 

KSS»KS*(K-1J»(NSHANK*N8.  AOE*NULTI 

MY  (KSS)mMX  (KS) 

NY (K  SS)  *NV (KS ) 

M«VIKSS>«M»(<S)»MV(KS1 

ISFACE(KSS)«2J*(K-1»»J 

KMORM1  (KSS)«1 

KMAKES(KSS>*1 

aanuamu 

SIS  CONTINUE 

£ 

S99  CONTINUE 
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NSBOOY»NSBOOYM 

*nftwar> 

NX (KSS)  *NX  (KSI 


n Ml  *.'4.-4 *r 


NXV(KSS)«NX(*S>>NV<KS> 


KNORML  (K$SI*1 


X MOTORS (KSS)«1 


IFCKSYMHZ.NE.  0>GO  TO  9U 

MsaooT»MsaooTMi 

KSSI ■ KSS+1 


NY(KSSl)*NYIXSS) 


i til 


ISFACE<KSSL>«Z0»(K-1>»3 


KNORMUKSSO»-1 


9 tO  CONTINUE 

n 


999  CONTINUE 


MRITE  <6,00(21 


KIT  I Bini 


1100  NELEN'NELENYNXYIIS) 


0001  FORMAT  <//5X,*NELEN  « •,!>/» 


0002  FORMAT (/10Xf AHOATA) 


NSSTOT«B 


NT(IS)«NSBTOT 


NRITE  (6,00(01  IS, NX  (IS)  ,Y<  ( ISI  ,NXY(I  S)  , ISFACE(IS)  , NT  ( IS)  , 


c 


c 


c 

£. 


J. 

-fl. 

c 


tUMDUTTi*  rinngT  lufi  rm  u«m.or  yntr r , inr  rmi>nn 

r.n»tM<m/?77i/MX(ifc),Mt(»)  ,nxy<tni  ,kctmm>,ics¥mnz.mstnmt,mstmmz 

CONMON/ZZZ?/NSFX,NSBOOY'  iS,NT(3%J  ,1  SFACE  { JN>  , KNOX  ML  ( 3X>  , X BAXES  (3  A) 

cnMM0NiZ»3AiABYU0M.M9nan  ^nnoom,  mb  not  t.nxtxil.  mimaki, MMua.NSHAN  it. 

1 N81A0E 

r.nmtouj777k/\iMkCM.n**r.L  iin.iafTi 

C OHNONX ZZZ9/XPYCTR,  YPYCT  X > ZPYC TR,  XX  PYL  ,XYP YL , RZPYL 

COMMON/ 777BXIM0Sf  .XBOX  ,HM»XIAII 

C0MH0N/ZZZ7/YN0SE,Y801,O0Z,YTAI|. 

nniuinM/77ZH/ZMn^ff.T«nt(fiM.7Txii 

COHHON/ZZZ9/tYBDl,XZB01.  TYB02.RZBD2 

COMMON/ 277 IB/ BSHArt,.  I Sit  :T.^«!HllKfl  SHAM 

C0HN0N/ZZZ1I/  XHU8CR,  YHU1 CR  , ZHU8CR.  XXHUB,  XYMUB,RZXU3 

common/ 7zzt7/BBotna,  board.  nmt,>L»» 

C0HW0N/72Z13/T  NET79,  TMCTlC.THETlS.CONINS.OZlNUTN 

rj\MMnuj777tujwai  inr,r.y  ;».t.yTc.;«ai  ■-r«afcrwMiaiHki 

C0HN0N/ZZZ15/VSPANtXLEZY  ,XTEZ¥, YAXLEY,  TANTEV.T AUY ,ZVT AIL 

CQhmM/ZZZIft/MMAXPT.NMAna.NMAKSE.MMA'cai 

C0MN0N/ZZZ17/M*KLPY>W*K.4S>NAKLSX.M«KlIL 

r.OHMOB/ZZXlltf  MANGPY.MAMT.AB.MABSSX.MABCa 

CONHON/ZZZ1V KPXINT  (lit,  NXE  AO,  NMXI TE,  KPEAO 

rnmnM/7?T,«/DT 

CON*ON/ZZZ21/KPY1,<PY2 

oonMflM/777»/XMxri  xfxn<r/P.xM<Tr» 

CONNON/ZZZZ3/XBOELE.X8nS/P,XBOTY» 

COMMON/77I2NZKTNftf  , xtm;mb.xtmtyb 

C 0NN0N/Z7225/XPVEL£t  XPYSBp .KPYTY* 

uumau^tzuuuuxtLS^uat  a»-^uxcu 

COMMON/ ZZZZ7/XSMELE,XSK?  IP  .XSHTYP 

COMMON/ 777  7/1/ XHflf  I F , KHIKHB  »XH8  Tti 

;ohnon/Z7229/kskele,xsk3/p,xsxtyp 

.aixina/ 777.iaxxfli.fi  f , xbls x b » xbl  ty? 

COHMOM/7ZZ31/NSTAG,NYOXr, MSPIXAL«$PXXAL 

rnMMnM/77TH/MMAXC 

CONHON/ZZ73WYCUT 

CQHHQH/Z7ZS0/  XCTZ.TCTX."-.Ta.BX,XY.i7 

C0NH0N/ZZZ51/HSPAN, XLEZ1  .XTEZH.TAXlEH,  TAhTEH,  TAUM.ZP 

connon/777iv  iTMisr.aBircu.Tin’ST 

DIMENSION  ACIOXOIZO)  ,MN>P(Z>  ,MNPN(2> ,MNMP<2>,NNNN(Z> 

otmemstom  XKU.mmnri.mxr-TmxfB.iAi.Gtizai.x’aAcz.m-i 

DIMENSION  0P»P<Z>,0FPNI’)  , OFNP  (21  ,DFNN  «Z) 

fiTNfMSTtin  HCMoanm i.haxisi n i.mcbqboi  ti  i.xaxishh 

OXMENSION  HCSI  (lll.HCTAflll 

XEAI  LSMAfT.LSMAMX. 

XtlteM 

KXEXOM 


IP(*XEAO.EO.O»SO  TO  9 

»fAO<imu,ui — inxix,,  xoaE».K«i,»,iMaot«i,MMooi» 
MRITE  INWUTE.25)  t(XK  <K,r  Y00C»  ,X«1 , J) , IN30E»1.  NNODEI 

MAAimiuja 

SO  TO  MO 

CONTINUE 


I 


REFLEM1 


19  FORMAT (10Ft<3) 


PI«S.  1A199269 


1 •MOFCT**t  X91 


OTV«l./Ht(KS) 


NYP»NY  <«>♦! 


IFCLPYELE.NE.  I)  SO  TO  l:  J 


WRITE  ( N WRITE,  29  ) (HCSIITO  ,IX*1,NXP) 


WRITE (NNRITEt 291  (HETA(IT)  »IY»1*NYP> 

PAUR*  Ul I* 


00  lit  IX*1,<(XP 


xx«<ix-i>  *ox« 


IF<KPYELE<EQ<  0)  XX.HCSI  (I  O 


K I * I I 1JT 


SO  TOC  131 » 132 < 133) t LPYiLE 


60  TO  115 


csx«xx*xx 


153  CONTINUE 


THETA  »PI»ETA 


f- » rr tiiii  i 


RY*CSI YRYPYL 


Tt«RT»SIN  (THETA  > 


J/n  >J  If  (I 


IN(SE«N0FCT(IXV,KS1 


n 4*^rTiT»  4m 


XK(Z,XNOO£l-Y«»YPYCTR 


MRXTE(t,ltttXX,IV,IXY,(r;i,ETA,IH3DE 


NRITE  (6. lVZtXBf Vt.ZB 


I lal  W I'tll  4 HP*4TiliI  Jl 


1M  F0RHAT<l«X,»IX»*,IS,ZXt»tY*»rIS,2X,»IXY»»,I9,  ZX,»CSI«»iElZ.6»2X, 


Ml  F0RNAT(1IX,*THETA«»,£1Z.  9,ZX>-RX>-,ElP.6.ZX,>RY>«,clZ.il 


1*3  F0RHAT(19X,»XK«»,3(E12.<>,  2X)> 

rnuft  UllC 


199  CONTINUE 


IF(  XBOCVl.EQ.lt  SIGNZ-’l. 


ii jii  <iit*n j»^jv44’v»rnF 


NXP-NX  (XS1  *1 


OXX»l ./NX (XS> 


00  990  IX« 1 *NX° 


IF(KNSEIE.EQ.1ICSI*XX 


X«XN0SEMXS01'-XN0SEI  »Cr 


IF (KNSSHP.EQ.lt  RR7>RZSOl>SQRT(CSII 


• LV(ti44liJ 


IF(KNSSHP.£Q.  2tRRZ*RZ301*  (CKI1°°.3333Y33 


I F (KNSSHP  . EQ.  3 IR9Z  «R  ZB01  • ( CS II  ZS 


THET-(IV-1I»TINCR 


I 4 4 #TTfr 


i ■ i 4 4 


ETA-COS  (THCTIt 


ZZ-SX  G PZ’RRZ*  SI  N (THE  TA I 


OONHZ-OELZ-OEl  Z»CSI 


rrf4Ci'ti  pJ 


ZZ-ZZ»ZB01 


IF  (IBOOT1.EQ.  It  IS-KS 


fec^aii  r. in** 


■PfJKrl  441  4jW  i 


MRITE(G,M(lIXtXV,IXY,CTCtETA,NOFCT(IXV,ISI 


NPITE<6,l6Z>X,YYtZZ 


t Mftrt  a 


XKIEtNOFCTIIXr,ISII*YY 


699  CONTINUE 


*99  CONTINUE 


SOOT- 


IF  (NB00Y2*EQ*9  HSO  TO  69 


00  696  I800Y?»1  ,»JlT 


J‘WH.Urr>/TlF 


I F (IBOOYZ.  EO. .?  >SI6YZ*-l. 


NYP*NY  ( KSI  ♦! 


OXX»l  ./NXIXSI 


XXXa(IX-l)*WX 


IF  (K80ELE.E0.  Z>CSI*XXX**«X 


Az««zm»(*z9ot'<tZ90i>*:si 


ZZZ*OELZ*CSI 


00  6St  IY»1,XYP 


YY-«1»C0S(TMET*I 


ZZ«Z7*ZZZ*ZB 01 


IFdBOOYZ.EQ.  1 >IS«« 


XK<1.  NOFCT  CIXTfISI ) « XX 


XX<3.N0FCT(IXV,IS>  )»ZZ 


696  CONTINUE 


>t*il  nose 


IFINSOOYJ.EO.  9)  GO  TO  990 

IF  A • 


00  999  I BOOTS*  It  WILT 


I Ft  IBOOY3.  EQ.  ?> SZCNZ »-l 


NYP»NYOCS>»l 


TINCR*9.S*PI/NY(KS) 


XX*(IX-1)»0XX 
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83® 

HI 


xx>t.-xx 

LE.IXTXEIE.EQ.1  WTSTiXT 


IFtKTN£LE.EQ.ZlCSI«XX»XX 
IKXTX3HB.Fa.ll  RRYaRYBIgf  3QRT<C«> 


I F IX  T NSHP  . EQ.  11 MZ  «R  Z8D?  • SQR  T ( CSI » 

HHTUIMI  tn 


IF  (KTNSHP.EQ.  Z)  RRZ'RZBOt  *CSI 

m«i  .west 


X»XBDZ*tXTAIL*XBOZ> *csx 

aa  bib  »ii.x» 


THET-«rr-l)*TIHC« 

Ufn.I.IMT.TUtf' 


ETAaCOS  (THETA) 


ZZ«SIGNZ*RRZ»SIN)THETA> 


on  T»ZT«IL«7»nz 


OONNZ*CELZ»CSI 
JJ»7ZtnnHMZ»7BOZ 


XXY«IXMIY-1)  *NXP 
iInifliwi.ni.nK.rt 


ifubooyj.eq.  z*  is>ks»kfa:e 

ltn.anfrTlT»,Ki..» 


XX(Z,NOFCT(Ixr,XS>  >-YY 
EE(T.N0EfeIA»T.I3>>«ZZ 


894 
-C 


CONTINUE 

continue 


CONTI NUE 


VERTXCA..  TAX.  ........... 


IF(N»TAIL.E0.®>CO  TO  39i 
r3KMannYi*yannyj>iianf»vnYmu  LAitatmmu 


00  398  IVTl*1,1 


IF  (KVTTVP.  £Q»  1)  GO  TO  31*8 

^FAn  ( nbfao,  13) — cog  i.kkzi.hhbp.iuibh.  mmib^mnni. 


■ltio 


NRXTE(NNRXTEt  ZS)G6(1),CS(Z)  ,mnpp,hnpn,  hnmp.hnhh 
■ OOMTINIIE 


SISNZ-M. 

i eh  i mttMMiZJa 


tailno-tailln/z. 

■RL£»«. 


RTEM. 

■ nxxai.^inuxxi 


OYY.l./NY<XS> 
MXP«MX(X3»  «1 


NTP>NY(KSM1 
-4.  XIFLFaXVTELf 


IFILVTELE.NE. B ICO  TO  3118 
i am  r.i 


REA0(NREA0,1I)  GG(t)fG3(ZI»  (HCSI  (XX)  ,XX>1,NXP) 
MBTTF  IHNBXTF.Z3IB0U  > ,fir.(  ? ) . (HC3K  I»)  . Hal  ,MXP) 


R£A0(NREA0,1S)  EG(  1 ) , (T.  ( Z) , (META  (I Y)  , IY»t,NYP) 
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WRITE  (NNRITE,  25)  GG<  1)  , GG  ( 2)  , (HETA  (IV)  , IT *1 ,NYP) 

f*  nu  »»  LAI  1C 


so  sit  ix«i,*xp 


XXXIX-l)»OXX 


IFIKVTEIE.EQ.  I)XX*HCSI<1  x» 


I fi  4 * trim 


■ ■ w * 4#  ■ 


60  70(3131,3132,3113, 313%)  ,LVTELE 


60  TO  313 


CSI«XX*XX 


60  TO  3 13 


60  TO  SIS 


CSI*XX«XX 


313  CONTINUE 


Y»T*n.wo»rr* 


XTE«XTE2*AT»NTE»*(T-RTEl 


VAXIS(IY|»XLE*VCHORO(IY)  *1.5 


IFOCVTTYP.NE.  2160  TO  31  3 


T1»(NNPP(2)-WNNP(2)  )»CS.  *WNNP(2) 


««  I’T fill 


Y2«(NNFH(2)-4NNN(2)l«CSl»tfNNH(2t 


Y3»<WNF»(2>-MNPN(2> I *ET1  ♦ WNPN(2) 


T4XNNNP<2l-*Nl1N(2>>*ETl»«INMN(?> 


ill  jn>  < n fjut  »o nm  »a<  sm in  ria 


I F (0ET.EQ,  (•>  GO  TO  3147 


ii  (it  m j it x j 


YX(Y1-Y2»»CXV»Y3-X3*Y4I  -CT3-T4)»(X2»T1-X1»T2))/0ET 


3167  CONTINUE 


3151  CONTINUE 


FOR  CIRCULAR  8IC0NYEX 


PL«0.3*CXTE-«LE> 


ATAU»t.»M 


iF<H.ca.f..M.tY.ea.Hvp.}R.ix.ea.i.ait.ix.ca.NXP)so  to  jh 


I'M  J JBT  X I 


IF«£TA.GE.ETACRlt*?*SQ«T(l.-«<ETA-ETACR»»YSPAN/ATAU>  ••?! 


315?  CONTINUE 


?«SIGNZ»TAU8AR»SQRT(CSI>  •ll.-CSI»»SQRT(l.-ETA«t> 


3153  CONTINUE 


315  CONTINUE 


IXY«IX*«V-1>  »NXP 


XK(?,NOFCT(IXY,KS)  l*Z 


318  CONTINUE 


■ h»i7i  ■I7T»i  ■ 


399  CONTINUE 


SNAFT 


IF1NSMAFT.EQ.01GO  TO  ?9» 


0YY«1./NYIXS> 


<1*1  1 ••  i 


NYP*NY<KS>»1 


IFILSHELE.NE.  0)  GO  TO  XT  Li 


LSHCLC-1 


NR1TE  (NNRITEt  ?5I  (HCSK  IYI  *IX*1.NXP) 


T wJ>  Tf’  111 


NRITE  (NNRITC.  ?9)  (META ( \<  ) . IV*l,NfP> 


00  tit  IX-l.RXP 


XX*UX-l)»OX« 


irTii  ail 


IFlttr -Lt.EQ.tl  XXmHCSKIXI 


GO  TO  (Xll.tSt.tni.LSHEuE 


CONTINUE 


CHHM3B 

ETA-TY 

60  TO  MI 

in  CONTINUE 



ETMTY 

i — cnMTiNur 

TNETA-FI*CTA 

MUNIR 

XB»R>COS(THETAI 

TB«BE9IM(TMETAI 

ZS-CSIM.  SHAFT 

— mnmintiiwa 

INOOE»NOFCT(IXV,KSI 

imi.IMOOSIIBAXBXCTP 

XKI1, INOOEI-T  8*YPYCTR 

Xl<l,IMOOri»tBTCTa»CSHA.T-IBABIBrL 

UNITE  (6.61IKS  ,IX,IY,  IXY.  I NODE 

U BITE  <4.411  »,»!>» 

MRITE<6,63MXK(LiIM00C>'.*1.3> 

4L-E0BB4H1BX.BKS  — ,IS.B (■»,»,» — — »»*»,», 

l • INOOC»*,IS) 

41  4,  » »H,I1I.4I 

43  FORMAT (19X.*  XX-  *.3  (E12 . 4 , 2X1 1 

21 A CONTINUE 

290  CONTINUE 

£ 

IF(NHUB.EQ.t)  60  TO  499 

W<NaOOTl*HAOOT|4MEOOWI  TNULItNU.tr*  (MBTLOlUMIfTA  IU  BMULTT* 

1 nsmaft»nulty 

KSISAt 

00  *90  INULT-l.NULT 


IFdNULT.EQ.il  SI6NZ-  1 
rfdNim.iq.il  ti6Mi— i. 

c 

niiNi  ./mxiksi 

OYY-l./NYIXSI 

NXONNX (191*1 

NYF-NV  (1ST  *1 

L NOEL  t IMPEL  t 

IF (LHBELE.NE.  9IG0  TO  ABB 

LMBELEnI 

REA0(NREAQ,15)  (HCSKPO  .IX-l.NXFI 

MBITS  ( MMBITSrlBI  INC9I  ( IT  I .IT-1.NXAI 

BE  AO  ( NNEAO  * 19 1 (HETAdYl  , IY-1.NYF) 

MBITS  (MNPITS.  WmeaWH  fOMMX* 

ABB  CONTINUE 

00  411  IX«1,XXB 

00  41B  IY-l.NYF 

XXN(IX»UtOXX 

YY-dY-ll  »OYT 

IfdMBELS.  Efl.  El  TXiMCSIUXI 

IFIKHBELE.EO.BI  YY-HETA([TI 
r.n  rn  i m rA.i>,  mi , ■ -»  r 

431  CONTINUE 
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WHITE  (6,51.)  I*r,  IS 


XKI2,N0FCT(IXY,JS>  >«-XK(?,NOFCT«IXT,ISM 


912  COMTIMUE 


isi«nss«mn.ty*nsh*ft*nu.  TT*1 


00  916  IS*IS1 , IS? 


NXP»NX«IS>  *1 


00  919  IX>1.NXP 


ixy*ix«uy-i>  »nxp 


XK  (1 1 NOFCT  (IXY, JS) > * XX(  l,NOFCT(IXY,IS>) 


i j m m h * a m 


t*  iMMll 


XK(I,NOFCT(IXY,  JS>>*  XKI  3.  NOFCT  (I  XT  ,XSM 

r'ntiTTuiir 


916  CONTI NUC 


KSKN*OOY1»NBOOY2»NSOOY3I  •WLT*NULTY*  ( NPYLON*NVT»  IL)  •HULTY* 


00  T9«  I#»1,X9L»0£ 


IF(I9.C0.1ISI6NY>1. 


IF(LSKELE.NE.S»CO  TO  TOO 

• r m 4 


*E«0(N«FA0,15>  (HCStdXI  ;IX»1, NXP) 


*f»fl(NQF*n,l?*  (KTMBt  sIYsl.NYFJ 


T9*  CONTINUE 


00  71*  IY»1,NYP 


i n a 4 n ii  x i 


YY*IIY-l)*OYY 


L 4, 


IF  (KSKEIE«E0>  It  YY*HET9(IXI 


CSI-XX 


*#««•  COS(TtCTX) 


4-n  k.Jk-»4,T4  <M||V 


Y0«CSI*tSM*NK 


ixr>ix»trr-n*HXP 


iHooe-NOfCT  (ixr.is) 

xJKt»ma«>«xMmcaxx< 

XK«,iNOOt»«si6Nr*frHOB:<*Rrw;8*LS»  «nk-v«> 


1(0*  continue 

no  1B1A  IXal.MXg 

00  till  IY-l.NTP 

mumumma 

VT«(IV-ll*OW 

immU.M.Il  iWNCtiVIXI 

IFIKBLCLE.EQ.il  VY*HET Af ITI 

60  in  Ma.Ti,inz.iavT.r tm  ,i.«lele 

list  CONTINUE 

want 

ETA*VV 

60  IQ-  .101 T 

1112  CONTINUE 

ciimm 

ET»«l.-(l.-YT»**C 

60  

1131  CONTINUE 

MOW 

ETA«».-C1.-WI*»* 

60  TO  tan 

1136  CONTINUE 

cnm»»»» 

CTA*VV 

1111  continue 

aCH0B0*XBTE-X5L6 



irnxin 

mwuim 

TAN|TE*T  ANTES 

TBAIUI.WWOTOO 

▼*<M0T0*-VCJT)  *ETA*VCUT 

»i  tn*  i f 7.T1M<I  f>l».»eilT 

XTE*XTEZ*TAN8TE*(  V-VCUT) 

ACMOin-KTC.lLC 

X«XLE*fCMOKO*CSI 


c 

60  ta..llM1,lH?(HS31«XAl  SHI 

1951  CONTINUE 

£ 

C FOB  ClftCULA*  BICONVEX 

■C 

XCM.S*(XLE»XTEt 

8l«A.  HIXTE.XLE) 

TAU8L-.S5I 

TAlf  TAUIL 

H«TAU*PL 

■ ATAIlMl.11  ...  ■ - — ■ 

ETACIt-t.-ATAJ/lPAn 

XXC»l«Xf.  , , . 

ZM. 

If  <0.60.1.. flK.Zt.iO.KTO  TA.TX.EQ.I.OA.  I«.fQ.NXH60  TO  Ml 

Z»3I6NZ*<S0BTC<<Pl*PL»»P-n/<*.*Nll**Z-XXC***>-<Pl*PL-H*H>/<C.*HM 

HIHA.6I.  ATACBI  I»II  TOOT  A 1 ■ -U«IA.1TA^A>  11BAN/ATAU>  I III 

60  TO  115 

murmur 

TAUBAMTAU*.  T 5 *501(1  ( 3.  I*  ( XTEZ-XLEII 


141 


Z*SIGNZ»TAUBAR»S0RT<CSI|p«1.-CSI>»SQRT  «l.-ETA*»l) 


1(93  CONTINUE 


11S  CONTINUE 


XXY*XXMIY-1)pnXP 

TCaVt 


IMINULT.EQ.7)IS»KS*1 


XX II » INODE) *X 


XK(3, INOOC>«Z»ZHUaCR 


AZXN-AZINUTH»PI/l(t. 


PITCMR«RPITCH»PI/18(. 


• m!T1  Pit* 


XXaXK I ltINOOE ) 


ZZ*XK CStINOOE) 


TIE  O tJF/I  in 


CZ«(VY-X1>/(XZ-X1) 


PITCH»C1»PITCH«*CZ»PIT0IT 


3(Z  CONTINUE 


PIT7S»(THET7i  *TMET  1C»C05  t AZIN)  *TNET  1S»SIN  I AZIN) ) »PI  71(8. 


XZ*.75»RR0T0Z 


rrm  «<■ 


YY«XKI7,IN00E> 


r I fl  4 llfli 


C1*(YV-XZ)  / CX1-XZ) 


3(3  CONTINUE 


C 3LA0E  TNI ST 


COSPXT»COS  (PITCH) 


XKJ-XXI J.INOOE) 


XKIJ,INOOE)»-XX1»SINPIT*XXJ»COSPIT 


C CONING  ANGLE 


CONE«CONING*P1/180. 


XK1«XK(1,IN03C> 


ini  » manna  un»u.nK*siiu  mum 

xx(?,  xncoe)«  xki’sincan;*)  »xk?»cos<  «ngp» 

«« 

WRITE  <6,5*>  19,  XNULr,IXT,IS«X,Y,r 

1 

mu 

L •XYZ»*,3IE12. 6,2X1) 
rnarfMiir 

1M? 

j $aj» 

CONTINUE 

CQffTXMlIf 

1999 

r 

CONTINUE 

nrx re («•««> 

c 

fltriiay 

ENO 

143 


COMMON/ZZZ1/1XC38>,NT(»>  , NXT  <J») , K STINT , KSYHNZtNSYMHT  .NSTMMZ 


C0MHON/ZZZ3/MPYL0M,N8O»l  ,N0OOrZ,MaOOY3,NFTAIL,NSM»rT,NHUB,NSHANK, 

UOI  3 FIT 


COHHON/ZZZA/UNACH.ONCCA.  ALFA,ASETA 


I.  i»T  *±***4^1  1 If  HI  1 i A il  »i  Til 


C0MM0N/ZZZ6/XN0S£,X801,<  JOZ.XTAIw 


l.  til ' T * 0 m m Mi  I riV-  M / 1 >TI 


COMMON/ ZZZ8/ZN0*E»ZS01«ZB0Z,ZT  Alt 


COMMON /ZZZ10/BSHAFTVLSM  - T,PSMAN<tLSW4NX 


rii/iirii'  II 


COMMON/ ZZZ12/  RIOTORt  BCH  *0,  TAUBL,  At  FAT 


C0NN0N/ZZZlA/KaLA0E,TANLiB,TANTE3,x8LE  ,X aTE.XROTORS ( J*» 


C0HN0N/ZZZ16/* 


r.NNMiS.X 


C OMHON/ZZZlb/  MAMGPYtMAM» 18  yHAMGSXt NAN*  BL 


i j in  in  iv 


COMMON/ZZZ20/PI 


C0NH0N/ZZZ22/KNSELE,XNST1P,ICNSTYP 


LjJiI,  r##8lAfi  i:li  il 


C0MM0N/ZZZ2b/KTN£LEt  KTN5  IP  .KTNTYP 


111  il 


0NM0N/ZZZ26/RVmE,KmiP»KVTTT» 


COMMON/  ZZZ28/  K HBELE,  KMB5  -IP , KMBTY? 


I 4.-1  i-lljf  4-i  cCI'j 


COMMOM/ZZZ30/  KB  LEI  E,  K8L5 1 P , K8L  TT» 


Til  F/##Airl'lTI<’M  11 


REAL  L SHAFT tLSHANK 


00  IBB  IW.1S 


IIJI  I lli.II 


H M iHH  II 


IF  (NT  (IS)  «6T»  NT  MAX)  CAl_  OEBU6(3tO) 


IF  (NCLEM.Cr. NTMAX ) CALL  OEBUGIMO) 


80  FORNATUOX,p<—  END  OF  CHECK •) 


THE  CO-OROXNATES  OF  MOOES  IN  THE  GLOBAL  NUMBERING  SYSTEM 


NUMBERING  SYSTEM  FOR  IWIYIOUAL  SURFACES • 


OXMENSION  NOFCT<NXYNP, JWI 


COMMON/ ZZZZ/NSFXt  NS800Y,  TS  t NTC  3AI  f I SF  ACE  ( SA)*XNORML  ( 3 A)  f KMAXES  ( 3A  > 


III  wrn  1*  A Fi  11/ 


COMMON/ZZZ5/XPYCTR.  YPYCT  R , ZPYCTR,  RX  PYL  ,«YPYL  ,RZPYL 


COMM  OM/ZZZT/T  NOSE,  T901,riOZ,YTAXL 


COMMON/ ZZZ9/RY901,RZB01,  RYBOZ.RZBDZ 


C0MN0N/ZZZll/XHUBC.R,YMia;R,ZHUBCRtAXHUB,RYMU8,RZ-lU3 


COMMON/ZZZ1S/  TMET75  » TMETlCt  TMETlStCOMI  NG  » AZIMUTH 


rT»rr***A 


C0NH0N/ZZZI5/VSPANfXLEZVfXTEZV»TANLEV,  TANTEY,  TAU»ZYT  AIL 


COMMON/ ZZZ17/  MAXLPYf  MAX  1 B , NAKLSK.  M AKL  BL 


COMNON/ZZZ19/KPRINT(10».  iREAO.NNRITE.KREAO 
C0MM0N/ZZZ21/  KPYL1 . KPYL.* 


COMMON/  ZZZZS/KB  DELE,  XBtK  NP ,XBOTY» 
CONNON/ZZZZ5/RPYELE,KPTS/P,XPYTYP 


C O PNOM/ Z77ZT/  KSHELE , XSM?  -I  P,  XSM  TY» 


COMMON/ZZZ  29/RSKELE,  XSX>  RP, XSKTY® 


COMMON/ ZZZJl/MSTAO,MYORTtMSPIRAL,  SPIRAL 


CZA  MRITEC*»ZJA1NX<II,NYCI>,YT(I»,*H0RNL«H,ISFACEII> 


NTBODY»NBOOYl»NBOOY2»N85  0Y3 


cc 


Mpnft.il 

NOOEU-I 

1 KMT  toot  60  TO  Xi 

NOO€X«1 

oo  ?■  i»i.MT»amr 

28  NOOEX<NOOEX+YX(I> 

II  CONTINUE 

cc 

iKMunnr.MC.  »> 

lNOOEU*NOOEU*NOOEX»(NYlNT800YI*tl  -(NTBOOY-ll»NYINT90nYI 

X£ 

CC  NOO£X*NUH0ER  0s  NOOES  AL04J  AIRCRAFT  CENTERLINE 

xc 

CC  N0DEU.NUN8ER  OF  NODES  04  UPPER  RNS  FUSELAGE  BODY 

xc 

cc 

XC 

IS*0 

MOJif  »a 

NLAST.8 

X 

IFtNBOOYl.EO.  8>G0  TO  99 

IS»IS*1 

NXP.NXCISIU 

MYP»MYirS>»1 

00  >(  IY»1,NTP 

no  n xx«i,mxb 

IXY>IX»(IY-1)*NXP 


IF(IX.£Q.t.ANO.IY.GT.tlSl  TO  S5 

Mont«»Moo£a«i 

89  CONTINUE 

inooe.mooci 

IF (IX.EQ.1)IY00E*1 

MflfCT  l ITT,  I SI  <1X001 

89  CONTINUE 

_aa unNTTMiir 

c 

MLAST.MQOEI 

IF  (NBOOY2.EO.  II  GO  TO  19) 

»«nti 

NXP<NX(ISI  <1 

MY  Pint  (IS1  *1 

00  111  IYM1.NYP 

00  III  IT.l.MXB 

IXY«IX*IIY-1I»NXP 


IF (IX • EQa 1)G0  TO  1S5 

MonriauonriAi 

119  CONTINUE 

IMOM.NOOCQ 

IF(tX.EQ.l>INOOE«  (MX  <11*11  ♦<IY-l>»NX(ll 

XflACT  (I XT. Ill  . IMOOE 

118  CONTINUE 

X99 r.nMTIMIJf 

c 


146 


NLAST.N00E8 


klji  ! Tilik*«i>WlTJi»<.»I" 


NYP»NY(IS>*1 


00  288  IX-l.XXP 


IFdX.E0.NXP.8N0.ir.GT.llG0  TO  295 


285  CONTINUE 


IF<IX.E0.1)H00E» 


2 NXC2>*<IY-1>  »NX(2) 


1 MX(  1)  »<NY<1)  *1) 


NYP*NY  CIS)  ♦! 


00  38 1 IX>1,NXP 


IF1IX  .EQ.NXP)  GO  TO  385 


185  CONTINUE 


IFdX.EQ.lircOOEa 


IFdX.EO.MXP)  INODE*  NX (1 1 • <NY<  1>»U  ♦1»NYP-CIY-U 

unrc'T  t 


388  CONTINUE 


NCAST  «NOOCO 


IF(KS YNHI.NC.  9)  GO  TO  894 


IS«IS*1 


NYPaNT(IS)  ** 


00  89  8 IXat.lXP 


689 


IF(IX.EQ.l)  SO  TO  889 
1flIT.EQ.HYPl  cn  Tfl  h«« 


1M 


NODE 8-NOOE 8*1 
:OkHNUE 


INOO£»NOOE8 

tmuaumantix 


If  (IX.NE.l.EYD.IY.EQ.NYMINOOE* 

i mnniiriii.T. 


m 


MOfCT  (IXY.IS)  »IMOOE 
CONTINUE 


EM  CONTINUE 


NL*STmNOO€« 


IF(KSTHMZ.NE.O)  SO  TO  H 

ixuiaonTT.fa.8tfo  ro  >aa 


tnim 


NXP*NX(IS>*1 

MTP.MT  lf<l 


00  588  IY»1,TYP 

no .si a ixhi.nxb 


lxY»ix*dY-i»  *MXP 


IF (IX . EQ.l) G)  TO  989 
tFtlY.EQ-NYPmn  Tn  ««« 


188 


NOOE8-NOOE8M 
CONTINUE 


INCOE  aNOOEO 

icux.Ea.t.axn  iT.Nr.NnuNontp 


l NLAST-NX(1>*(NY(1>-T>  -(YX(  l)-l)-t»  (lY-t>  »NX(l> 
IFUX.Ea.l  .Mn.tT.fg.NT’IINfllga 


1 NX ( 1 ) • (NT ( 1)  »1)*1 

iF(ir.ia.Mr8.«Nn..x.NE.iUN00E»- 


l NX  ( 1 ) • (NY (11 
T NXIT1TNY(T>- 1*IX 


98fl 


999 
X 


NOFCT  (IXY.IS)  *INOO£ 
CONIXNIJf 


CONTINUE 


NLXSTeNOOEO 
IEtXSYHNT.NC.9IS0.IQ  40a- 


IF(N800V5.EQ.8)GO  TO  89) 

. »»«»( 


NXP>NXtXS>»l 

NTB«HT(IS)*t 


DO  888  IY»1,NYP 
—S  »ij,M8 


IXThIX*(IT-1)*NXP 


IF(IX.EQ.l)  SO  TO  889 
irdT.ia.NMico  to 


8X5- 


IFdX.EO.NXPISO  TO  889 

NonsapNonraxi 


continue 
INflOflNOnFI 


IF(IX.  EQ.I.EnO.IY.NE.NT:  • INOOC* 


148 


1 NL8ST-NXtZ)*!NY(Z)-Z)-<  9X<ZI-1»-1»  «IY-ll  »NXCZ> 


1 NXIll»INYIl»*t>*l* 


5 <NX(l)-l)*NY<3>»l*<IX-tl 


1 NXtl)*INY(l)*lt*l» 


3 NX  t 3 ) 


1 NX(ll*<NY(l)»ll»l*NX<Zt  • (NY(Z>U» 


MB  CONTINUE 


NOOET»Nt*ST 


IS«IS*1 


NYP-NYUS)  ♦ ! 


00  7BB  IX*t,NXP 


INOOE»NOOE8*Nl*ST 


78 B CONTINUE 


IF(KSVHNY.NE.O)  GO  TO  8)9 


IS*IS»i 


NVP*NY (IS) ♦! 


00  880  IX»1,NXP 


IF(IY.EQ.l)  GO  TO  885 


865  CONTINUE 


IF(IX.EQ.1>H0DC»1 


NOFCT  (IXV.ISKINOOE 


NL*ST»NOOEB 


IFCNBOOYZ.EQ.  8 >GO  TO  999 


IFIXX.EQ.IIGO  TO  985 


NOOEB»NOD£0»t 


CMSiHll  I 4 


IN00E«N00E8 


1 NLXST-NXCllMNr(ll-l)-<TX<ll-l)-l*(I»-l>»NX(l) 


2 MXIl)*(NTCl»*l»*l»CIX-tl 


1 NX  ( II  ♦! 


95#  CONTINUE 


XFIXV.EQ. 1.8X0.  «X.LT.IO.O«.IX.GT.X2I)  SO  TO  951 


951  CONTINUE 


IFtIX.EO.il  IXOOE»NL8ST-(XUI*INr{l> -I  > -<NXC11-U-1M  IY-1I»NX  (II 


IFUX.EO.l.AiO.IY.EQ.lir'IOOE* 


ifcix.  eq.nxpi  inooe* 


Z NX(21»(NY(Z)  »11»NX(3> 


1111  CONTINUE 


IFIKSYNHY.NE.Xl  60  TO  11)9 


IS*IS*1 


NYP*NY(IS1  *1 


DO  1180  IX*1,NXP 


IFUX.EQ.l.OA.dX.EQ.NX’llGO  TO  1185 


NOOEO*NOOC»*t 


INODE  -NOOEO 


IF(IX.EQ.11IN0DE»N00EU*1 


IFUX. EQ.NXPI  INOOE  «Nt  AST  • (NX  IS  )-l»  *NY  ( 31  -NX  (2)  »NY  (Z)  ♦l-(IY-l) 


IFIIY.EO.NYP.  ANO.IX.EQ.MP)INOOE»'IX  (1)  »(NY  ( 11  ♦!>♦!»  1 


1180  CONTINUE 


IFIKSYHM.NE.  0)  60  TO  1199 


NYP»NT(IS)*1 


00  1288  IX*  1«  NXP 


IFItX.E0.1163  TO  1289 


NOOEO -NOOEO*! 


ri  u 


INOOC*NOOCO 


IFCNPYLON.EO. 0! PSS»1 . 


IFdX.EQ.il  IN00E*1 


2 (NX(%)-ll*(NT<tMl>  *IX 


1 MX (G  > • (MY  (VI  •!)  • (MX  ( 3>*  1)  *NT  < 31  • 


GO  TO  1291 


IFUX.EQ.il  INOOE*t 


1 NOOEUMIX-U 


lMOOEO*NXm»MT(l»*MX«2»»XT(2)*(MX«J)-t>»MT«  3)« 


1291  CONTINUE 


1211  CONTINUE 


IFUX.EQ.1IC3  TO  1319 


NOOC(>NOOE1»1 


INOOC«NOOEt 


If  (NATION. EQ.  |)P99*1. 


IFUX.EQ.l.ANO.  UT.NE.l.m.ir.NE.lTPI  I INOOE* 


IF (IX .EQ.1.1NO.  IT.EQ.tlt  900E* 


IF  IIX.  «T.  I.  MO.  IT. EQ.lt  I 900E* 


2 NXCl)»NT(ll*IIX-ll 


1 Nt»ST-NXCU»<NTa>-l>-N<(9»»CNTC9>-ll  - 


IF(IX.6T.t.M0.IT.EQ.NT»l  INOOE* 


2 INX«6l-U»CNTUI-l>-<t*t 3I-U*NTIJ1- 


60  TO  1391 


IFCIX.EQ.I.XNO.  (IT.NE.l  0H.IT.NE.9TFM  INOOE* 


2 ♦NX  111 


1 NOOEU*NX(t) 


IFdX.Efl.l. MB.  :».EQ.*r>)  INOOC* 


if>ix.ne.i.*no.it.eq.nt=>inooe« 


139*  CONTINUE 


1 Jo*  CONTINUE 


Nl*ST»NOOE* 


I 


NTP»NT(IS>01 


IFIIT.E0.1I63  TO  10SS 


IF (IX. EQ«NXP>GO  TO  10*9 


10S9  CONTINUE 


IF(PSS.EQ. 1.160  TO  109* 


i NX(E»ocrr-a •nx<ei 


■ j'niirMii 


1 NOOEU*(NX(9)-1>*(NTU>H>*1* 


IFdX.NE.l.SOO.IT.EQ.lH  NODE* 


? NXd>*NTd>MIX-l) 


B<  IHMl'n  1 


t NC*ST-NXd>»(NT(r> -!>-«<  (1>MNT(1>-1>  - 


1 NL»ST-NX(Z)0(NT(Z»-l>-4<  (l)MNT(l)  -1)  - 


IFUX.EQ.NXP>  INOOC" 


GO  TO  109* 


IF (IT.EQ. 1) INOOEa 


IF  (IT  . EQ.NTP1  INOOEa 


IFTIX.EQ.NXP>  IMOOC* 


1F(1X • EQ.l) INOOE* 


tv  it  n » * fin  * 


irjp  icfi>  iPiErrfVJ 


nooeb»moo€b+i 


aim — continue 

INOO£*NOO£t 

uLm.ia.niMOOK- 

1 NB00Y»IX 

t«m.ca»MYmMonn 

1 N800YMNX(IS)»1)*«Y<IS  *IX 

NOgCHlIT.X««XNOOE 

17a»  CONTINUE 

axm — continue 

c 

nlastnuookb 

IFIHHUe.CQ.ltCO  TO  1888 

»!»*» 

NXP«NXUS)*1 

amwui 

oo  taao  iy«i,nyp 

no  laaa  ixbi.mxp 

ixy*ix«(iy-u  »nxp 

< 

if  a* .EQ.i.i'to.  it.st. t >■, d to  taaa 

ttftKa»ttt>06a»t 

taaa  continue 

XMOaCaMOOU 

if<ix.eq.i>hode«nlasth 

Moacxxxxrinamoot 

taaa  continue 

AM* — CONTINUE 

c 

M.ACT«»oo«a 

IF <KS  YNNI»NE*  a)  GO  TO  1TJ8 

ixiMMua.«a.a>Go  to  taaa 

IS«It*l 

NXfcM>x«xw*a 

NYP«NY<IS>  *1 

no  xaaa  iy«i,ny8 

oo  i«aa  ixm.nxp 

XMWftUaMl  KM 

c 

x«<xa.«a.t.*Na.xa.GT.4A.x  to  ioaa 

iFtls.eo.Nxnco  ro  na» 

Nonfanonraai1 

taaa  continue 

UUObJlQOU 

IFCIX.E0.1)H00E*N800YP 

1 INfif  > Altai  NXtli  »Al)N 

l <NX118)U>»<NYU8)-1>* 

i aiu»fi«»inmtiu 

IF«X.EQ.NXP)  INOOE-NBOOT* 

1 mail  Hall  TINT  1 1H  a IX* 

t CNxita>*t>*(Nm»>-i>» 

XIMlIXUUMXI.UlNIlffl  

NOFCT  f XXTt  IS)  ■XNOOE 

iiu continue : 

taaa  continue 
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NYP«NY  (IS»*1 


hum-J  1 1 1 


00  2««I  IX«i,NXP 


IF<IY.E0.1.0*.IY.Ea.NYPl  iO  TO  Z«S9 


70*5  CONTINUE 


8 IX 


imx.8Q.wxa.  »>ta.ix.i«.tnM8P«-maB*» 

1 (NX(15>»1>*<NTU5>«1>8 

I WUWHWMHWHM- 

3 NX(17»»(NTC17»M>*1* 

8 <MX<t»-8>M»XU»>»<« 

XFUX.EQ.NXP)  IMOOf  «N800»* 

1 tMXM.51  tUNMT  IIIHIH 

7 (NX(16|81)»(NTU1>-1>* 

1 MXU71XtM»<U71*U8X8 

8 <NX«  1«)-I)*(HX<HI  *11  *t  ► 

» miuudi^iiwiui 

IFtlX.EQ.NXP.XNO.IT.EQ.MTPlINOOE* 
i MBnnx* 

7 <NX<lS)8l>*tMTfl5>M>  * 

» AMXU»8lM<MXtm»-U8 

8 MX(ir»*<MXCI71*l»*l 

ixax.8o.Mxa. «nn.ix.8Q.i>  inept- 

l NBOOT* 

t-ouumtuonutmu 

J <NX<  16)*1>  »<NT  (16>-l>  ♦*<  (17) 

8-*l 

NOFCT<IXY,IS1«INOOE 

-U»4 — COMTIMUC 

7199  CONTINUE 


NLASTaNOOEO 

IXtM9MAMK.E0.XIC0  TO-  ?»» 

IS»IS*1 

MXB«MX  (I9>  81 

NTP*NT  (IS)  81 

no. 7 7X1  IX»1>MTP 

DO  7780  IX»1,NXP 

IXX»IX.8tIX«HTMXg 

C 

if iix.Fo.Mmcn  in  wax 

N00E 0 ■NODE 0*1 

77X5 CONTINUE 

INOOEaNOOEO 

IXIIT  .f.Q.MXBl  lM0nriMt,A9T8lX 

N0FCT  tlXXf  IS)  8IN00E 

I2U — COMTIMIIE 

7799  CONTINUE 

£ 

NLXSTsNOOEO 

IXtMacA0C.EQ.X160  TQ-33H 

I$8tS»t 

wtxMMuaui 

NX P»NT  «IS)*1 

no  »axa  itmi.mxb 

00  73X1  IX-t'NXP 

IX9»IX>IIX«.U.XMKa 

c 

MoncaMMoocCAi 

1 NOTE* NODE  0 

mopct  »rmrT  to  »mnnr 

73X0  CONTINUE 
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7399  CONTINUE 

£ 

Nl«ST*NOOEB 

It(KSXMBI.B£.0>  60  TO.  2»B» 

IFINOLMJE.EQ.  8)  GO  TO  29  9 

»■»»! 

NXP«NX(IS)*1 

NXPOUUTS)  tl 

00  2988  I T-l,  NYP 

DO  2MB  I«»1,NXP 

ixr»ix*irr-ii  »nxp 


IFCIX.EQ.l.OX.IX.EQ.NXR  JO  TO  2985 

IFI  TT .Efl.NTP)  fin  Tfl  2M9 

NOOC8*NODE8*l 

?U.l — continue 

INOOE*NOOE9 

TEflX.En.niNflOEi 

1 NL»ST-NXP*NTN*in-l>»l«» 

1—U 

IF(IX.EQ.NXP)lNOOE*NLAST-NXP*NTP»t»(ir-l)  »NXP  NX  (IS) 

IFirT-.FO.NTB)  TNOnFiNI  OST-NXPm 

NOFCT  ( I XT,  IS)  »INOD€ 

rnuTTunr _ 

7999  CONTINUE 

£ 

NL*ST»NOOE8 

ntenbnnooeb 

IFINSH8NK.E0.  0)  GO  TO  79)  9 

ISaltoi  ■ 

NXP>NX(I$}»1 

NXPnNX  (IS)  *1 

00  7949  IT«1,NTP 

no  7900  IXal.MXP 

IXT»IX*(IT-1)  »NXP 

r 

IF(IT.EO.NTP)  GO  TO  2989 

MnnKB.Nonra*) 

CC 

-2503 — CONI  INI  IF 

INOOE«NOOE8 

IEIIX.EO.NXB)  TuntWiNlXgoIX 

NOFCT  (IXXt  IS)  »IN00E 

-2S0B CONTINIIF 

7999  CONTINUE 

-C 

HLtST'HOOEf 

TtINBI  il)g.in.fl)CO  TO  2000 

IS-ISM 

wxPNMomitt 

NYP»NT(IS>»1 

nn  kh  tv.<|M»p 

00  2889  IX*lt NXP 

IXX»IXO(IX  .1)  fNXB 

C 

nnngfl«Mnnra.t» 

INOOE*NOO£B 
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NOFCTim,IS»  -INOOE 


canimut 

2999  CONTINUE 

-C . 

NIASTaNOOKt 

W1W9t.AOE.Ea.0160  TO  »1B» 

is«is»t 

ufciuuaui 

NVP*Nr  usm 

op.  »t«»  

00  2799  1X«1,NXP 

HBHIMIMl  ANEW 

c 

wux. ea^i.oa. it. ta.Nxa  ;o  to  2799 

IFIIT.EQ.HYPIGO  TO  2799 

NOOEQaHOOEONt 

2795  CONTINUE 

INOOEaUOOEO 

IF(IX.eQ.l>IWOOE»NLAST-NXP»NYP*<IY-l>*NXP 

1— 

IFCIX.EQ.NXPI  IN00E=NLAST-NXP*NTP»1*«IY-1)»NXP*NX(IS> 

WUT.fiQ.NTBI  IN006iUL1g-WXB*IX 

NOFCT < IXY.IS)  «INOOE 

-2-ZAJ CONTINUE 

2799  CONTINUE 

NAA9T«Nonra - - 

NOOE-NLAST 

■CC 


IF1KSYHNY.  £Q.  9)  NSYNNYat 

lEUSTHNY.NC.  IDMSTWMYa? 

IF  (KSYM1I.E0.  0 INSY1W?al 

■ W-USXWIU.WE.  9>NIYXN2«2 — 

WNOOSaMCOEO 

WRITE  (9.9900) 

9990  FOOWA  T l/tOX.EUOREBROI 

WRITE  <9.99911  NNOOE 

• 891  EORHAT  (OX.EMNNQOEa  ,19) 

00  100  IS*1,YS 

NX 0«MX  <!«■*!  _____ 

NYPaNT(ISI»l 

no  109  IYnI.WTB 

00  199  IX*t.WXP 

IX.Y«IX»-<IY-1>  tMXO 

WRITE  C 9 . 290)  IS  » IX  , IT  , IXY  , NOFCT  < IXT,  IS ) 

199  CONTINUE 

209  F0RNAT(2X,*IS«*tU,2X,PiXa>,I9r2X,>tTa*,I5>2X>«IXT<«,X9.2X> 
1 *MOff.T»T,in> 


SUBROUTINE  GEONET  (HELEN,  NXYHP, NOFCT  ,NOO£,  XK.NTNOOE,  KMAKE,  KROTORI 


COMMON/  ZZZ?/NSFX,NS30OYf  BS,NT  <3%>  .ISFACE  { 39>,KNORHM  39)  ,KNAKES<  J9I 


C OMHON/  ZZZ5/X  PY  CTR  , TP  YCM  , ZPYC  TR,  RX  PY  l ,R  TPYL  ,R7PYl 


COMM  ON/ ZZZ7/Y  NOSE,  YB01,f)0?,YTAIL 


COMHON/ZZZ9/RY0O1,  RZSOl, T YBOZ,  RZBO? 


COMM  ON/ZZZ 11/  X HU8CR,  THU)  I R , ZHUBCR.  R XHU8.  RTHU8  ,RZMU9 


C0HN0N/ZZZ13/THET7S,  THETICtTHE  TiS, CONING , AZIMUTH 


COMH ON/ZZZ 15/ ¥SPAN,XLEZV,XTEZV»TANLEY,TANTE¥,TAU*  ZVT III 


1 f ) rli 


C0MN0N/ZZZ17/NAKIPY,NAK-  B 0 , MAKLSK,  MAKL  09 


COMMON/  ZZZ19/  KPRINT  ( 1 0 It  BREAD, NNRITE,  KREAO 


C OMHON/ ZZZ21/KPY1,  KPT? 


47Vi  ii  Ji  4.  lx 


COMNON/ZZZZ3/K9OEIE,KBO3BP,K0OTY<> 


1 ill  r**#r  J T J 44  1 


CONN  ON/ZZZ  Z5/KPYELE,  KPT5  BP,  KPTTT* 
C 

C OMN ON/ZZZ 27/ KSH£LE,KSHS BP, KSHTY9 


IBf.l.  TtT'WBBMBJ  .TM  4 I i 41-M  4 1 


C ONNON/ ZZZZ9/ KSKEL E , KSKS BP  , KSKTT» 


liFi.  ri1'F*##4i(i] 


COMNON/ZZZJl/NSTAG.NTORr.NSPIRAL,  SPIRAL 


DIMENSION  XK(3,NTN0DE>  ,UFCT(NXYMP,  39)  ,KNAKE(NELEN)  ,NOOE (A.nEIEh) 


NXXsNX (IS) 


NTT-NTIIS) 


00  999  IT*1,NTT 


IF(KNORHUIS)  .EQ.-1IC0  TO  98ft 


I X MM* I X 


UUMttt 

IXHP-IX 

XXBB»IX»1 

IYNN*IY 

IM»t» 

IYPP«IY*1 


c 

— iuum  ...  - - — - 

Tfl  ft«*» 

906 

CONTINUE 

c 

• + — 

r 

**  *•  - - - - 

c 

TO*MnT» 

ixNP«rx 

'«w«r«»  - - 

IXPN-IXM 

imumi 

IYNP-IY 


Q«7 

XYPH*IYn 

c 

TXYPP»rxpP»frYPP»f i>ytf 

rxYPH*ixPH^crr«>H-i)»MXp 

nfY<ifl»t¥MP*  ITVMP-J  1 »y»p _ 

IXTNP*»iXNN*(tTHN-U»NXP 

NODE  ( If  I£LEM)  *N0FCT C IXr  » 1 1 S) 

MfinF  19 tm  PM  .wnrrrr  y*Y-  mv  jo 

MQOE  (3  yZCLCKI  *M0PCT  C IXYi  ? \ IS) 
t *rrrif m 

KMAKC ( I ELEN) * 0 

IF(KMAKES(IS)  • GT#0*ANO«»  X«EQ«NXX)  JCMAKc  <IELEM)*t 
rrr^0nTAB«<|«i.mr.ni  r».  rne  t jri  f«i  * \ 

iZHt  16*50)  IS*  IELEN*  NOtT  ( 1 * IELEN) .NODE  (?,  IELEN)  , NODE  ( 3*  IELEN)  , 
« Mflnrik.Tfi  cm  .yuayorg  :m  . ronna  t rn  f*> 

31 

FOKNAT  (5X»*IS***I5*£X**TEL£N*#»I5*2X**N00Ei**tl5f  ?Xf  •NODE?*  **I  *, 

c 

q«M  _ 

AAMTTIUIP  _ 

1999 

£ 

CONTINUE 

MRITE (6,6*1 

RETURN 

END 
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c 


COMMON/ ZZZ 1/NX  < »*)  ,NY<  J.  I , NXT  « J%>  , KSTNMY , KSTMNZ,  NSTNNT,  NSTNMZ 


COMMON/ ZZZS/NPYLON, NBOWi  , NS OOY2, NS OOYJ,N*TAIt.tNSHAFT,NMU0,NSHANK, 


COMMON/ZZZN/UNACH.OHSSA.  tlFMHTI 


17d  i HT17y*1 


:onmon/zzz*/xnose,xsoi,«)02.xtaii 


COMMON/ZZZS/ZNOSEt  ZBOl  tl  SOZ.ZTAIl 


COMMON/ ZZZ 10/  RSHAFT  t LSM*  ?T  t RSHANK  * l SM  A M< 


C0NH0N/ZZZ12/  RROTOR*  DCHlROtTAUBltALFAD 


F* rm  M 


COMM  ON/ ZZZ  IV  KBIAOE*  TA)t.i8iTANTEB*XBL£»X9TE , KROTORS  C 3A) 


C OMriOM/  ZZZ 1 6/  NW AKP  Yf  MNA<  49  » NWA  KSKt  N WAK9L 


COMMON/ ZZZlS/WANGPr.NAIC  4B»WANGSK*WANGBL 


COMM  ON/ ZZZ  20/ FI 


C0MM0N/ZZZ22/ KNSELE  , KNS?  4 P , XNSTYP 


T|T ll 


C ONNON/  ZZZ2A/  KT  nELE  , KTNS  4P , KTNTTP 


COMMON/ZZZ  26/KVTEl  Ef  KVTS4P  f KVTTYP 


C ONNON/ ZZZ28/ XHBELE*  KMBSMP , KW TT» 


riT  r*>i*i:ri  « il a 


COMNON/ZZZJO/KBLEI.E,KBLS4P,KBl.TYP 


DIMENSION  XKCS.NMOOE),  W0E<«.,NELEN> 


■ ■;  i * 1-4  ■ fy  l 


MEAL  L SHAFT SHANK 


ALFAR*ALFA*3. 1%199/lBO* 


cosalf*cos<alfar> 


XK1»XK(1. INODE) 


I 4 fAMI 


XK  (1,  INODE  >*(  XX  1»C0SALF*XK3*SINALF>  /BETA 


WRITE  16*61)  INOOC»X  XI  fXK3 


10*  CONTINUE 

mi 


XKl-XKU,  INODE) 


SUBROUTINE  PRINT 


SUBROUTINE  PRI  NTA  (NELEH,  NXYMP.XK.NNODE.PC.NOFCT.NOOE  ,NPRINT) 


ommon/zzZ2/nsfx,nsbooy.  <is,nt(39),isface(  j9),knorml<  39>,knakescj9> 


4»rrT.T'r***Ji»r 


COMMON/ZZZS/XPYCTR.TPYCYR.ZPYCTR.RXPYL.RYPYL.RZPYL 


OMMON/ZZZ7/YNOSE,YB01,M02,YTAI. 


COHNON/  ZZZ9/RYB01,  RZBOl,  RYB02,RZ302 


Tirfl##* 


C0MM0N/ZZZ11/XHU8C  Rt  YHU9SR  , ZHUBCR,  R XHU  B,  RYHUB.RZHU3 


C0MM0N/ZZZ13/THET7S,THETIC,  THETIS,  CONI  NS,  AZIMUTH 


I ' TiT’»*  ***■ 


COHMON/ZZZ1S/  YSPAN,XLEZV,XTEZV,TA*LEY,TANTEV,TAU,ZVTAIL 


■mVT.i;  r»i  'TBmmm*  ■ 


conmon/zzziz/maklpy,mak  tb , naklsk, nakl bl 


COMMON/ ZZZ19/  KPRINT  ( 101,  TREAD, NMRITE,  KREAO 


COMMON/ZZZZ1/KPY1.KPYZ 


■ MtKI-ltl’r#**' ■*•>■**  47V» : 


COMMON/ZZZZ3/X9DELE,KBOHP,XBOTY» 


Uil  I * 


C0MM0N/ZZZZ9/KPYELE,  KPYS  RP,  XFYTYP 


C0MN0N/ZZZ27/KSHElE,KSWRP,KSHTYP 


COMMON/ ZZZ29/KSKELE,KSIG  HP,  KSKTYP 


C0MM0N/ZZZ31/NSTAG,NV0RT,NSPIRAL,SPIRAL 


real  L SHAFT, taSHANX 


GO  TO  (1,2, 1,9) ,NPRINT 


MRITECNNRITE,  111) 


00  112  II*l, NS 


Nin  iliiil 


NRITE  ( NHRITE,  1 1 3)  ISF  ACE;  1 1 ) 


IT1  ■ i -M  I I I-)  • 


112  CONTINUE 


NRITE  (NNRITE,  115)  NELEN,<3YMHY , KSY9MZ,  REFLEN.SPAN,  TAU  , 


119  F0RNAT(2X,»NX«»,I2,/2X,»1T«»,I2/) 


1 2X,*REFERENJE  LENGTH  »• , F6.2/2X,»SPAN  «»,F6.2/ 


1 ?X, *NACH  NUMBER  ••,F7.»//2X,*IZZZ»»,  I5//> 


116  F0RMAT(2X,»TANGLE«,,F6.»/2X,»TANSTE*»,F6.2//2X,*CH3R0  «»,F6.2 


2210 

221 

WRITE  (NWRITE,  221)rN00E,(  <K(K, INODE)  fK*l#  3) 

FORMA  T (4XfTL,3M»rl'l  1 

?z« 

CONTINUE 

BP  TUB  It 

3 

continue 

MBJPf  ( MlIB  J TF  r tkAl 

JM 

F0RNAT<////lX,*ELEM«*f  M,  •XPC*.7X,*YPC*»7X#»ZPC*I 
no  Tkk  T«1  fMPl  PM 

344 

MRITEINMRITE,  Jn5»I,  (PC  «,  I),K*1,J) 

tkl 

CO  DMA  T t 1 Xrl  Y*  t 7Ci  fL‘S* 

RETURN 

kf 

WRITE (NNRXTE,  41) 

FQfiMAT  * MAflAl  Ml  IMA-'  3 TMt  FAB  *UBr  4 OF  * • V /» 

00  44  XS«1»NS 

uVE.yvqii»i 

NYP«NYCIS> *1 

UVfO.UVPfUYP 

WRITE  <NNRIT£,4  6)  ISF  ACE  (.  3) 

nn  kk  t*«i-m<p 

WRITE  < NWRITE,  45 MN0FCT(  XY,IS>  , IXY«  IX  , NX YP , NXP> 

OOMTfMIlP 

*5 

FORMAT  (1116) 

kA. 

WRITE  CNWRITE*47) 

FORMAT  f//<U..»MOnAl  MUMA*  2TM&  FOB  PMFMl**/ /I 

WRITE CNNRITE, 49 > 

nn  kk  rrir»«<.Mri  m 

48 

NRITE  (NWRITE,<*5)IELEH,  M30E  (I CORNU,  IEL  EM)  , ICORNR*  1, 6 ) 



RETURN 

F MO 
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C0MM0N/ZZZ19/KFRXNT ( 10 >t  9R£ADfNmiTE,KR£AO 

— turrf.OiMairt,  ti  r. 

I FORMAT (/2K,»  ERROR  C03:  • •,!»/» 

— RftUEJI 

CMO 


X AR1BRRRBBA  SUBBOIITTME  3ITHTB  IIIUIIIUIIIIIIIIIIIIIHIIIIIIIIIII 

c 

X 

SUBROUTINE  PAINTBINELEM,  (ESO,  A I.SOURCi  ,NPRINT) 

COMMON/TlXl/MX-HH  , Mill.  ) , MXTdBt.KtYMM*.  KSTMMZ , MET  MM  t , MSXMMZ 

COMHON/ZZZZ/NSFX,  NS800T,  X S , NT ( SMI  ,ISF1CE(3M,  XNORHL  ( 36)  > XNA  XES  < 36  I 

COMM OM/Z7X3/ MBY LON. lOftm  1 .MB 00X1,88 OPTS.  MMTAU.MSBAAT. MMUO, MSN AMK. 

1 HBLAOE 

COMMON/ ZZZfc/UMACM. JMtOA. ALEA.ABEtA 

C0»N0N/ZZZ9/XPYCTR,TPYCTA,ZPYCTR,  AXPTt.  .AYPYL.RZPYL 

coiM/m/mB/mow.  natu.no>.  xim- — 

C0hMON/ZZZ7/YNOSE,YBOl,M02,YTAIL 

commoh/zzz*/zmo«.  zooi.r  ioi.itail 

C0NH0N/ZZZ9mB01,RZa01,  »YB0Z,RZ90Z 

caMMOM/zm«<gsigRT.ctN>iT.asMAMit^tSM>m 

COHHON/ZZZll/XMUaCA.YMUlIR.ZMUBCR.RXNUB,  AYHUJ.RZHU) 

C0MM0N/»zu/BaoToa.8CMT6o,  lAuaurAaCAa 

conmon/zzzi3/thettb,thetic,thetis,:oni.n;,»zimutm 

n0HH0M/7Z7t»/xai Anc,  riM  - fuiiMTEa+xai-: .kite mannas i tli 

C0NH0N/ZZ2H/*SPAN,XLEZV,XTEZV,TAXIEX,  TANTEY,TAU,Z/T»U 

C0MM0M/ZZ7l6/MMAKBT,HHA/<8,HMAt:SX»HMAXlL 

C0HNON/ZZZlT/MAXLPY,MAK_3B,MAKLS<,MAI<LaL 

rriMMOM/mtAf  M«Mf.BT.M4M.<8.MAM6S*.MAMr.a 

CONHON/ZZZH/KPRINMIO).  IRE AO,  NMRITE,  BREAD 

common/ mia/Bi 

COHNON/ZZZZl/XPYl.XPYZ 

rnM«nM/f;m/itM«LE.icMs;4»,KMtTT» 

C0MN0N/ZZZZ3/XB0ELE.  XBtPMP , XBDTT* 

C0MMCM/Z7Z16/XIMEI.E,  KIM'./B  ,XTMTT* 

CONN ON/ZZZZVKP TELE,  XPY33P,XPYTY» 

COMMON/  Z 77  26/ XX  TEL  E,KMTHB,XXTTT7 

COMMON /ZZZZZ/KSHELE* MSN?  3P,KSMTY0 

C OMMQM/ 72228/  MM8ELE , XM81MB  , lf.M8  IT  A 

C0HH0N/ZZZ29/KSKEI.E,KSIGmp,XSXTyP 

C OMMQM /ZZZ3  l/»IHi,ti  XBi  ;<B,K8LTH 

C0HH0N/ZZZ31/NSTAG,NX0Rf , NSPIRAl, SPIRAL 

X 

DIMENSION  S0URCE(NELEH).IA(N£SQI 

MmHCB-aiBAUWMUMWW 

real  l shaft, l shank 
mT6«6»ibX(I>«i> 

60  70(1,2,3,6,9) ,N PRINT 

J COMTIMIIF 

MRITE  (NNRITC,  not 

111  E08BATI/Z/33,  tQlSTRZBUT  M OB  AAlUJAY/) 

00  It  I»i, HELEN 

BBITf  IMMBIH..UIU 

Nl«l 

MF.MPI  PMPMFl  :■ 

111  FORMAT (2X,*IN0£X  »p,I2> 

MBITf  (MMBITF,  1121 IAA  IKI.  X»MI,M  Z.MEi  EMI 

HZ  FORMAT  (IX, (El'S. 61 

II.—  CONTINUE 

RETURN 

X CONTI  MIN 

MRITE  IMMRITE,  ZZ1) 
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221 

225 

FORMAT {///2X, *THE  DISTR  3UTI0N  OF 
CONTINUE 

SOURCE*) 

NSBTOTM 

na  22  4 ISsl.VS 

WRITE (NVKITE, 22  3>  I SF  ACE<  I S > 

2 2i 

IND=NS°T0T 

dS5 TQT-NS9TL LAtiXU-SJ.  »NY  f5) 

IFIN* NSBTOT 

ALXXsNJUIS) 

00  226  IX=1,XXX 

MRTTE  (NHPTTF-  7?«l 

IFINPRINT.GE.5 .ANO. ( IX.EO . NXX) ) GO 
_nm=iNO*j 

TO  *26 

<*RITE  ( NWRITEt  227)  ( SJURC  ( KK ) , KK=I NO  ,1  F IN , NXX) 

-CONTINUE 

222 

FORMAT  <lXf8E15.5) 

-22*i 

FHDMA  T tf\ 

229 

CONTINUE 

RETURN 

3 

CONTINUE 

JARXXE  (NW.°UE,J30) 

330 

FORMAT  (///2Xf  *THE  OISTRiJUTION  Or 
r.n  Tfi  ??*. 

THE  VELOCITY  POTENTIAL*) 

4 

CONTINUE 

-WRITE  t KHRlXE-tA^-DJ 

440 

FORMAT  (///2X,  *THE  DISTR'JUTION  OF 
X-Q-TO  -225 

CP*  > 

5 

CONTINUE 

BFT'JRN - - 

ENO 
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C0HH0N/ZZZ19/KPRINT  ( ltl.  4REA8,  MtUITE,  KHEAO 


S OHHON/ZZZZi/ KPYt, KpY? 

C0M40M/Z77ZZ*«tt«LE.*ltS;4B.rNtTTA 

COHHON/ZZZZJ/ICBO£LE,KBtmP,KBOTY» 


C ONH  OH  / ZZZ  75/  KP  YEL  E , KPT*  4 P , KPYTYP 


[lUJiITf*#**,  TJ  4 i 4jIJ 


iirfl Jl 4 tM 


COHHON/ZZZ27/ KSHELEf  KSH54PfKSMTYP 


W0&J  rn*'  t im 


C OHHON/ZZZZV  KSKEIE , KStG  4P,  KSKTY» 


COMHON/ZZZJt/NSTAG,NYOW  , MSP IR4L,  SPIRAL 


III*.  M 


mjI  f 1 j 1 a .1 


REAL  L SHAFT ,LSH ANK 
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0IMENSION  KtM(l) 


00  S L»l,f 


1 1 Villi 


IF<T8B-PIV>3,J,7 


start  elimination  loop 

I «T_« 


00  IT  K*t,N 


TEST  ON  SINGJLARITV 


A IF<IER>7.5,7 


I Jff  Jt  H ITBlWI 


A IER»K-t 


i i € » ^wm  Jm 


J«<I-1»/M 


DO  A L*Kt NMiN 


TB«PI»I»RILL» 


A R<Li»T8 


IS  ELIMINATION  TERMINATES 


9 LENO»LST»N-K 


00  IS  LsZSTfMM»N 
UH.  J 

A(LI*A(L>«PIVI*A<LL> 
IBB*  1B14KLII 

ir<TM-*tm5»is>t* 


C ENO  OP  ELIMINATION  LOOP 
X 

c 

x PACK  LuatTITlJTUlM  Attn  WIN  INTERS MA MSE 

IS  IP  (H~l)E3f  EE»  19 

1-8  tSIUMMI 

LST«M*1 

no  ai  i»m 

II-LST-I 

ttltllhUT 

L-IST-N 

lUlllt.E 

00  El  J»II,NH,N 



LL«J 

no  n mm.iw.M 

LL*LL  *1 

It  IO»It«AINI»KtLH 

M»J*L 

ouimouo 

El  R<K>-TS 

- EE  RETURN 

c 

Jt 

C CROOK  RETURN 


E3  IER»1 
a£XUUL 

END 
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**********  SUBROUTINE  C1EFF  ft*********************************** 


SUBROUTINE  COEFF(NELEN,=  : , NNOOE.XX.  NOSE,  KNAKE,  NE5Q,  II  .SOURCE, 


CONHON/ZZZZ/NSFX,NStOOY,  IS, NT  (391,1  SF4CE<39>, KNORNt  till, KttlKEStJS) 


Tilill'l  W I IMt.il 


»l  WW’U  11  Jn  TTIII.WTL'T  :F  I’ll 


1 NBLIOE 


rT  #•<  WH  #W  TT  ill 


COHNON/ZZZB/XPYCTR, YPYCT1 ,ZPYC  TR,  RXPYL  ,RYPYL,RZPYL 


C0NN0N/ZZZT/VNoSE,VB01,V10Z(VTIIL 


CONNON/ZZZ9/(YB01,RZB01.  *YBDZ,RZ80Z 


CONNON/Zmi/XHU8CR,YNU»:R,ZHUBCR,RXMU8,RYMUa,RZHUB 


CONNON/ZZZ  IS/  THET79,  THET  tC,  THETIS, CONIKi  , IZIHUTM 


0A ITI  riw  * 


i # * Hi  ^1  III  ;•»  41  lEiifi 


COHNON/ZZZ19/VSPIN,XLEZ/,XTCZ¥,TINLEV,TINTEV,TIU,Z«TIIL 


CONNON/ZZZ  1Z/NIKLPY,  NIK.  I8,NAKLSK,NAKL  8L 


CONNQN/ZZZIB/KPRINT  (181,HREIQ,MNRITC,KREIQ 

r » /immau  / T9T a* 


C ONNON/ ZZZ  Zl/  KPY1 , KP  YZ 


C0NN0N/ZZZZS/K80ELE,  K905  IP,K80TYa 


C0NN0N/ZZZZ9/KPYELE,KPVS-lP,KPVTY» 


4 e*  riri  JW  Hn  4-J 


CONNON/ZZZ Z7/KSHELE,  KSW  IP,  K»TY> 


*»*  I'l  i 


CONNON/ZZZZ  VXSXELE,  KSKJ  IP  ,KSKTYP 


C ONNON/ ZZZSt/ NSTIG , NYORT , NSP IRIC (SPIRAL 

r.mmnM/zzzsz/csTiG.cYtiRr 

C ONNON/ ZZZ33/ UN  IKE 

nrMKKTOM  mtifumci  rat 

OINENSION  QVISt 

. OINENSION  XKtT.NMOOEl  (NflJEtfc.NEURt  , 

1PC<  J.NELEH)  (SOURCE (NELEH  .KNIKEINELEN)  .AIINESQ) 

-OUgM-SION  H(3.?I.I?t3(?l  .IIIUTUAZIZ<7U0C«H<3),0r.BlZt31 

OINENSION  111  (3) ,IZV  (SI  ( / NORN  (31  ( UN  (31  ,SNCS»,SNUNm  .IIIZIZ.Z) 
OINENSION  RZC  II  (0«S,N»>aiCgMI3,Z>,0Nr.RIZI.3tZl,IIUS.ZI,/UZ<3>Z> 
OINENSION  YNi  (3  ),  VOOC31*  YNO  (31  , YTEPZC  31.  YTENZI3)  > SISNPT(S) 

OINENSION  l¥lt(  1)  ,IMII<3>  .I1CRIZOI  ,QNII,I) 

OINENSION  NP<>  131  ,NPN  (31,  4NP  (31  ,MN<(  I 31  , XPC I 31 

— OINENSION  Br.IMFLFN I 

OINENSION  YP*  (31,  VPN  (31,  TNP(3t  ,VNN(  3> 
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DIMENSION  PCJIJI 
» EU.  L SMAfT  SMAMK 


aaiwoiiMMMi.n,^  .uMMumam 


a>wxiHi.wi,Hi,na,rn,m,iin,>»i,imMinmiM.>»ini)iwi 

ZZZ>»YYl*(XX2»TYJ-XXJ»YTZ)»m 

iBaanni.nr.nz.tY^Trfzi  ■ otyez.ctyoi _ 

xproy  (ox,oy,iz,ex,ey,ed  «oz»ex-ez»ox 

»am  ia»,BMi,«,»,i>  »a»«a»H«B» 

It(UNACM.6I.A.  


9ETA«SQRTC1.-UNACM»»ZI 


sircar  mu  .a 


XOUNT »tt 

aa  iax  jiUMfifN 


KStUALat. 


(.SPIRAL  *1 

jxuaaaut^Mi  ksumau  mm. 


IP  IICMECX.EQ.V  > ISPIRAL'  YSPIRAl 

oo  ta 4.«a»ir«saiaAt. 


00  It  I NSP«1, ISPIRAL 


eo  ro  (ia#,z#»,  Jto. rcMEcx 


CONTINUE 


00  9<»  KM,J 

«uutoBUM>attMu 


YPM<K)*XK(K,Y00£(2, J)» 

xjiau  > «xmN,»oor  <!>■»»> 


YMM(K I «XK(K>N0  0E (•»,.!>> 

acjonacis.jt 


CONTINUE 

r.n  to  saa 


ZtO  CONTINUE 
NST.AB«I 


IFOCOUNT.EO.l)  return 
af«Q<MBEAft>ZOI.>  (TORCH. 


NRITE (NNRITEt  ZBZICYPPOO  , K«l,3> 

aumuoiMM untua  ,oi,ii 


NRITE < NNRITC*  Z»Z)  (YPNOO  , 
araniNflfAn»3ai) >tju»io 


NRITE ( NNRITEf  2B2)  (TNP(I0  , N»l,  J> 
arnuktnruzan — (tmmiik  ,k«i,d 


NRITE (NNRITEt Z«Z> (TMN(K 

afl-ZBJ  ui»i 


ZBS  PCJ( X ) ■ ( YPPCXI  ♦ YPNOCI  PYX  POO  ♦YHN(K)  ) / % * 
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WRITE  INWRITE,  282)  ITNPC  X)  , K*l,)l 


WRITE (NNRITE.  2021 (VMNTK  ,K»1,I> 

AM  V«l  Wm  4 * 


30$  PCJ(*l»«TPP(<l*TPN(K)"T'«MK>»TNl||Kll/‘., 


201  FORNATI  10*0.3) 


SO  TO  500 


N SPIRAL*! 


IFCKROTORI  Jl.EQ.OISO  TO  180 

nn  LAI  IF. « ■ 


TPP(K)»XK<KtN00E(l,J)) 


14j.II  4 1#  i <1 


TWFIKt  ■XK(K>M0  0C(3>  Jl» 


VI  4*0  <414 


A 9 A CONTINUE 


RATI0*UWAKE/}NE6A 


rl  >*l«Tl'Jl4iiVl*'l 


RATI0*-1. "RATIO 


liWliT  nUJ  ifi  U il  M 


00  180  lSP«l,ISPIRAl 


IWAKE«NSPa(ISP-II"NSPIR». 


THCT  A 1«  II  WAKE- 1)  "FACTOR 


COSHCOS  (THETA  11 


TPP2»-TPP(  1>"3IN2"TPP(2>  *C0S2 
XNP1*  TPPI1I"C0S1*TPP(2*  "SIN1 


TNN1>-TPH(1>"SIN1"TPN(21  "Cost 
XPM2«  TPH(1>"C0S2»TPNI21  "SIN2 
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TPH?»-TPN«1»»SIH**TPH»?»  »COS? 


RPP?»SQRTCXPP?«XPPZ*YPP?»TPP?> 

aBBl«$QBT(XPBl*XNPl»TBPl»»NBl> 


RNHi>saRTtxNNi«xNNi»YNM*TNHii 

aBn3«tw<xami«pm»»pm»vpm> 


if  msB.sT.ti.  so  ro- $oi 


ZPP?«ZP*aOTIO»TM£TO?»RP=?/RROTOR 

T»auM*a*Tn»TaBx«tmfci/aaoxoa 


ZHN1»ZP*R0TI  0*THET0  1»«5N  1/RR0T0R 
ZPHZaZEfam»THCTtZPgPlZ/aaOTOB 


so  to  $oa 

r.ai  continuo 


ZPP?«ZP*ROTIO*<THETOZ-n:TOR>*ROTIO»f»CTOR»RPP?/RROTOR 

tufli.XBAaxTxati  TM£Tii.ia:ioa>-*a«riotf4CToatRMBwaaoToa 


ZMNl»ZP«-RBTlO»»THETOl-r.CTOR»  *00110  •P*CTOR*RmU/RROTOR 

IBB? « ztoao  rxa»  i TME.iAz.f^iroauaaria»f  acro*tBPHB/^BOT  oa 

$04  CONTINUE 

TPP(1)»XPPZ 

XBBUUtBBl 


TPPCJ I *ZPPZ 


TPNC1 )*XPM? 
tMHUtBU- 


TPN1JWPN? 


YMPUJ-XHPl 
tmpi  ?utmpi 


YHP15WNP1 


TNN { 1 ) *XNN1 
mil  »«»“■» 


VNH10I *ZHN1 


IF(ZPPZ.LT.(ZB01*RZB01»- ZPYLDCO  TO  la  0 

if (ZBBQ.LT.izaot*azani»iZB'ic>*-  so  to  ua 


00  BQB.,X«1.X- 


50$  PCJlBI*«TPP«)*rpN(X>*T-  P(K)*TNN(KI  1/4. 


SOB  CONTINUE 


00  f BO  X»trl 

001CK,U«0.5MYPP(M-TN>>(ICI> 

00?IK,l>»B.S*tYPP«K>-W<K>» 

uttMUMUMim  msojua 

*1(K, 1>>001(K,1> 

01  IX.?>«001  IX, ?> 

0?<K.1)*00?«,1I 



100?  CONTINUE 
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1003  A2A2(L>*00TPT0(A2(l,l>  .A!  <2 ,U  , A2<3  , L)  ,A2(1,U , AZIZ,  L > . A2(3,L)  > 


00  II (A  K>1,3 

»<AlAlMl.fa.a.lA4(K,4>  iAAIK.31 

IF  (A1A1(2>  .E3.0.  >A1(K,  a «A1  <K,  II 

If  (A2A2IH.f3.il. >A7(K,J»  «l?(l(. 21 

IF (A2A2(2)  .EQ.(.>A2(K,a  >A2(K,  1) 

uuu — r.nxriHUF 

c 

an  (aqa  lh,? 

A1A1<L>»00TM0(A1(1,L>,A1(2,L>,AH3,D  ,A1  ( 1,U  ,A1  (2,  L I ,A1  (3,  L)) 
1016 1717(Ll«nOTP7niA7H  .1  .A7Hrl  I I.I7H.I  II 

c 

do.  iiai  i*i,2 

00  1117  N-1,2 

1»AT  A1AI(l,llUnnWOAA.l<l»A.l,Al<a>U,ti(3.  U ■ Aail.io.  A3<2,M>.A2(3,MH 

c 

AlCaA7(1USMT(AlAl(71t7A?Mlf  A(  A?  (7.11*771 

AlCRA2(2)>SQU(AlAl(l>*a  2A2(1)-A1A2  (1, 1)  **2> 

umAuusauaiAuimaiuuAiuiMUAa 

A1CRA2«,»>SQRT(A1A1(2)*1ZA2(2)-A1A2  (2,2)  *»2> 

C 

00  Ilia  K>1,3 

AKA1  I1HI.S’(A1  (K.H.AHn.21) 

iaaa  ava2(io*o.5*(A7(k,i> *A2(<,2>> 

C 

VNORN  (1)*AVA1(2)*AVA2(3  - A VA1  ( 3)  *A¥  A2  ( 2) 

uflAmatAAimawaAut  oaAumwwi 

YN0RH(3)aAVAl(l)*AVA2(2l  - AVA1  ( 2)  *A»  A2  < 1) 

C 

SN  ( 1)  ■ YN0RK1 ) /BETA 

SM(2)»YMf)BN(2) 

SN(3)*YN0RH<3) 

A<iM»^Qg,itfii(nTf?«m(2)fr  nauimi 

AYN»SQRT(Y NORMt  1)**2*Y!OTH(2)**2*7NORH(3)**Z) 

on  tm  K»if3 

UN(K>*YN0RM(K>/AYN 

3M»»»)ilA(«l/AW 

iaia  CONTINUE 

c 

c 

no  17a  i«i»NfiFw 

c 

an  1AA  ISTNUral.MTTMNY 

00  161  ISYNNZal.NSYNNZ 


Q(K,  J)*TMP(K)-PC<KfIl»ST  GNPT(K) 

atn,h»»Tmaic>  -pc  iK.  ix±sl  mmi 

lilt  CONTINUE 

00  MU  K»l,.l 

QM<Ktl)*0.9«(q<K>l)»a<K,  »>) 

BMU.  2>»i.8T<Q(«.2>4fl<fc.  I>> 

QN(K>3>>S.5*(Q(K,3)»Q<K.4)> 

HHnt.i>i»n.^T<Qnt(i>i«a<K.t)> 

till  continue 

X 

00  1112  K«t>3 

KB1»FM 

KP2«K»2 

IF<*P1.ST.TUB1»KB1«-I 

IF«KP2.GT.3XP2«KP2-3 

QHCB11  <Ktl>»aN(ICBl,21t«llKB2.,l).QA>ltP2,2>.FAl.<*Bl.  11  - 

0NCRAl<K.2l*3N(KPl,4>*Al(KPZ,Zi-3i(KP2,4>*AlfKPl,7> 

3NCBA  7 < If  » 11*04  I KB1, 1 1*17  t ICB2, 1 > -,1Mt  KB7..1  I »» NBA . U - 

QHCRA2(K,2>«3M(KPlf3»*A>  ( KP2,2)  >04  ( KP2  >3)  *12<KP1' 21 
rnurr  mir 

C 

RC*S3RT<DOT«01BI<l>->B;t<21  ,HI  II  ,BZ  <1  > ,B7<  2>  ,BI(  3>  1 1 

no  n*  TC0Bta»i,t 

SO  TO  (5582.5584,5586, 5588). ICORTT 

CONTTNIIF 

SI6N12—1. 

IC5I«1 

IETA-2 

60  TO  Till 

5504  CONTINUE 

5ICM17»*1. 

ICSI«1 

IETi«* 

SO  TO  5510 

1U1 — CONTINUE 

5ICN12  — 1. 

1C»U 

IET*«1 

CO  TO  EH  a 

5508  CONTINUE 

IICNAB-*!.. ... 

ICSI*2 

IETA»? 

5518  CONTINUE 

X 

00  5528  K»l,3 

q8lltl»0«».rCOBNB>  

A1V(KI*A1(K,IETA> 

A7WIICl4*2tK»tCSH 

QCRAMKMQNCTU  (K.IET8) 

Qcaiz  (r,  ir.su 

5528  CONTINUE 

X 


QQ«SQRT(OOTP*0<QV(1)  ,QV(2>  ,QV(3>,QV<1>  ,QV  < 2) ,Q« (3 > ) > 
CALL  L0G(IC0RNR,QQ,QV,A1Y,QCRA1, 11061.  U 


c 


IFIQOOTUN.EQ.6.IGO  TO  T- 5 » 


i f (oe  noh. he.  a.  )t»np**t*'(>  (hmuher.be  noh) 


FiilHiTwr  VMt-lM 


QCRA2<3>> 


COEFF 2*00TPRB  ( UK( 1 ) , UN  <21  . UN  ( 3 ) .BCRA2  1 1)  ,QCRA2(2)  f3CPA2<3>> 


SRCINT*-CONST*SIGN12* 


691  CONC-  1.6 


862  CONOCSTAG 


6 63  CONOCVORT 


60%  CONC* 1.0 


0BTINT*-C0NST*SIGN12*TA,P*C0NC 


SGNINT*!. 


nnn«xmj-d*nelen 


AA(NNN>*AA(NNN> -SGNINT  «J8  TINT 


UFREE*UHACN*111Z. 


0ME6AY*0« 


RX*PC (1,J> 


R7*PC(3,J> 


8CY*X  PROYIOHEGAX,  ONE  GAY,  ONEGA  Z.RX.RY.RZ) 


8CR*OOTPRO  (OCX.  8CY  ,BCZ,JX(1),UN<2),UN(3>> 


8CB"-l.*UFR£E*UN(li 


l I j-iii  ffirlY'f  \ WirWr*  \ i 


198  CONTINUE 


186  CONTINUE 


SUBROUTINE  LDG(ICORNR.(a,  QV,AV,QCRA  V,  »L3C,  ICOORI 

aimsiM  anti)  

OIHENSION  AVI  31  fQCRA  V<3) 


00TPR0<X1,Y1,Z1»X2.  Y?,Z’I  «X1»X2»Y1»Y2*Z1*Z2 

z 

AVAV-00TPR0IAY11),  AV  < 2> , A V ( 3)  , AY (1 ) ,A  V (2)  , AV(3)  t 

a0AlUnoifB,0IIWHI>QV<2).3M(MtAV(l)  ,A-V-IZ1^-AVH ) ) 

0XA*SQRT  (OOTPRO  (QCRAVC 11  • QCRAV  (2)  , OCR*  V(  J)  , 

1 icmu ti . orjAmzi.arRAVMin 

alog«asinm<oqav/qxa»  /so? t iavavi 

RFTIIRM . 


FUNCTION  ATAMPtMNUNCR,  0E  HOH) 


»T*NP»*T*N?«HNUNCN,*mtNDi) 


i p 1 1# ri  ] 


SUBROUTINE  A/ERACISOURC:  »SOR, PC, HELEN,  NOOE,NOFCT,NKYhF,NNCOE, 

1-AttU 

DIMENSION  SOJRCE  (nELEM),  ? OR(  NNOOEI  , ABC  (NNOOC)  ,PC(3,  NELEMI 

OIWNSIOM  MQ0GU.MfcLSM».H08CTlMMYNB,H> 

COMMON/ ZZIl/MKH)  ,NY(  JVI  , NXY  (»>  , KSYNMV,  KSVNMZ.NSYMMV.NSYMMZ 

r.nMMnM/7777/Msrit.nsannv.  ms-utiimi . r i ,n,r„,nni.i 

COMMON/ ZZZJ/YPYLON,NBOD/l  ,N800Y2,N800YJ,NirTAlL.NSMaFT,NNUB,NSNAMK, 

1 M8UUM 

COMMON/ ZZZL/UMACH, OMEGA  ALFA, ASET* 

cnnMOM<ZT7s/YBYrTB,YBYga,ZBYcni. i«BYi  .AYBtL.azan 

COMMON/ ZZZt /THOSE, YBOl ,1102,  XT  AIL 

cflimnN/zzzz/iNasf.Yani ./  mz.TTin 

COMMON/ZZZI/ZNOSE, Z801.Z 102, 7TAIL 

r ONNON  / T7Z8/BY  801 , B 1801.  axanz  ,82802 

COHMON/ZZZlS/RSMAFT.LSHt  ;T,BSHANK,LSHANK 

irONNnN/Zmt/«NU8CB,YMlN  •’■B.THMBCa.aXMlia,  BYMII8.8ZNIH 

C0HM0N/ZZZ12/RR0T0R,  8CM0AO,  TAU8L,  ALFAR 

COMMON  /ZIZtl/.INETZl,  THET1C, I METIS, CONI NS, AIINUTM 

C0MMON/ZZZ1L/K3LAOE,  TAN.  r 8 , TANTE 9,  « 8L  £ ,XRT£ , KROTORS  < JA1 

COMMON/ 11118/YSBAM.XLtZ/,  XTG28.TAMLG8.  TANTEN,  TAU,Z»T  A XL 

COMMON/  ZZZ18/  NNAKPY , NMA<  -48  , NMA  KS<»  N MA  X8L 

COMMON  / ZTZ1 2/  MAJCLBY , MAN,  AE.MA81  tlf.MAKl.8L 

COMMON/ ZZZU/MAN6PY,  MAN, 4B,MANCS<,  HANG  8L 

COMMON/ 77718/ KBRlMTi  HI.  YBf  AO.MMaiTE.XBEAn 

C0NM0N/ZZZ2B/PI 

COMMOM/77771/KBM.XBT2 

COMMON/  ZZZ22/  RMSElE,  XNS*  /P.XNSYYP 

COMNOM/77171/AaflELC.  ABgj/B.ABOTYa 

COHHON/ZZ22WRTnELE,KTN>  /P.KTNTYR 

CQHMflM/7777VKFTFl  F.  KBYSHB.EPTTTP 

C0MM0N/ZZZ28/LYT£LE,KVTS'tP,KYTTY> 

COMMOM/71177/XTMGlE.XSIC  /B.XGMTXB 

C0IW0N/ZZZ28/XH8ELE,  XNS5HP,KM8  TY° 

~0MM0M/77778/ESan  r.XlXI/B+JtSXTia 

COMMON/  ZZZ38/R9LCL  E , X8LE iP.XBLTYP 

f-niiMnM/7TTYt,MMTtc,imn«r  , ,',r.n 

INTEGER  FM.PP 

RIAL  LSHAEI.LtMAMx: 

C 

00  88  IHOBEmI i NMOOE 

AVGIINOWImS. 

-88 SOB<INCn»>«  a. 

00  188  1ELEN>1, HELEN 

BB«MOM  It.  IELEMI 

PMmMOOC (2, IELEMI 

MByHOOE  IT,  IELEMI 

MM.MOOEC4,  IELEMI 

U8XM1AME081  * SOURCE  11  -LfMI 

S0RIPNI'S0R(PNI»30URCE(.  EL  EMI 

SQAIMBImSORIMBI  aSOUBCF  /IELEMI 

S0R(NHI*S0R(NNI  »sou«cf  celemi 

ABCIBBImAMCOBIai, . 

AVC(PHI*AVG(>MI»1. 

A/fclMPHAMCKMIM. 

AVGINNMAVGtNNIM. 
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DIMENSION  SOI  (NNOOE)  ,PH!  : CNELEHI, ’Mill  IEIEHI,  PHIZ  (HE  LEM1 


DIMENSION  NOOE < A»N£LENI 


C0NM0N/ZZZ2/»SEX,NSa0DY.  MS, NT  (34) , I SF4  CE  (3L>, KNORNL  ( SHI . KMAKES  ( T4l 


OHHON/ZZZS/XPYCTR.YPVCTR.ZPYCTR.RXPYl.  ,R¥PYL,RZPTL 


C0NH0N/ZZZ7/YN0SE,YBD1,<SDZ,YTAI 


^TB  * kl 


UliInAim 


;0MN0N/Z7Z9/IYBD1,RZBD1,  RYBDZ.RZJDZ 


COMMON/ ZZZU/  XHU9CR,  YHU1  C R,  ZNUBCR,  R XHU  B,  TYHU3  ,RZHU3 


OMMON/ZZZt  VT  NET/ 5. THE  1C,  THETIS,  CON  I NS  , AZIMUTH 


r wm  tm  f TJ  4:1m  » 


OMNON/Z7ZH/¥SPAN,XL£ZV,XTEZY,TAMLE¥,T4NTE¥,TAU,7VTAU 


COMMON/ ZZZ17/NAKLPY,  NAK.  TB,  W4KLSK,  N AKi  SL 


C0HN0N/ZZZ19/KPRINT  ( Ml.  TREAD, NNRI  T E,  KRE40 

rnNNAy <¥¥¥«■/■¥ 


C0MM0N/ZZZZ1/ KPY1, RP¥? 


COMN  OM/ZZZZJ/KBOElE,  X3TT  MR,  KBOTY* 


CONMON/ZZZZVKPYflE.KPYSTP.KPYTY* 


comhon/zzzz7/ksmele,ksm!mp,kshty» 


common/zzzzvkskele,  ksk?tp,kskty» 


C0MM0N/ZZZ1|TNSTAU,N¥0RT , N SPIRAL, SPIH l 


MEAL  l shaft, i. shank 

no  II  > I«HM.|,MEH» — — 

PP-MOOE  ( 1,  IELEN) 

■— 

NP-MOOf  (J.IELENI 

— unnf  u.tei  m» 

PHIC  C IE  LEMI*  (SOR  IPPI  »SO-  (PHI *SOR  I MP  I *SOI  < MM  I I/A. 

bmH4  miAUi  sna(Rai*soi(RNfso«i«ai-soAiMNiwL. 

PHIP  I IELEHI*(  S0RIPPI  -SO--(PNI»SOR(RP|»SOR(MMII/N. 

BNlUIAHNI»4SARlRA|yI0i(BNlMt0R(»AJ*T0RIMNIt/A. 

rst  continue 

NPTTR  (Pf  YA  •» 

»l«  FORMAT  (//AX, •ElEN»,9X,  "MIC*, 111, 

l»PMf,  PPM  , YMMT  99  TPX  - PPMYt  9 1 

DO  ?ZB  IELE(*»1, HELEN 

Ui  HMtUJ UUHI  ■ pmi  Yirxi  m tmnni  f»> 

Ttl  FORMAT (5X.IB, 9X,BC19«5I 

mm 

END 
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1 HELEN) 


01  ME  MS  ION  MiLXt  HELEN)  .PHI  i CMELEM)  ,BHI?  (M  EL  EH) 

DIMENSION  PHI 3 ( NELEM) 

DIMENSION  HMmX (HELEN I.  Iff  » I CMELSMA.  

DIMENSION  PCC  J, HELEN)  , PI  ( 3.MELEH)  , PZCY, HELEN)  ,PJ(3,NELEN) 
CQMMOII/777t/<IXC3fcl  .MYfXLt  ■innHkl.»WI<«rif<»i.l»,w.li»t,ln,iMy 

CONNON/ZZZZ/NSPX,NS800T,  HS,NT(3%>  ,1  SPACE  ( SA)  , KNORNL  C SA)  . KNAKES  C3A) 

r.nnMn»i/777s/«aYi  ny-uaniwi  . mniYt.Mnim, mhtaii  r«..rt,.Un. 

t NBLEOE 

r.nHMOM/mh/i)H*rM.oHEaL  he  a.abeta 

C0NH0N/ZZZ9/XPYCTR,YPVCn,ZPYCTR,RX  PYL.RYPYl.RZPYL 

camin»i/77ZE/«Mnspf  «»m  t«%n»T«T»ti 

C0NH0N/ZZZ7/YN0SE,VB01,Y]0Z,YTAIl 

r.aMNnM/zzzE/ZMosi.ZBni.nni.mii 

CONNON/ZZZ9/AY8DI , RZBD1,  RY80Z,RZ9DZ 

C.OMMOMEZZZlBXBSMAEI,ltM>f  T, aSMAME. L SHANK 

COMMON/ ZZZt  1/  XHUSCR,  Y MU*  OR,  ZHUSCR, RXMUfl,  BYHII8 ,RZMU3 

CONN OM/miZy«a Dina,  BCNDRD,  TAUBL.ALEAB 

COHHON/ZZZ13/TMET79,TMCTIC,TMET1S,CON*NJ,AZIHUTM 

COHHOMXZZZtA/KaLAOE.  TEN.  IB ,TAMttE,XBti  ,« EYE .KBOTORS  1 YA) 

C0HM0N/ZZZ19/VSPAN,XL£ZV ,XTEZV,TANL EV, TANTEV,  TAU,  ZVT  AIL 

COMMONS  imt/MNEKBX,  MMAC  4»  .MHAKtK.HMA  It  EL 

C0MM0N/ZZZt7/N*KLPY,N«K.  4B,HAKISK,HAKLBL 

r.nimM/7Z71  A/HA  MOBY,  MAM.  i)B,MAMGSX.M  AM  i.BL 

COMMON/ ZZZ1 8/ XP RINK  10),  HREAO, NNRITE, KREAO 

COMMON  XIZZZBXBI 

COMMON/ZZZZl/KPYl.KPVZ 

aumatuinsa*  owu , hues  a a , kns  w 

COHMON/ZZZZI/KBOElE ■ KBD5  CP ,K8DTY* 

CnMNnM/mZN/KlMfl  F.KTMUR.ETMTTB 

COMMON/ 7ZZ2V XPVELE,  KPYS C P , rPYTY» 

COMHOH/7I77A<KMTELS,IUn<B.*MTT7a 

C0MH0N/ZZZZ7/KSHELE,KSN>  CP.KSHTYP 

COMM  CM/ 777 18/  KMBELl,JCMBIiiB  ,KM8  HB  

OOMMON/ZZZZ9/  XSXELE,  KSKC  CP , XSXTY8 

cnMMnMy777xii/KHiPi  r.xa.  ;-iP.rai  i» 

COMM 0N/ZZZ31/NSTA6, NXORT  .NSPIRAL, SPIRAL 


OOTPRO CXI, XZ, XS,Y1, V2, Yl)  «X1»Y1*XZ»YZ» XJ»Y3 

MErBBxcoi.ni,  m,Et„Ez,ts>  .mumhh 

VECPRYCOttOZ,  03 , El,  E 2,  ESI  «03*E1-01*E3 
XrrJBZIOUaZtB3,El,EZ,E3l»nttEZ-DltIl 

t* 

CM 

00  BOB  JELEMmI , HELEN 

C 

— AtXCPlIU,  JEAEME 
A1 YC«Pl  (Z,  JClZN) 

uzcrtimam 

AZXCaPZIl.JELEN) 

AZYCaBUZ.JELEM) 

* A2ZC*P2 (3, JCLENI 


A 1XA2X»VECPRX  (A1XC,A1YC*  A 1ZC, A2XC, A ZVC  ,A?ZC) 
-AlXA?X«t)ECPRX  1A1XC.  ACYC.  AC  7r.A7xr..A7VP  rA?rr.l 
A1XAZZ»VECPRZ  (AtXCfAlYC,  A t ZC  .AZXO, AZYP  , AZZC) 
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VCUTtJELtH»«PHIHJCLCMJ»«ltU*fHIHJEU£W  »A?TU 


bOO  CONTINUE 


SUBROUTINE  C®LINR(CP,VE.  X, BELT ,»ELZ  ,PHIC  ,NELEN,PC,<ROTOR» 


DIMENSION  VE.Y(NELEN>,VE.Z(NELEN» 


COMMON/ ZZZ1/TX (3%) ,NT(IW , NXY(J4» .KSTtMY, KSYNMZ, NSYNMV , NSYMMZ 


■ atm  4 rs  4 i.4 mil 


OMNON/ZZZ3/T  PYLON,  NBOW  l ,NBODYt,TBOOYI,imAIL,NSH*FY,NMUB,NSHANK, 


COMMON /ZZZ6/JM ACM,  OMEGA,  ILFA.ABETA 


«iT>F4ifi  i Aip* F-m f 


COMMON/  ZZZ6/XNO  SE,X801,<  102,  XTAIl 


Mi  III 


COMNON/ZZZ  A/Z  N0SE,ZB01,?30?,ZT4IL 


.T  ill  rmtJrL  Tlk 


OMMQM/ZZZH/RSNAFT.lSH.  rT ,RSHAN<,1 SHI N< 


HAIW  >i  Ail 


C0MM0N/ZZZ1Z/ RAOTOR, BCMD 1 0 , TAU8L, AL FA  3 


0MN0N/ZZZ14/K9LADE,  TAN.  EB , TANTE3,  X9L  E ,X3TE,KR0T0RS  ( 34) 


FFAiil/l  u jri 


COMMON/ ZZZ16/NMAKPY,NMA<-(B,NNAKS<,NHA<BL 


COMMON/ZZZlo/UANGPY,  MAN*-  RB , WANGSK,  M ANGBL 


COMMON/ ZZZ  20/ PI 


'Bfol.l.  Ul  V*M*A-lVA  MiM  * X*M 


COMMON/ ZZZ2Z/KNSELE,  KNS3RP  ,KNSTYP 


m A l 4 


COMMON/ZZZZ4/  KTNELc, KTNSHP, KTNTY* 


COMNON/ZZZ 36/ XV TELE,  KVnRP,KVTYY® 


COMMON/  ZZZZS/  KHBELE , KMBS  -tp  , KHB  TY® 


COHMON/ZZZ30/K3LELE,  K8LSRP ,K8LTf ® 


OOTPRO(*l,Yl,Zl,*Z,TZ,rt  ®X1®X2®YI®Y2*Z1®Z2 


XPROY  (OX,OYt3Z,EX,ET,eZi  *OZ*EX-EZ®DX 


8ETA2A»ABS<1.-UMACH®»2> 


UFREE»UMACM*111T 


OMCOAT»B. 


ONE6AZ>OME6AZ®2.®PI/6I, 


00  IBB  IELEM®!, HELEN 


RT«PC(?,IELEM) 


VX-VELXC IELENJ 


VZ*VEL  Z (IELEMI 


■Mj-.no-j 


i LA  .mi  k LA  kM.  rAJ  111! 
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